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PEEFACE. 



My object in the following treatise is to apply to the Kinetic 
^Theory of Gases a method of analysis different from that 
generally employed. It has been treated always on a certain 
fundamental assumption, namel y^ that i^hpi mnlppn^p of a gas 
are, as regards their relative motion, inde pendent of one 
another. As a consequence, we may say as the expression, 
of that independence, the law of distribution of momenta 
assumes the exponential form €~^^, and, so far as concerns 
translation velocities, 

ra being mass, and u, v, w component velocities. v/From this 
independence and from thi s form of Q ar e deduced Boltzmann's 
theorems, namely the H theorem, and that of thja-equality 
of mean kinetic energy for each degree of freedom. 

I propose to give to Q the more general form of a quadratic 
function, namely 

Q = 2w (w' + v^ + v^) + Xtb (uu' + tw' + mv'). 

Here 6 is a negative function of the distance r at the instant 
considered between the two molecules whose velocities are 
u, u\ etc., which/Tuuction is inappreciable except for very small 
values of r. ^I shall endeavour to prove in Chapters IV., V. 
that without the b coefficients the motion cannot be stationary. 
It has been proved abundantly that, assuming the independence, 
the motion is stationary with the usual form of Q, I question 
the axiom, not the demonstration. 

The con sequence of attributing to Q the p roposed new form 
is, that mole cules near to each other have on average a motion 
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i n the samft dj rfto±ion. They tend to form streams. That 
result, if it can be established, is worth investigation. 

For ordinary gases under ordinary conditions the b co- 
efficients are probably very small, and their eflfect negligible 
in such investigations as those of Tait and Boltzmann con- 
cerning diflfusion, viscosity, etc. But I think that the law 
a-hQ \y\ jtifl ql^-^roH form will express the state of the system 
without restriction as to density, except as follows. A physical 
limit there must be, when the gas liquefies under pressure, 
if not before. For it will not be contended that the dis- 
tribution of momenta among the molecules of the liquefied 
jJr gas is represented by the same exponential form as in the 
gaseous condition. An analytical lim it there is, when Q in 
its altered form ceases to be necessarily positive, that is when 
the determinant of the coefficients ceases to be positive. It 
"can be shown that this determinant does diminish as density 
increases, or temperature diminishes. But I have not cal- 
culated its value. It is therefore no more than a conjecture, 
^though perhaps a plausible conjecture, that the vanishing of 
the determinant may coincide with the physical change in 
the substance. 

It appears to me that the law of equality of mean kinetic 
I ^ energy for each degree of freedom cannot be reconciled with 
my proposed form of Q; that in fact the law holds only for 
the limiting case of a very rare gas. 

It is no light thing to question a conclusion maintained 
by^oltzmann, if indeed he does maintain this conclusion for 
^11 substances, or for all gases irrespective of density. I can 
but state the objections to this theorem, and to a certain aspect 
of the H theorem, as they appear to me. The reader will 
judge what weight is to be attributed to them. 

S. H. BURBURY, 

10 May, 1899. 
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CHAPTER I. 

OUTLINE OF THE THEORY. 

1. A GAS according to the Kinetic Theory consists of a 
great number of molecules in rapid motion. And the object of 
the theory is to explain on this hypothesis certain of the 
physical properties of gases. 

Any quantity of gas which can be isolated for the purpose 
of experiments is to be regarded as containing a number 
of molecules practically infinite. It is not possible to control 
or to observe the motions of individual molecules. But the 
theory assumes that such motion is subject to the usual 
dynamical laws. Also that if the gas, as an aggregate of 
molecules, be at rest, no dissipation of energy takes place. 

2. A molecule may consist of one or more than one atom 
according to the chemical constitution of the substance to 
which it belongs. It may be that hereafter we shall be able to 
explain on dynamical principles the chemical relations of atoms 
as constituent parts of a molecule, and of molecules inter se. 
And some progress has been made in this direction. At 
present the theory is concerned not with the chemical pro- 
perties, but with those properties of gases which may change 
without any change taking place in the chemical composition 
of the gas: for instance, density, pressure, and temperature. 
And as depending on these latter, it is concerned with the 
phenomena of viscosity, diffusion, and conduction of heat or 
electricity. 
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3. A molecule moving as such, one and indivisible, has 
three degrees of freedom in respect of its motion of translation 
in space. And it may conceivably have three other degrees of 
freedom, namely, motion of rotation about its principal axes. 
In respect of its internal constitution it may have many degrees 
of freedom, but we know nothing of the internal constitution 
of molecules except what chemists tell us. 

4. Any finite quantity of gas may, if protected from dis- 
turbance from without, retain for an indefinite time un- 
impaired its physical properties, for instance, its pressure and 
temperature, and also its chemical constitution. If therefore 
these properties are not changing, and are to be explained 
on dynamical principles as resulting from the motion of the 
molecules, such motion must be stationary motion. 



OF HOMOGENEITY. 

6. The gas, or system of molecules, may be in a certain 
sense, i.e. if we regard all the molecules at any instant within 
a finite space, homogeneous. For instance, a cubic centi- 
metre of air has sensibly the same properties in all respects 
whether taken from one part of a horizontal table or from 
another part. It is true that it will generally not have in all 
respects the same properties if taken from near the floor, as if 
taken from near the ceiling. It will, namely, be heavier, i,e. 
will contain more molecules, in the former than in the latter 
case. That is because the air in the room is in a field of 
uniform force, gravity. But this force will not necessarily 
affect any other property of the gas, at all events under 
ordinary conditions to which it is subject. 

If we assume the infinite divisibility of matter, or that any 
space, however infinitely small, contains an infinite number of 
molecules, the condition of homogeneity may apply to the 
molecules contained in any such infinitely small space. We 
may in that case say without ambiguity that the gas at a 
point P has the same properties as at another point Q. 



2; 



I.] OUTLINE OF THE THEORY. 3 

6. Modem physicists however teach us that the molecules 
have finite dimensions, or that there is a finite distance, how- 
ever small, within which two molecules cannot approach each 
other. That being the case, the homogeneity can be asserted 
only in the following form. If we consider ^ quantity of gas 
within a sphere of radius r described about P as centre, 
and a quantity within an equal sphere described about Q as 
centre, either sphere containing a number of molecules practi- 
cally infinite, then r being within certain limits, the firot 
quantity has the same physical properties as the second. 

OF THE DENSITY OF A GAS. 

7. The same doctrine of the .finite size of molecules causes 
coDsiderable difficulty in defining certain other properties of 
the system, for instance ihe density. 

The density of the gas is usually defined to be the sym of 
the masses of all the molecules in unit of volume. But if the 
molecules have finite dimensions, that may make the density 
at any point relative to that at other points depend on the unit 
arbitrarily chosen. It is not therefore a complete definition. 
We can however say, and it is generally sufficient for our 
purpose, 

(a) If there be N molecules within a definite volume 8, 
the average density at any point within 8 is NI8. 

(6) Again, we can define the density at a point which 
answers some general description, or as we may call it the relative 
density. For instance, suppose / to be a continuous function 
of the position of a point, then we may define the density 
with reference to /, as for instance at points where /= 0, 
provided that there exists an infinite number of such points. 
For consider any n of them, and about each a small sphere 
of volume <». Let n<o be the unit of volume. Let N be 
the number of molecules supposed of equal mass which, or the 
centres of inertia of which, are within some one or other of 
these small spheres. Then the density where /= is Njnw in 
the limit, as n is increased indefinitely, neo remaining constant. 

1—2 



4 KINETIC THEORY OF GASES. [CHAP. 

To define density at a point, as a function of the coordinates 
of that point, is with finite molecules a more difficult matter. 

Let r be a line of any arbitrary fixed length. As a first 
step to our definition, let us say the density at a point P, or 
Xyy, z ia the number of molecules, supposed of equal mass, con- 
tained within a sphere of radius r described about P as centre, 
divided by the volume of that sphere. If n be that number, p 
the density, then our definition is 

But that makes p generally a function of r as well as of 
X, y, z. How are we to get rid of r ? Boltzmann and other 
writers after or before him would complete the definition thus, 

p at P = J— ^ when r becomes infinitely small, any space how- 
ever infinitely small being assumed, for this purpose only, to 
contain an infinitely great number of molecules. This defi- 
nition does not materially di£fer from that which we should 
give if our gas, instead of consisting of discrete molecules, were 
a continuous fluid, any portion of which may be conceived as 
infinitely subdivided. 

But it is inconsistent with the teaching of these same 
writers, who maintain that molecules have finite dimensions, 
that in fact we can ascertain limits between which those dimen- 
sions must lie. It may be admitted that no ill consequences 
have yet arisen fix)m the use of this inconsistent definition, 
and if it be either necessary or useful, the inconsistency need 
not be regarded as a serious objection. Inconsistency may 
indeed be an advantage in mathematics as it often is in politics. 
It may be admitted also that this definition occurs to us by 
instinct when we are asked to define density at a point. That 
may be because we have passed through a former state of 
existence in which matter was supposed to be continuous, and 
therefore the definition strictly accurate. And our instinctive 
resort to it is due to unconscious memory of that former state. 

8. We may avoid the difficulty by making a further 
assumption with regard to the properties of our system. Let 
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To be the radius of a sphere described about P as centre which 
contains a very great number of molecules. Let rj be greater 
than ro in a finite ratio. Let as before n be the number of 
molecules within a sphere of radius r. Now let us assume, as 
a property of our system, that for all values of r between ro and 
ri, the ratio w: r* is constant. Then we may define p at P as 

r—^ using this constant value of n/r*. 

This method implies either that the state of our gas does 
not change from point to point in space, or that it changes 
so gradually as to be sensibly constant throughout the space 
between the spheres of radii r© and rj. 

9. Assuming that the state of our gas changes so rapidly 
in space that we cannot make the above assumption, and 
assuming that we scruple to base our definition on an assertion 
which we elsewhere contradict, it seems to me that the ex- 
pression "density at a point" or "number of molecules per 
unit of volume at a point " has of itself no meaning whatever 
when the molecules are of finite dimensions, but a meaning may 
be given it in an arbitrary way. One way would be to take an 
arbitrarily chosen volume for the unit. I would suggest the 
following as the best form of definition. Let/ be a continuous 
positive function of r which is equal to unity for all values of r 

less than a certain distance a, and for which -^ is always 
negative when r>a, and such that | r^fdr is a convergent 

J a 

series. Then define the density at P as Sm/, r being measured 
from P, m being the mass of a molecule and the summation 
including all molecules in space. For with molecules of finite 
dimensions, density at P must depend upon the distances of 
many molecules from P, but cannot be afiected by those whose 
distances from P are very great. 

THE STREAM VELOCITY. 

10. As with the density, so precisely the same difficulty 
presents ilisiglf if we try to define the str^eam velocity, or in 
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Professor Tait's language, mass motion, of the gas in given 
direction at any point P. M. Ladislas Natanson in a paper, 
to which I shall refer later, makes use of this function 
without giving any definition of it except by calling it 
" vitesse moyenne ou apparente." He is however, as I under- 
stand, willing to assume that any infinitely small space contains 
an infinite number of molecules. As the same difficulty arises, 
so the same solution mutatis mutandis may be resorted to, as in 
the case of density. 



OF THE PRESSUfiE OF A GAS. 

11. The motion which we attribute to our molecules cannot, 
assuming the gas to be at rest relatively to the surrounding 
space, be related to any particular direction in space. That is, 
at or near any point P there are on the average of time as 
many molecules moving with given velocity in any one direction 
as in any other. For if otherwise the pressure of the gas would, 
as shown by the next paragraph, be unequal in different 
directions, which is known not to be the case if the gas be at 
rest. 

12. Assuming the property mentioned in the last paragraph, 
and assuming the distribution of momentum among the mole- 
cules to be known, we can calculate the quantity of momentum 
transferred across an imaginaiy plane within the gas per unit 
of area and time. Let m be the mass of a molecule, p the 
number of molecules per unit of volume, to be defined as above 
explained, and pf(V) dV the number per unit of volume whose 
velocity of translation lies between V and V + dV. Then by 
our assumption the number of these last whose direction of 
motion makes with the normal to the plane angles between 
d and d + ddis 

y.f(V)dV8ia0d0. 

Their velocity normal to the plane is Fcos 0. Therefore the 
quantity of momentum normal to the plane transferred across 
the plane in one direction per unit of area and time by 
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molecules having velocity F . . . F+ dF is 



ir 
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iwpF» f 8in^cos»eci^ = im/)F». 

And the quantity of momentum so transferred by all the 
molecules is 

_^ " UNIVERSITY 

= ^p F« suppose, X^=,H^ 

F* being the mean value of F* for all the molecules. 

13. If further the mplecules exert on one another no forces, 
twice the quantity of momentum so transferred through unit of 
area of the plane is the normal pressure of the gas per unit of 
area of the plane. See Watson's Kinetic Theory of Gases, 
second edition, pp. 56 — 58. 

For we may suppose the gas confined in a vertical cylinder 
under a heavy moveable piston from which molecules striking 
rebound as elastic bodies. The momentum normal to the base 
of the piston, i.e. vertically upwards, which if it were an 
imaginary plape would be transferred through it, is in fact, 
the mass of the piston being very great, reversed, so that twice 
that quantity is transferred to the piston, and tends to make 
the piston rise in the cylinder. That is, the piston receives 
per unit of area and time from the impacts of the molecules 

momentum equal to ^mpV^ vertically upwards. It also acquires 
by the action of gravity per unit of area and time momentum 
Mg vertically downwards, where g is the force of gravity, and 
M bears to the whole mass of the piston the same ratio which 
unit of area bears to the base of the piston. For equilibrium 

^9 = ^P F* = Pi the pressure per unit of area. 

If the impacts of the molecules be suflSciently numerous 
per unit of time the piston will remain sensibly at rest. 
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OF THE TEMPERATURE OF A OAS. 

14. We have then in the absence of intermolecular forces 

p varies as p F". 

If the molecules do exert mutual forces on one another, the 
above value of the pressure is not accurate. But there are 
sufficient reasons for believing that for ordinary gases under 
ordinary conditions the pressure is not materially a£fected by 
the intermolecular forces. For such gases therefore, accurately 
or approximately, 

©Varies as p F*. 

But for such gases/accurately or approximately, the laws of 
Boyle and Charles hold good, whence p varies as p0f denoting 
the absolute temperature. It is thence inferred that, accurately 
or approximately, 6 varies as F". And that statement we 
must accept provisionally. Temperature measured from the 
absolute zero is represented by the mean square of the velocity 
of translation. On this theory temperature is a quality which 
cannot be attribiited to a molecule, but only to an aggregate of 
molecules. 

OF INTERMOLECULAR FORCES. 

15. It is however certain that for sufficiently dense gases 
intermolecular forces do exist, and we must now make a further 
assumption regarding them. We shall assume, namely, that the 
intermolecular force R between two molecules 

(1) Acts in the line joining their centres of inertia, or 
within very small limits of error may be regarded as so acting. 

(2) . Is a continuous function of the distance r between 
those centres of inertia, which becomes evanescent for all values 
of r exceeding a very small finite distance ro, which distance 
again is negligible compared with the dimensions of any mass 
of gas upon which we can make experiments. 

(3) The force must for sufficiently small values of r become 
repulsive, and should become infinite as r is indefinitely dimin- 
ished — but save for these restrictions it may be that the force 
is attractive for some values and repulsive for other values of r. 
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JB shall be taken as positive when repulsive. 

If c be the smallest radius beyond which the force exerted 
by a molecule at the centre on another molecule is negligible, 
c is defined to be the radius of action of the central molecule, 
and a sphere of radius c described about the molecule is its 
sphere of action, or its effective volume, 

THE FUNDAMENTAL ASSUMPTION. 

16. As above stated, the number of molecules of our gas 
per unit of volume may be diflferent at diflferent points of the 
system considered. Representing by p that number, we may 
have to treat /> as a function of x, y, z. At present it is con- 
venient to treat only the case in which p is constant throughout 
the system, and the distribution of velocities among the 
molecules is the same throughout the system. That being the 
case we may define as follows. 

At any instant let pfx{u)du be the number per unit of 
volume of molecules whose component velocities in the direc- 
tion taken for co lie between u and u + du. Similarly let 
pfy (v) dv be the number per unit of volume whose component 
velocities in direction y lie between v and v + dv, and pfg (w) dw 
the number whose component velocities in direction z lie 
between w and w + dw. 

We might say, and it would be precisely equivalent to what 
we have already said, that /«. (u)du is the chance that a given 
molecule shall have velocity in x between u and u + d% and 
fy (v) dv,fg {w) dw are the corresponding chances for y and z. 

17. These definitions, however apparently unambiguous, 
may be based upon either of two diametrically opposite 
assumptions. 

We may make namely, assumption A, 

A. The chance of any molecule having velocity in x between 
u and U'\-du is independent, not only of its position in space, 
but also of the velocities v, w, which it has in directions y and z, 
and further except in the case mentioned below, it is inde- 
pendent of the positions and velocities at the instant of all the 




10 KINFflC THEORY OF GASES. [CHAP. 

* 

other molecules of the system. The excepted case is when the 
two molecules are so placed that they are, or very recently have 
een, within one another's sphere of action. The force of this 
exception, and the necessity for it, will appear in the considera- 
tion of the H theorem. Chapter ill. 

On that assumption it follows that the chance of any given 
molecule having at any given instant velocities in a?, y, ;? 
respectively between u and u + duy v and v-^dVyW and w + dwy 

is fgn {^)fy {^)fz {^) dudvdw. 

Further, if we consider any pair of molecules, and denote 
by the suflSx 1 quantities relating to the first, by the suffix 2 
quantities relating to the second, then the chance that they 
shall at any instant respectively have velocities 

in CD between t^ and ^ + dt^, 

T^ and ii2 + dii2; 

in y between Vi and Vi + dvi, 

V2 and Vz + dv^; 

and in z between Wi a,ndwi + dwi, 

W2 and W2 + dw^t 

is fix (wi)/iy (vO/w (wi)f^ (t^a) . . . dvqdvidwidv^ . . . dw^^ 
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This or its equivalent is the assumption on which the kinetic 
theory of gases has hitherto generally been treated*. I shall 
reler to the state of things assumed a& j^ondition A . Without 
jfloubt it may be assumed legitimately with respect to the 
^limiting case of an infinitely rare gas, that is one in which 

* Dr Watson (Kinetic Theory of Gases, 2nd Edition, p. 3, equation B) 
assumes the independence of the chances for two spheres approaching collision. 
And that is sufficient for his purpose. Professor Tait (" Foundations of the 
Kinetic Theory of Gases," Transactions B, S, Edinburgh, 188S, Art. 21) makes 
an equivalent assumption. Dr Bo l tzmann (Vorlesungen Uber Gas Theorie, 
Part I, p. 23) assumes that Jhe_.motion_is^ and ever jiontinues to be, 
** molecular-ungeordnet . " I shall consider later (Chapter in.) what may be the 
effect of that assumption. In the meantime the weapon that Dr Boltzmann 
draws from that armoury is precisely the same as the one used by Watson and 
Tait, namely, that for all pairs of spheres or molecules approaching collision the 
chances are independent, as stated in my condition A. These writers do not 
define the exception, the occasion for definition not having arisen, but they do 
not assert the independence in the excepted case. 
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the dimensions of any molecule, or if it be a centre of force, 
its radius of action, are negligible compared with the mean 
distance between any molecule and its nearest neighbour. 
It is probable also that all the known gases under ordinary 
conditions of pressure and temperature, approximate closely to v^ 
the limiting condition of infinite rarity. And therefore the 
calculations made on the assumption of condition A or its 
equivalent by Tait, Boltzmann, and others, of the rate of 
di£fu8ion, conduction of heat or electricity, &c., for such gases 
remain unaffected. 

18. Or we may make assumption B, 

B. The chance of a given molecule having at any instant 

assigned velocities is not independent of the positions and 

velocities of all the other molecules at the instant. On this 

assumption B, instead of deducing the chance of the members 

of a group of n molecules having respectively at any instant the 

velocities 

1^... iLi + dvq^ &c. 

from the assumed chances for individual molecules, we must 
reverse the process. Let the chance that at the given instant 
the coordinates of the n molecules shall lie respectively between 

Xi and Xi + dxi, 
yi and yi + dyi, 

Zi and Zi-^dzi] 
for the first molecule, 

0^2 and x^ + dxi, 
&c. 
for the second, and so on, 

and their component velocities between 

ibi and 1^ + dTii, 

Vi and Vi + dvi, 

Wi and Wi + dwi] 



for the first molecule. 



i^ and v^ + dUi 
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for the second and so on, be denoted by 

Then the chance that the x velocity of the first molecule 

shall lie between 

Til and iii + dn^, 

whatever be the positions and velocities of the other n — 1 
molecules, is 

r+oo 

{dx^dx^ . . . dzndvidwiduz . . . dwn) F(x^i . . . t^j . . . Wn). 



J J J —00 



19. By way of illustration let us suppose, under condition A, 
where C is the usual constant, chosen so that 



of. 

J —00 



e-A«'dtt = l. 




Then on assumption A the chance that the members of the 
group of n molecules supposed to be of equal masses, shall 
respectively have their coordinates and velocities between the 
above limits is of the form 

C€~*^ctei • • • dZndUi , . . dWn , 

here Q = 2 (w* + v^ + w^), and h is constant ; 

and in dealing with assumption B, let us assume, the molecules 
being all of the same mass and the same structure, 

F{^i . . . 1^ . . . w») = Ce"*^, 
where Q=t^^ + V + Wi' + t^*+ ... +w^»; 

+ 6ia ( i*i^ + W2 + WiWa) 

+ h^ (^^8 + Wz + WiWg) + &c. 

Here the h coefficients shall be functions of the distance r, at 
the instant between the pair of molecules to which the suffixes 
relate, which functions become evanescent when r exceeds a 
certain very small distance. In that case omitting the constant 

rrr+oo 
f(ui)=^ III €''^^dXi,,.dZndVidWidU2*''dWn 

= € 2i>u , as shown in the Appendix (g). 
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Here D is the determinant of n^ constituents 



2 


6m 


hu .. 


i« 


2 


&S8 .. 


• • • 


K 


2 



2) = 



and Du is its first coaxial minor. 

The form o( f{u) under assumption B is the same, except 
for the factor of A, as under assumption A/But the physical 
state of the system may be very different. /For instance, under 
assumption B, two molecules at a distance r apart, so small 
that the corresponding h is not negligrihle^ ar^^ if h hft nftgativft^ 
on average movinpf in the same direction . According to Boltz- -^ 
mann (Vorlesungen, p. 21) the motion in this case is not 
" molecular ungeordnet." 

20. Of these two assumptions A and B, B includes A as a 
particular case. It is therefore better, if and in so far as it 
may turn out to be possible, to establish a proposition on 
as sumption B, ra ther than on assumption A. My object in 
the present treatise is to show that it is generally possible. It 
will appear also that A requires us to assume our gas to be 
infinitely rare, whereas B imposes upon us no stringent con- 
dition with regard to the density. The only condition in fact 
which appears to be required at present is this : the coefficients 
b must be such that for all possible values of the u*s, v% and 
w*s, Q shall be positive. Mathematically that is expressed by 
saying that the determinant D and all its coaxial minors must 
be positive. 

D will appear to be generally a function of density and 
temperature, and if ever D = the mathematical treatment 
changes, whether or not the physical system changes its state, 
and if so whatever the nature of that change may be. 

It will appear also that if u, v, w be the component velocities 
of a molecule of mass m, u\ t/, w' those of any other molecule, 

^ (W + iw' + tvw') 

is related to the Virial defined in Chapter ii. of the intermolecular 
forces acting on m. See Arts. 58, 85. 



CHAPTER II. 



CLAUSIUS' THEOREM. 



21. It will be necessary to consider at length later the law 
of distribution of velocities or momenta among the molecules 
for any given value of the mean kinetic energy — that is the 
form of the function which we denoted above, Art. 12, by /(F). 
The following proposition, due originally to Clausius, is inde- 
pendent of that law. 

If m be the mass of any one of a system of molecules in 
stationary motion, a?, y, z the coordinates of its centre of 
inertia at any instant, we have, ^ denoting summation for all 
the molecules of the system, and Smo;, Smy, Xttiz being zero 

or 



2m(^)' + 27n..^ = 0, 



with similar equations for y and z. 
It follows that 



X, F, Z being the component forces, whether external or 
intennolecular, acting on a molecule. 
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22. We will now however reserve X, Y, Z for the com- 
ponents of external force only, and let R be the intermolecular 
force between two molecules acting in the line joining their 
centres, and let R^, Ryy Rz be its components. So that Clausius' 
equation becomes 

-S2(i2^+i2yy + J2^) (1). 

Half the expression on the right-hand side of the last 
equation is called the Virial. And separately 

-iiXxA-Yy + Zz) 

is the Virial of the external forces, and — J22(-Ba^«? + Ryj/ + Rzz) 
is the Virial of the intermolecular forces. 

23. Let one molecule be at P or x, y, z, another at P' or 
of, y\ /. Let r^ = (a? - i»')2 -h (y - y7 + (5 - 2!)\ Then for the 
force on P due to P, R being taken as positive when repulsive, 



J2« = 



"« 



a? — a?'^ ^_y~3/-- ^ z — ^ 



Ry Ry=^ 



Rf Rg = R. 



T " If* ~ ip 

Similarly for the force on P due to P, 

J. x-af p p y-y' p p z-z' 



R. 



and therefore for two molecules at P and P' 

And so the term in Clausius' equation due to the inter- 
molecular forces is 



.(2), 



22 (Ra^x + Ryy + R^z) = 22 J2r 
each pair of molecules being counted once. 



24. We have next to consider the external forces, whose 
components are X, Y, Z. I shall follow Van der Waals' 
method. 

Let the whole system be enclosed by a surface 8, through 
which the molecules cannot pass. Let p be the normal pressure 
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per unit of area on S, directed inwatds. At any point P on 8 
let X, fi, V be the direction cosines of the normal « to /S at P, 
I, m, n those of a the line joining P with the origin. Then for 
a molecule at P 

— {Xx + Fy + Zz) =p {l\a + mfia + nva) 

=jpacos(«, a) 
da 

if ds be an element of the normal, measured outwards, 



= 4p 



ds ' 



Let U8 now assume p to be constant over the whole surface 
S. Then 



-jj(Xw+Yy + Zz)d8 

da^ 



-i.ii- 



, dS 
ds 



= iP III ^*a^dxdydz 

throughout the space enclosed by S, by Green's Theorem, 

= Spa>, 
if © be the volume enclosed by 8, 

We see then that the Virial 

orby(l) f 1)0) = 2im (i^« + v* + w») + iSSiJr (3); 

or, if T be the mean kinetic energy of the system of molecules, 

T 

and — represents the mean kinetic energy of a molecule. 

r 

If we take the volume co for a generalised coordinate, 
supposing it to change in magnitude but not in shape, the 



mmma^mBmmmmmmmmimmmmmmmmmmmmmmf^mm^a 
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other coordinates of the system may be a;©*, yar, &c., where 

T dT 
X, y, &c. are numerical. Then | - = — , and if x ^® *^® 

potential of the intermolecular forces ^SXRr/(o = — -^ , and our 

equation assumes the Lagrangian form p — '-j -r^ , p being 

the applied force which maintains © constant. 

25. If the surface 8 be an imaginary closed surface, the 
system of molecules in the external space being homogeneous 
with the system within 8, the pressure per unit of area of 8 
is given by the above equation, 

^pto = Siw (u^ + v^ + w^) + ^tlRr, 

The whole Virial or fp® — jS2-Kr is thus equal to the 
whole mean kinetic energy of the molecules within 8, It 
does not follow that the two parts of the Virial, namely ^po) 
and — 22 J-Kr, are separately proportional to 2 Jm (u^ + v* + w^). 

26. If 8 be the external elastic boundary of the system 
through which molecules cannot pass, we might calculate p by 
the same method as we calculated it in the case of there being 

no forces, namely p = ^mp'V^^, where p' is the mean density, and 

^V^^ the mean kinetic energy, of molecules near 8, This 
result is not inconsistent with Art. 24, because the mean 
potential energy near 8 is, in the case now supposed, different 
from the mean potential energy throughout the enclosed space, 
T^ and soT'^ 4= V^ and />' + p. 

The value of 22-Br /or Elastic 8pheres. 

27. We have assumed i2 to be a finite force. In the 
limiting case of so-called elastic bodies, the force which acts in 
collision between two such bodies is, according to the usual 
convention, an infinite force which reverses the velocities of the 
colliding bodies normal to their tangent plane at collision in an 
infinitely short time. As a consequence of the time being 
infinitely short, no change takes place in the coordinates during 

B. 2 




* 
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collision, and no body can ever be in collision with more than 
one other body at the same instant. 

It is of course admitted that no such body exists in nature. 
Bidt the behaviour of gas molecules in regard to their mutual 
yaction is probably, at least in the case of rare gases, approxi- 
y mately the same as it would be if they were elastic spheres of 
this form. And the assumption made in the case of the 
sph eres that the kinetic energy of translation of two spheres in 
collision is in the aggregate unchanged by the collision, must 
if molecules are stable be true on avfi rflgft fnr any fo rm of 
molecule. 

28. It is then possible to calculate X^Rr for a system of 
elastic spheres, regarding the force as finite, and suflBciently 
great to produce the reversal of the normal velocities in a 
finite time t, and treating the coordinates as constant, to proceed 
to the limit when the force becomes infinite, and r zero. 

Let the colliding spheres be each of mass m, and of 
diameter c. Let p be the number of spheres per unit of 
volume. 

Let T be now the mean kinetic energy of a sphere. And 
let u, V, w, u\ t/, v/, be the component velocities of any two 
spheres, V their relative velocity. Then 

and mV^ = 2m (w* + 1;* + v^) on average 

= 47. 

Let be the angle between the relative velocity and the 
line of centres at collision. Then for each sphere the normal 
component of relative velocity is JFcos 0, and it is reversed in 
time T. The assumed force is then J? = mF cos 0/t, and it acts 
during the time r. The number of collisions per unit of 
volume and time is, given V, irc^Vp, And the number for 
which lies between and tf + d^ is 

2ir(^Vp^\n0cos0d0. 
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Also r = c. Hence for all collisions of a pair of spheres with 
given V 

XRr = 27rm(f^V^p f^sin cos^ 0d0 

Jo 

And for all collisions of a pair of spheres with whatever 
value of V, 

XRr = ^ird'pmV^ 

or ^tEr = f TTC^pT. 

But each collision is now counted twice, namely, once from 
the point of view of each of the spheres concerned. Counting 
each collision once, or from the point of view of one sphere only, 

i2iJr = |7rcVr (4), 

whence we obtain for each sphere 

^ iSiir = §7rc»r, 
and for all spheres 

and therefore ^ 2S-Br = if 7rc*/> . pT. 

I shall write k for ^irc^p. If a be the radius of a sphere 
which on average contains one molecule, evidently 

or K is four times the aggregate volume of all the spheres in 
unit of volume. 

It follows from Art. 24 that 

f 2)6) = 2im {u^ + 1;= + vj') + ^^kT, 
or p = i(l + ^K)pT (5). 

On a system of Elastic Spheres in Vertical Golmnn. 

29. The following investigation depends upon the usual 
assumption made with respect to collisions between elastic 

2—2 
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bodies, that a certain inomeDtiim is instantaneously transferred 
through a finite space. The responsibility for that assumption 
rests on the inventor, whoever he was, of the theory of rigid 
elastic bodies. Consider a vertical column of gas, or material 
system, whose molecules are equal elastic spheres, under a 
constant vertical force/. 

If p be the pressure per unit of surface, x the height of a 
point in the column above a fixed plane, m the mass of a sphere, 
p the density, p^ the density at the base, we have 

dp _ 



'^fp> 



also by (5) p = § (l + ^k) pT. 

Whence, if we make 



we find 



or 



T = constant = 77-, 



dx 



dx 



Again, consider N spheres crossing the plane a? = with u 
for vertical component of velocity. Of these some, sayiV— JV', 
will reach the plane x=dx without collision. N' will undergo 
collision before reaching dx. But for these N' there will in 
stationary motion be substituted, as the result of collisions, iV" 
other spheres with the same vertical component u. 

Now, if the impact were direct, i.e,, the line of centres at 

collision vertical, the eflfect of the 
collision would be to substitute for a 
molecule with vertical velocity u at 
height X a molecule with the same 
vertical velocity u at height x + c, 
where c is the diameter, or, as we 
may express it, the substituted sphere 
would gain a vertical height c, with- 
out losing in respect of that distance 
any kinetic energy to the force/. We might without aflfecting 
the general motion suppose the two spheres to change places 
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after collision, so that one would lose, and the other gain, the 
height c without change of kinetic energy. If, therefore, all the 
N' collisions were direct, the average height of the N spheres 

doc 
or their substitutes at the end of the time — would be, not dx, but 

u 

dx+ j^c, while their loss of kinetic energy would be Nmfdx. 

But all impacts will not be direct ; we must consider then 
the result of indirect impacts. For this purpose consider two 
classes of collisions, in one of which the sphere A has vertical 
component u before collision, and in the other A' has vertical 
component u after collision. The eflfect of a pair of collisions one 
from each class, is to substitute A' for A as the sphere with verti- 
cal component u. Now let I (see figure) denote the vector line 
of centres at collision, and cos (ul) the cosine of the angle which 
I makes with the vertical. Then in the first of the pair of 
collisions the centre of A is below the point of contact by 
\c cos {vl). In the second, the centre of A' is above the point 
of contact by ^c cos (wZ). There is no reason why the point of 
contact should be higher or lower in one case than in the other. 
It will be on average at the same height. Therefore on average 
of all pairs of collisions substituting A' for A with vertical 
velocity u, A' is above A by 



c cos {ul) = r, suppose. 

Let q be the relative velocity of the two colliding spheres. 
Then considering q as radius vector of a cone of axis w, and I as 
radius vector of a cone of axis q, we get 



r = c cos (uq) cos {ql) 



= |c cos iuq)y 

because cos {ql)= | cos^ 6 sin 6 dd P cos ^ sin 6 dd = ^. 

Let V be the absolute velocity of the sphere whose vertical 
component of velocity is u, so that 

cos (uv) = - . 



22 KINETIC THEORY OF GASES. [CHAP. 



Then ^ = f c cos (uq) 



u 



= |c - cos iyq). 

Let -^ be the velocity of the other colliding sphere, E the 

angle between v and ^, Then' 




« 1* v — ilr cos E 
^=§^- ^I 



V 

We have to multiply this by the 
number of collisions which N 
spheres having velocity v undergo with spheres of velocity 
'^ .,.'^-\- dyjr, making with v angles E .,,E + dE in time dt, or 

QJJC 

— , and then integrate for all values of -^ and E, 

Let pf{'^)d'^ be the number of spheres in unit volume 
with velocity '^...'^ + (i^. The result is 

Nird'p j^ dytrfiir) J^ ^ sm EdEq^c ^ — 

= Ittc'/) . Ndx 
= icNdx, 

Therefore at time dt the average height of the N spheres or 
their successors above the plane izr— is (1 + k) dx. 

But the energy which allowing for substitutions they lose 
in the ascent is Nmfdx, The loss takes place only during 
free path. It follows that the loss of energy due to the ascent 
dx is, allowing for-substitutions, m/ctz?(l—ic) per sphere. 

Now suppose that at 4? = the number per unit of volume 
of spheres having Jj^mt^^ . . . d {v?) for energy of vertical velocity 
is 

^g-Amua^(^2) ^ (^)^ 

where K is constant. Then, by what has been proved above, 
the number which at height 9a? have i^vJ^ .,, d{y?) for energy 
of vertical velocity is 
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and the number which at height 'dx have 
for energy of vertical velocity is 

g— afcn%fa»g— 9ie ^g— A (mtt«— am^pa?) (i-ic)) ^|^8 

that is, K€-^'^*du' (6), 

because, neglecting ic^ &c., 

and therefore ^hmficdx = — 3ic. 

The two groups (a) and (6) are equally numerous, and 
therefore either can by ascending or descending, allowing for 
substitutions, exactly replace the other. Now this is the 
reasoning by which in the ordinary case, when /c = 0, we prove 
p = pq€-^'^^. It now proves p = po€"^'^*€-*. 

I shall refer to this problem later for illustration. 

30. The pressure per unit of surface is increased in the 
ratio 1 : (1 + ^k) as the molecules, from being material points, 
become spheres with finite diameter o. 

But the pressure per unit of surface is the quantity of 
momentum which is carried through unit of surface in unit of 
time. Art. 12. Now, so far as this momentum is carried 
through the surface by molecules during their free path, it is 
not altered in the least by k acquiring finite value. The 
increase of the transfer of momentum consists in the case of 
elastic bodies in the process above explained, namely, the 
instantaneous transfer of momentum through a finite distance 
which occurs on collision. 

The result of Art. 29 may also be explained thus. The 
mean effect of all collisions on any sphere is equivalent to a 
force k/ acting from the greater to the less density. 



CHAPTER III. 

THE T>TRTRTP7TXTQ|<[ Qf MOMTTArTA A MONa TRTC MOLECULES, 

ASSUMING CONDITION A. 

31. The simplest form of molecule that we can conceive is 
an elastic sphere, such as described in Arts. 28, 29. 

An elastic sphere, so conceived, may be considered as having 

three degrees of freedom only, namely, motion of translation 

in space. For ^the spheres being supposed perfectly smoot h. 

mot ion of rotation about an axis, if non-existent, will not be 

V^ produced, and if existing will not be altered, by collisions. 

Let Uy V, w be the component velocities of a molecule, m its 
mass. Then the molecules being elastic spheres of mass m, let 
the number per unit of volume of molecules whose velocities 
in direction a? lie between u and u + du be proportional to 
g-*wua^^ This is called Maxwell's distribution.vM[t has been 
proved by several writers that assuming this distribution 
to exist, and assuming condition A (which is necessarv)>^the 
difitribution i^ffl riftf b^ Higf.nrKftrl hy r^nllfqinng bctweeu the 
spheres. 

My object in the present chapter is not to give another 
proof of this proposition, but to show in what manner condition 
A affects the truth, as well as the proof, of the proposition. I 
shall therefore follow the general method elaborated byJBoltz- 
mann and Watso n, as being more convenient for my purpose 
than other known proofs, e,g. Professor Tait's. 

32. Consider two sets of spheres, one set having mass M, 
the other mass m. Let the number per unit of volume of M 
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spheres be p^, of m spheres p^. Let us form a ^agram of 
velocities, the velocity of any sphere of either set being repre- 
sented in magnitude and direction by a line drawn on the 
. diagram from the origin to any point P or x, y, z. At P 
suppose an element of volume dxdydz. The number per unit 
of volume of spheres M whose velocities are represented by 
lines drawn from to points within that element shall be 

pM^i^y y> z) dxdydz, 
or shortly p^Fdadydz. 

Call these the class F (dxdydz). 

Similarly the number per unit of volume of spheres m whose 
velocities are represented by lines drawn from to points 
within the element of volume d^d/qd^ at the point f, 17, (;" shall 

be pmf{S>V,OdSdvd^, 

or shortly pmfd^drfd^, 

and we will call these the class 

fd^dr)d^. 

We may also express the above definitions as follows. The 
Amance that a sphere M shall belong at any instant to the class 
^ Fdxdydz is Fdxdydz. The chance that a sphere m shall belong 
at any instant to the class /dfdi^df isfd^di^d^. 

33. Now assume condition A. 

Then the two chances are independent. And therefore the 
chance that a pair of spheres, M and w, shall belong, M to the 
class Fdxdydz, and m to the class /dfdiydf, is 

Ffd^drjd^dxdydz. 

On our assumption of condition A this is true even if the 
two spheres are on the point of collision, and therefore very near 
each other. For that case is not within the exception defined 
in Art. 17. If we assume B instead of A, it may be true only 
if the two spheres are at a considerable distance from each 
other. But in this chapter I assume condition A. On this 
assumption the number per unit of volume of pairs of spheres, 
M and m, which have velocities between the limits aforesaid, 



26 KINETIC THEOBT OF OASES. [CHAP. 

and for which m lies within the element of volume dto^ so situated 
that M and m are approaching collisiony is Ffd^cUfd^dadydzdo). 
If the element dm be so situated that M and m are separoiting 
after collision^ the case is within the exception of Art. 17, and 
therefore condition A does not apply. Let us now define dto. 
Let, namely, i2' = (a? - f / + (y - vf -^(z-^, so that R is the 
relative velocity of M and m. 

Let c be the sum of their radii. About the centre of M 
suppose a circular area described with radius c in a plane 
perpendicular to R, and let 

d(o = TTC^Rdt 

Then the number per unit of volume and time of collisions 
between spheres M of the class Fdxdydz, and spheres m of the 
class /dfdiydf, is 

PMpmFdadydzfd^Tid^wc^R, 

And the whole number per unit of time of collisions which 
members of the class Fdxdydz undergo with spheres m of any 
class is 

pMPmfdxdydz 1 1 1 fd^'nd^Trt^R. 

Each of these diminishes the number of the class Fdxdydz 
by one. So the rate of diminution per unit of time of the class 
Fdxdydz due to collisions between its members and spheres m 



is 



PMpmFcLcdydz jjj fd^dvd^TT&'R. (6). 



The Effect of Collisions. 

34. We have now to find the number of spheres M which 
by collisions with spheres m pass into the class Fdxdydz per 
unit ofvolume and time. Let P,phe the centres of M and m. 

On a collision the only thing that changes is the direction 
of i2. And that change depends upon (1) the angle i/r which 
the line of centres Pp at the instant of collision makes with R, 
(2) the angle ^ which the plane through a line drawn through 
the centre of M parallel to R, say Pn, and the point of contact 
of M and m makes with a fixed plane through Pn, 
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The direction of R after collision is in the same plane 
through Pn, and makes the angle 2^ with Pn, That is the 
easiest way to picture to ourselves the effect of a collision. 

In order that a collision may take place, it is necessary that 
a line ps drawn through p parallel to R should cut the circular 
area above mentioned. We may divide that circular area into 
elements whose type is 2c" sin -^ cos -^d^d^. Assuming con- 
dition A, the line ps is as likely to cut the circular area in any 
one point as in any other. Therefore the chance that the new 
direction of R shall make the angle 2^...2'^ + d^ with the 
original direction of R is proportional to 2 sin -^ cos -^d-^, that 
is to sin 2yjrd'>lr. It follows that, assuming condition A, for any 
given direction of R before collision, all directions after collision 
are equally probable, and on average the whole kinetic energy 
is divided equally between the spheres. 

36. The direction of R after collision is then a function of 

^, y> Sly f , % ?, -f , ^. 

Let the new values of these variables be denoted by the 
corresponding accented letters. Evidently '^, ^ are uoaltered, 

or ^' = ^, yjr^^ylr. 

All those pairs of spheres for which before collision the 
variables a?... fare between 

and '^y between limits 

-^ . . . -^ + d'^y ^ . . . ^ + d(f> 

will after collision have these variables between limits denoted 
by the corresponding accented letters. 

Now let pj^'dx'dxfdsi be the number per unit of volume of 
M spheres which at any instant belong to the class F'daldifdz', 
Similarly after collision m passes into the class f'dl^dy]'d\^y 
whose numbers ^^f'dl^d/irld\\ 

36. If after collision the velocities of M and m were 
reversed, without change of their positions, they would by a 
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reverse collision pass into the classes Fdxdydz and fd^dtfd^ 
respectively but with reversed velocities. 

It follows that the number of pairs of spheres, M and m, 
which before collision belong to the classes F'dafdj/d/ and 
y'df di/'df respectively, and by collision with given values of 
y^ and ^ pass into the classes Fdxdydz and fd^dr)d^ respectively, 
is, again assuming condition A, 

PMpmF'da)'di/d/fd^'dvd^ 2 sin i^ cos ^dy^djx^R. 

37. Now afy y', /, f ', rf, f ' are functions of a?, y, z, f , 17, ?, -^j 
^, and of those quantities only. Also the system of two spheres 
passes with unchanged kinetic energy from one state to the 
other. 



^ 



Whence it follows by the general proposition proved later. 
Art. 50, and can be proved independently, see Watson's Kinetic 
Theory of Oases, second edition, p. 9, that 



dxdyd2!dl^dr)'dK' = dwdydzd^vd^ (7). 



Therefore the number of pairs of spheres, M and m, which 
per unit of volume and time pass out of the classes F'doddffdsi 
and f'd^dAf{di,' into the classes Fdxdydz and fdl^difd^ respec- 
tively is with given -^ and ^ 

PMpmF' fdxdydz d^drid^ 2 sin i/r cos yltdylrd<f>c^K 

And if now f, 17, f, •^, ^ vary, we find for the number of 
M spheres by which, owing to collisions between M and m, the 
class Fdxdydz is increased per unit of time the expression 

p^p^dxdydz \ \ I d^drjd^ P / 2 sin -i/r cos '^'^d^F'/'c^R 

J J J -<x> J oJ 

= pMpmdxdydz j I j d^d'qd^F'fTrd'R. 

But the number by which the same class is diminished per 
unit of time by collisions with the m's we found to be 

pMPmFdxdydz 1 1 j fd^drid^irc^R, 
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If therefore -tt dasdydz denote the change per unit of time 
of Fdxdydz, due to collisions of M with m, 

^ dasdydz = dasdydz jjj d^dvd^ (F'f - Ff) wc'R . . .<8). 

Similarly 

^df di7d?= d^dvd^ jjj dxdydz {F'f - Ff) ird'R (8a). 

38. If Ff = 75^ in all cases in which a pair of spheres M 
and m can by collision pass from the classes F and / to the 

classes F' and/' respectively, or vice versa, ;ir = ^* Similarly ^^ 

4^ = 0. And this being true for all values of F and /, the 

motion is stationary. Assuming then condition A, F'f = ^is 
a suflScient condition for stationary motion. 

The proof that it is a necessary condition is due to Boltz- -^ 
mann. and is as follo ws. 



Let H= III " dxdydzF (log ^ - 1) 

+jjf_y^dr,d^/aogf-l). 



Then f=///eWyd.flogi' 

+ fffd^d^^logf. 

JET Jf 

And substUuting for -tt and -^ their values above foirnd, 

^ = jjj^'^ dadydz jjfy^vdifji^'f " ^f) (log ^ + Hf) ■^<^R 
= ijj dxdydz jjj d&lvd^ {F'f - Ff) log (Ff) ■k<?R. . .(9). 
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But in the course of this integration x, y, z and f , 17, f assume 
all possible positions on the diagram of velocities, and therefore 
assume the positions corresponding to any F' and/', and there- 
fore F and / among other values assume the values F' and /'. 
Whence 

^ = j^l dadydz jjjd^dnd^ (Ff - F'f) log F'f'^d'B. . .(10). 

and adding together (9) and (10) 

^ = i JJI dxdydzjjjd^vd^F'f' - Ff) log ^ irc^R. . .(11). 



/ 



And this is necessarily negative, if not zero, and then only 
zero when F'f' = Ff for every case in which two spheres can 
pass by collision from the classes F and / to the classes F' and 

JIT 

/' respectively, or vice versa. Bu t in stationary motio n -^ 

must be zero . Therefore in stationary motion, given condition 
A, Ff = Ff, is a necessary condition. A distribution of 
velocities in which F'f = Ff for all cases in which a pair of 
spheres can pass by collision from the state F and / to the 
jj/ state F' and/' or vice versa shall be called the normal distri- 
^f\ hution. Any other distribution shall be called an abnormal 
distribution. The H theorem proves that when the distribution 

JIT 

is abnormal — ^ is negative. 

This is known as Boltzmann's minimum theorem, or the 
H theorem, and IT is Boltzmann's minimum fimction, or, as he 
-^ sometimes calls it, the Entropy function . The theorem is as 
here given, founded on, and depends for its existence upon, 
condition A, without which we cannot use the product Ff to 
express the frequency, or chance, of the simultaneous occun^ence 
in collisions of the two states denoted by /and i^. 

On the H theorem. An objection considered, 

39. If when the entire system of elastic spheres has, with 
H diminishing, reached a certain state, all the velocities were 
simultaneously reversed without change of the coordinates, the 
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system would exactly retrace its course with H increasing. It 
is therefore possible, when a system is set in motion, for H to 
inprease in that motion. This was at one time considered an 
fi(bjection to the theorem, because the theorem was supposed to 
prove that vmder no possible circumstances could H increase. 
It however does not profess to prove that H always diminishes, 
but that, assuming a certain condition H must diminish on 
a verage. For the assumption underlies the whole proof as 
hitherto given. 

Now to meet this objection, consider two spheres, M 
having velocities oc, y, z^ and m having f, 17, f. Construct 
a spherical surface of radius c about the centre of if, and 
consider the element of volume 

2c* sin ^ cos-^ d-^ d^ Rdt = dto. 




If m be within d©, the two spheres are either (1) approaching 
collision, or (2) separating after collision, according as dm is on 
one or the other side of M. Condition A asserts that the 
number per unit of volume of pairs of spheres, M and m, 
belonging respectively to the classes F and /, which are op- 
proaching a collision of that kind— i.e. a collision with the given 
yfr and ^ — is Ffdx...d^d(o, see Art. 17 and note. But con- 
dition A does not assert that the number per unit of volume of 
such pairs which are separating after a collision of that kind is 
Ffdx,..d^da), because that case is within the exception of 
Art 17. In fact if x, y, z and f , 17, f denote the precollision 
velocities of the separating pairs, x\ yf, &c. their velocities after 
collision, the separating pairs are the identical pairs which dt 
seconds ago had respectively velocities 

a?. ..5?+ dx f. ..f + df. 



i 



* 
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and were approaching a collision of the kind in question. The 
number per unit of volume of the separating pairs is therefore 
by our assumption of condition A 

or which is the same thing, 

Ffdaf, . . .d^dco, 

and not Fy'dx...d^'da). (We are assuming at present that 
Ff^Pf',) Now in the reversed course the pairs which in 
the original course were separating have become the approach- 
ing paii*s, and vice versa. Therefore in the reversed course the 
number per unit of volume of pairs, M and m, which, belonging 
respectively to the classes F' and/', are approaching a collision 
of the kind in question, is ^/da?'...df'da), and not F'f'da!,.. 

/rff'dtt), as according to condition A it should be. We see then 
that condition A is not satisfied in the reversed motion, and 
why not. 

Thft fact then that in the reversed motion B. increases with 
the time is no objection to the H theorem as a mathematical 
proposition, bec ause the reversed motion does not satisfy the 
condition on which the H theorem is based. 

40. We now see the necessity for the exception referred to 
in Art. 17. For if condition A continues to exist throughout 
the direct motion without that exception, it must also continue 
to exist throughout the reversed motion. Therefore by the H 

theorem -=- is negative or zero in both motions. But if it be 

negative, and differ from zero, in one of those motions, it must 

be positive in the other. Therefore "jr must, if we assume 

condition A without exception, be zero in both motions. 

In other words, if we assume condition A with the ex- 
ception, the mathematical consequence of our assumption is 
that ^is either minimum or diminishing, as the theorem asserts, 
f we assume condition A without the exception, the mathe- 
matical consequence of our assumption is that H is minimum 
already. This condition is satisfied by Maxwell's distribution. 




III.] THE DISTRIBUTION OF MOMENTA, ETC. 33 • 

41. We must now consider Boltzmann's own assumption 
that the motion is " molecular-ungeordnet ." It being assumed 
that in the direct course the motion is molecular-ungeordnet, is 
it molecular-ungeordnet in the reverse course, or not ? I think 
Boltzmann's answer to this question would be in the negative, 
yOn the ground that we, or Maxwell's corps of demons, have 
" ordered " (geordnet) the motion by the very act of reversal. 
If this be so, then " molecular-ungeordnet " has, as applied to 
this theorem, precisely the same properties as my condition A. 
See Vorlesunpen, p. 42. 

Let us endeavour to construct synthetically a system which 
shall without doubt be molecular-ungeordne t. The molecules 
being distinguished by numbers, I ask (say) Dr Watson to 
assign velocities to them according to any law he pleases. 
Then I, in complete ignorance of those assigned velocities, 
scatter the molecules at haphazard through space, and they 
shall start from the positions which I so give them with the 
velocities so assigned by Dr Watson. That is, primd facie, a ^ 
molecular-ungeordnet system — in fact it is as near an approach"^ 
to chaos as is possible in an imperfect world. 

Clearly also in this system condition A is satisfied, and the . 
system will, or in all probability will, move from its initial/ 
position with H diminishing. Equally clear is it that if in the 
initial position all the velocities were reversed, condition A 
would be satisfied, and the system would, or in all probability 
would, move oflf in the opposite direction with H diminishing*. 
It follows then that if the selected distribution of velocities be 
abnormal, the state in which we have placed the system is one 
from which the system can and probably will move with H 
diminishing, but into which it must move with H increasing. 

dH 
* This statement is not inconsistent with that of Arts. 39 and 40 that -r- 

at 

changes sign on reversal of the velocities. For in those articles the pairs of 
molecules having velocities x', |', <fec. are separating after an actual collision, 
before which they had velocities x^ |, <&c. Corresponding to the classes x, ^, 
their number is Ff, In the synthetic system no collision has taken place, and 
the nmnber is f V'* 

B. 3 
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If the system were to move through that state, -i— would be 

discontinuous. 

It m ay indeed^ if the distribution of velo rri^-^'*^" KftppAna f.n te 
the normal one, I^a^^ MavwAll'a rl^pt.y^l^yjjj^Ti It cannot be one 
of a series of states through which the system passes in its 
approach to Maxwell's distribution, as the H theorem requires. 
The system then which I have supposed to be constructed is 
not, as appears to me, in Boltzmann's sense molecular-unge- 
:^|^rdnet. It should seem that the molecular ungeordnet state 
must cease to be such on reversal of the velocities. 

What the ^HLjbheorem proves t hen is this, that the distri- 
bution of velocities expressed by the equation Fy = Ff is the 
only distribution which can be permanent consistently with the 
existence, and the continued existence, of condition A or its 
equivalent. It is assumed without proof, and as I hope to 
show in the succeeding chapters is not generally true, that 
condition A can and does continue to exist. 

The rate at which H diminishes has been calculated by 
Tait, Watson, and others for certain very interesting cases. See 
Watson's work above refen*ed to. Art. 15. 

42. Boltzmann (see Berlin Sitzungsherichte, February 
1897) considers that we have in th is behaviour of the system a 
trn ^ irreversible proce ss. That term is generally applied only 
to processes involving dissipation of energy, that is conversion 
of some other form of energy into heat. If however heat, or 
temperature is to be explained as the kinetic energy of mole- 
cules or atoms moving according to usual dynamical laws, 
every process, even if it does involve conversion of other forms 
of energy into heat, is theoretically reversible ; would in fact be 
reversed if at any instant all the velocities were reversed, as we 
assumed them to be in discussing the H theorem. If that be 
so, by asserting that any process is irreversible we assert only 
that by no means within our power can we reverse it, because, 
that is, we cannot control individual molecules. It is in this 
sense that, as I understand, the diminution of H in accordance 
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with the theorem would be called by Boltzmann an irreversible 
process. 

Bo ltzmann also shows (Vorlesunpen, pp. 58 — 60) that — H 
di ffers by a constant only from the logarithm of the chance of 
the system being in its actual state, and therefore as H ^^ 
diminishes in absolute value, the system passes from a less to 
a more probable state. See note at end of chapter. 

43. A solution of the equation 
for the two sets of elastic spheres is 



/-m' 



g~Aaf(«»+j/«+2a)y 



/= (—\ €-*"*(f«+V+^) J 



(12). 



That is found by treating the kinetic energy as the only 
thing which remains unchanged by the collision. It follows 

from this solution that for the same sphere a?* = y^ = z'^, and 

that for different spheres Mai^ = m^ &c., or the mean kinetic 
energy is the same for each sphere whatever its mass. 

44. The complete solution of the equation Ff^F*f\ as 
given by Boltzmann, y(yrUsmigeay p. 131, is 

where il is a function of the coordinates, and 

Q = m{(a-i^)2-|-(/9-t;)« + (7-w)«J (13), 

Uy Vy w being constant velocities ; F has the corresponding value. 

In stationary motion f = 0, that is 






Suppose now the system to be in a field of external force 
whose potential is j^ , e.g. the column of Art. 29. 

3—2 
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If no other forces act on the system than those derived from 
X, we have 

/da du\ dy /, .x 

and therefore the equation becomes 

(a-«)(^ + 2A4g)e-»Q + &c. = (15), 

the solution of which is A = e"2*x (16). 

t 

For instance, if the molecules be elastic spheres of infinitely 
small diameter, or material points, between which no collisions 
occur, there are no forces except those derived from x> ^.nd we 
deduce by the above method, A = e"^^. 

The expression now asserts that the chance of a molecule 

Jiving velocities a ... a + da &c., and being in a position where 
he potential of the external force is X"*X + ^X* ^ 

g-2ftx e-^{(a-«)«+&c.} dadfidydx. 

It follows that the distribution of the velocities is the same for 
all values of %, that is in all parts of the system, and that the 
number of molecules per unit of volume is proportional to 
g-2fcx. This agrees with the result of Art. 29, because in the 
case now supposed /c = 0. 

46. If the spheres have finite diameter c, so that collisions 
do occur, there are other forces acting besides those derived 
from x> namely, the infinite forces assumed to act for infinitely 
short times during collisions. And it is no longer true that 

d y V dv 

and we cannot obtain the result 

except by ignoring the infinite forces altogether. It is true 
indeed that they act in equal and opposite pairs, but no 
solution can be satisfactory which takes no account of them. 
See further as to this problem post, Art 103. 
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46/ If the intermolecular forces are finite, and have a 
potential -^j we have 

dt dx dx' 

And our formula leads to the conclusion that 

\ A^€'^^^^' I (17), 

expressing now the fact that as before the distribution of 
velocities is on an average the same for any molecule, in 
whatever position it happens to be, not only with respect to 
the external field, but also with respect to other molecules — but 
that the number of molecules per unit of volume in a position 
in which the whole potential is ^fH--^ is on average e'^^+^^^jOA 
I Rf^ft no reason why W f^ should hesitate to accept this result, or 
ho w we can avoid accepti ng jt^ ftYOftpt that Dr Watson (in the 
work above referred to, pp. 70, 71) hesitates to accept it, 
though I do not understand him to deny its validity. He 
puts the case that, with finite intermolecular forces, a molecule 
P may have a sphere of action, within which may be many 
other molecules. Any individual molecule outside of the 
sphere of action exerts on P a force which is negligible. But 
the infinite number of molecules outside of the sphere of action 
may exert on P a sensible force. And he shows that this may 
be included in the fixed centre forces. But he refrains from 
drawing any conclusion for the molecules within the sphere of 
action. 

« 

Boltzmann a nd WjaisoiCs generalisation. 

47. Our molecules regarded as elastic spheres have each 
only three degrees of freedom, that namely of motion of trans- 
lation. 

Let us_now _ attribute to our molecules the most general 
fo rm possible. 

Let a molecule m be a material system, whose position and 
.^ state at any instant are defined by n generalised coordinates 
jfi ... jn, and the corresponding momenta |)i ...|)n« 
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Let pm be the number of molecules of the kind m per unit 
of volume. 

Let pfn/(qi'"qnPi'"Pn)dqi...dpn. 

or shortly pmfdqi . . . dpn , 

denote the number per unit of volume of molecules m for which 
the variables lie between the limits 

?i ... qi + dqi 



?n ... qn+dqn 

Pi ... Pi + dpi 



/ 



(c), 



Pn -*-Pn+dpn 

or, as we shall express it, are in the state c. 

We may express this in a different way that i^ chance of 
an m molecule being in the state c isfdqi ... dpn» 

It is assumed that the molecule is, so far as its own 
internal forces are concerned, stable. That is that if it be now 
in the state c, and if no forces except its own internal forces 
act upon it, it will after the lapse of some time, greater or less, 
find itself in the state (c) again, except so far as the coordinates 
defining its position in space, which are included in qi...qn, 
are changed by its motion of translation. It is also assumed 
that/ will not be altered by reversing the sign ot pi..,pn 
simultaneously. And therefore /cannot contain odd powers or 
products of the p*B, 

In the same region let there be another set of mole- 
cules M, e ach having r degrees of freedom, whose coordina tes 
and momentaare denoted respectively by Qi... Qr and Pi ...P^. 
And about these make the same assumption as above stated. 

Let pjg be the number of M molecules per unit of volume. 

Let PmF{Qi... QrPi ... Pr)dQi ... dPr, 

or shortly pj^ FdQi . . . dPr , 
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denote the number per unit of volume of molecules M for 
which the same variables lie between the limits 



Q,...Qi + dQi \ 






> 



(0), 



Pr...Pr+dPr) 

or, as we may express it, are in the state C. 

48. Let ^ (^1 . . . 9n Qi • • • Qr) be a function of the coordinates y^ 
of two molecules m and M, such that when <f> is positive there 
is no mutual action between the two molecules. 

Let us now make an assumption equivalent, with respect to 
the molecules now under consideration, to j;ondition A . Let 

us assume namely that so long as if> is positive and ^ negative, ^ 

the two chances F and /are in dependent. And therefore the 
number per unit of volume of pairs of molecules, one belonging 
to the rns and the other to the if 's, for which at any instant 
the variables lie, for m within the limits c, and for M within 
the limits C, is 

pmpM^f^i ••• dpn^Qi ... dPf 

When ^ passes through zero the two molecules pass into 
their sphere of mutual action, and an encounter commences 
between them. The chances F and / generally cease to be 
independent. 

It may be the case that ^, having become negative, would 
if the two molecules are not influenced from without, remain 
negative for all time. The two molecules, that is, may, each 
remaining stable, form a permanent union revolving round 

/6ne another in a closed orbit. I shall not consider that case, 
but assume that after a certain time t, <f> will again pass 
through zero and become positive. The encounter ceases when 
^ again becomes zero. All those pairs of molecules for which 

at a given instant <f> lies between zero and ^dt, -^ being 
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positive, will within the time dt after that instant commence, 
and as I am now assuming, will in due time afterwards com- 
plete, an encounter. The number per unit of volume and time 



of such encounters must therefore contain the factor 



d4 
dt' 



/ 49. But the two molecules have between them only n -|- r 
degrees of freedom. We must therefore use ^ for one of the 
generalised coordinates defining the position of the double 
system. 

Let us substitute ^ for g^. Then all those pairs of mole- 
cules, m and M, for which at any instant the variables other 
than qn lie within the limits c and C respectively, and for 

which at the same instant 6 lies between zero and -^ di. will 

^ dt 

in the interval of time di after that instant encounter one 

another. The number of such encounters which commence in 

time dt is therefore 

PmPM^f^i • • • d^n-idpi . . . dpndQi . . . dPr -^ dt 

As a consequence of the encounter, after the time t the 
variables will be found between the limits 

?/ qi-^dqi \ 



for m, and 



q n—i • • • y n— 1 "1" dq jt_i 
Pi' Pi'+dpi' 

Pn Pn' + dpn ^ 



Q^.^.Q^ + dQ^ \ 



(c') 



J^f * . • JTif "t" CvJr^ 



.(C) 



for M, and <f> between zero and -^ dt, and the encounter 

ceases as <f> again passes through zero. The pair of molecules 
pass out of their sphere of mutual action. 
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According to our notation the number per unit of volume 
of pairs of molecules, m and My for which the variables lie 

between these limits c', C, with -^ positive, is 

pmpM^'f'dqx ... dq'n-~idpi ... dpndQi ... dP/ -^ dt 

Conversely if at the beginning of the time t the variables 

were in the limits c', C, ^\ -?- dt with reversed velocities, 

they would at the end of the same time t be found within the 
limits c, C with reversed velocities. We may call these reversed 
encounters. The number of reversed encounters less the number 
of direct encounters per unit of volume and time is 

pmpM^fdqi . . . dgfn^idpi . . . dpndQi . . . dPr -^ dt 

" pmpM^fdqi . . . dqn^idpi . . . dpndQi . . . dPr -j: dt 



60. I t can no w be pro ved (tl^e sy st em being conservativ e) 
that 

dqi ... dq'n-idpi ... dpndQi ... dP/ -^ 

= dqi . . . dqn-idpi . . . dpndQi . . . dPr -^ . 

For the accented variables qi, &c. are all functions of the un- 
accented variables and of the time t. And as r is the same for 
every system within the initial limits, we have 

dS , dS , , 
d^r^'' d^r^''^'" 

where / SLis Hamilton's Principal function for the moti on in 
question. Similarly 



ds ds - /#^L^S^ 

dqi ^ , dq^ ^'' (university 



£Lc>^liforh^ 
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Now writing d<f> for -^ dt and dff/ for -^ dt, 

dqi ... dq'n-\d^'dpi ... dpndQi ... dP/ 

= (ig/ ... d^fi^id^'dQi ... dQ; x %' ... dpndP^... dP/ 
= dqi ,.. d^n-\d<f>dQ^ ... dQr'dqi ... dqn~id<f>dQi ... dQr 

X S + ^' -^ ^-P/ 
dji rfg'j '" dQr * 

where the last factor denotes the functional determinant 

dpi dpi dpi 

dqi dq^ '" d<f> 

dpj 

dqi 



and since 



, dS dpi d^S „ 
^''dqr dq, dq^dq^ ' 



the Ainctional determinant may be written 

d«S d^S d^S 



dqidqi dqidq2 *" d^id<f>* 

d^S 
dqidqi 

In the same way 
dqi ... dqn-id<f)dpi ... dpndQi ... dPr 

= dqi . . . dqn^id<f>dQi . . . dQr x dqi . . . dq'n-id^*dQi . . . dQr 

and the functional determinant has, but for its sign, the same 
value as before. Hence the continued products of the differ- 
entials are numerically equal, that is 

dqi . . • dq[n-i d<f>dpi . . . dpn dQi . . . dPr 

= dqi . . . dqn-.id<f>dpi . . . dpndQi . . • dPr (1S)> 
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and therefore the excess of the reverse over the direct encounters 
per unit of volume and time is 

pMpmiFT"^/) ^1 ••• dqn-id<l>dpi ... dpn X dQi ... dPr, 
or 

JTjf 

:jt = pMpm dQi... dPr jJSiF'f - Ff) dqi... dqn^id<l>dpi ...dpn ; 




and is zero if in all cases Ff^ = Ff. 



A corresponding equation holds for -^. See further for this x. 

proposition Watson's Kinetic Theory of Oases, second edition, 
Arts. 8, 10, 12, 13. 

61. As this property is true for all values of F and f with 
the corresponding F' and /', we see tha t if Ff^ ==Ff in all cases 
wh ere a pair of molecules can by enco unter pass out of the 
classes F, f into the classes F\ f respectively, or vice versa, 
we have a s ufficient condition for statio nary moti on. 

The H theorem of Art. 38 can easily be adapted to prove 
that it is also a necessary condition. Further we may take 
for solution of the equation F'f^ = Ff, 

where T is the kinetic energy of the molecule M, and T that of 
m ; and if 7^ T' os\r\ Kp rpHnppH to thp '^"^ /^^ g^Tl— ^"-JlfJ^^^ 
v elocities, we get a result analogous to that obtained for elastic 
spheres in Art. 43, namely that, as a consequence of conffilidn -^^ 
A, the mean kinetic energy is the same for each degree of 
freedom. 

62. It is necessary here to point out what further 
assumptions, be sides condition A, we have made concerning 
our molecules in the above demonstration. We expressly 
assumed in Art. 47, that each molecule l eft to it^fi l f is sdnhU ,^ 
that is, that so far as the mutual forces between its parts are 
concerned, the coordinates and momenta, whatever values they 
have at this instant, will or may, at some future time pass 
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simultaneously through the same values or infinitely near to 
them, the coordinates x, y, z defining" the position in space of 
the centre of inertia being exceptedv/We also expressly assumed 
that if two molecules pass into an encounter, they pass out of 
it again, that is do not form a permanent union. We have 
further assumed tacitly that each molecule is stable, not only 
under the influence of its own internal forces, but also under 
the influence of the new forces which may be developed during 
its encounter with another molecule. Suppose for instance for 
an m molecule /= e"**, where « is a quadratic function of the 
momenta ^1 ... ^n with coefficient functions of gi...gn. A 
necessary condition of stability is then that 8 shall be positive 
for all possible values of the ^'s. And that requires that the 
determinant of the coefficients in s, which we will call d, and all 
its coaxial minors, shall be positive. See Appendix {q). 



If 



s = a^^ + K2P1P2 + &c., 



then 



d = 



Oi 61a 613 

V12 flf'2 ^^28 



Similarly in order that an M molecule shall be stable if left 
to itself, we have if 

ir'=e-^«, and S=A^?^'VB,JP^P^-\-&Q,, 



D = 



Ai B12 -013 

i)l3 A^ X>2S 



and D must be positive. And so long as the two molecules 
do not influence each other, the condition is that Dd shall be 
positive, which is necessarily true if it is true for d and D 
separately. But when the two molecules are within the sphere 
of their mutual action, the corresponding criterion of stability 
is that s + S shall be positive whatever values the velocities may 
have. If A be the determinant of the coefficients in 8 + 8, 
A, and all its coaxial minors must be positive. 
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But if we fofm the determinant A for the coefficients in 
8-\-8, it will generally during an encounter or when ^<0 
contain constituents of the form pP. The condition may fail 
for A, although it held when ^ > for d and D separately. 
So that the two molecules, though stable separately, may be 
unstable together. The encounter may be the death of the 
molecule. 

Note on the ff Theorem, Art. 42. 

/ The R Theorem proves strictly that, assuming the independence of 

the chances (Art. 17), -j- is generally negative. It is claimed for it, as a 

physical theorem, that the diminution of jET is analogous to the increase of 
Entropy in irreversible processes. To make that analogy complete, the 

dTT 
chance of -j- being positive ought to diminish indefinitely as time in- 
creases. I will state the difficulties that this view presents to me, without "A 
saying that they are insuperable. 

A finite number of molecules in a finite space are set in motion at an 

initial epoch at random, in the manner described in Art. 41, and then left 

for infinite time, undisturbed from without, under conservative forces. 

dlT 
Initially --j- is almost certainly negative. After time t the coordinates 

and momemta of every molecule are determinate functions of t and of all 
the initial coordinates and momenta. If (7, C be two successive states of 
the system, each in its turn a consequence of the same initial state, in 
what sense is C less probable than C ? It may be that the initial state, 
formed as we have formed it, is less likely to be C than C But the 
chance of state C at time t is not proportional to the chance of the initial 
state being C, It is then improbable that the initial state should be one 

in which -^ is positive. It is not proved to be improbable that it should 

dir 

be one, a necessary consequence of which is that -17 will be positive at 

time t. 

Again, the continued independence of the chances, if conceded, 
attributes to every subsequent state of the system the same properties 
as the initial state has, which would be right if the system were con- 
tinually receiving external disturbances. Can it be conceded for the 
isolated system ? Let aS' be the Principal function, a;, of type coordinates 

of two molecules at time t. It x=af nearly, are the momenta -r-, -rn^ 

correlated or independent? I think it can be proved that they are 

dF 
correlated. If so, the formula 8 of Art. 37, expressing -^- , is inexact. 

dS 
And the theorem fails to prove -^ negative, except when the factors 

F'f -Ff&re great^ that is when If deviates widely from its normal value. 



ass 



CHAPTER IV. 

CORRELATION OF VELOCITIES. FINITE FORCES. 

63. I PROPOSE to prove in this and the next chapter that 
in a system consisting of molecules of finite dimensions in 
stationary motion, it is not true for molecules very near to one 
other, that the chances of their having velocities between 
assigned limits are independent, ss condition A assumes : but 
on the contrary if the forces be repulsive, they tend to 
move on average in the same direction ; and therefore that in 
dealing with their encounters, or (in the case of rigid bodies) 
their collisions, we cannot legitimately base our conclusions on 
condition A, except in the limiting case when the density — 
i,e, the ratio which the aggregate volume of all the molecules 
in unit of volume bears to unit of volume, — is infinitely small, 
though conclusions based on condition A are approximately 
accurate for gases under ordinary conditions. 

In the present chapter I shall deal only with the case of 
finite intermolecular forces. 

Let m, &c. be the masses of a number of molecules in 
stationary motion. Let os, y, z denote the coordinates of m, 
oi^ y'y z* those of any other molecule as m'. Let r be the 
distance from m to m'. Let / be a function of r which is 
ever3rwhere positive, finite and continuous, which is of negative 
degree except for very small values of r, and such that 

^m^fdr is finite. 



Jo 



Then for any m, Xf is the sum of the values of/, referred 
to the position of m as centre, for all the other molecules; 
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and if p be the number of molecules per unit of volume, 
^TTT^pfdr, Also 2r -4- has corresponding meaning. 

Eoetendon of Clausms* Theorem. 

64. We can now put the Clausian equation in a modified 
form as follows, 

for the function within brackets has a determinate value at 
every instant, and cannot on average increase or diminish with 
the time, the motion being stationary. 

Again for the same reason, 

i {2m| 2/.'+ S«.| W + 2m|2A} = 0...(II). 

Here the first S denotes that every molecule successively is 
to play the part of m, and the second S includes all the 
molecules except m each with its own x\ y\ /, and each with 
its own /referred to the position of m for the time being. 

• Subtracting I. from IL we obtain 



-H^<ti*[%)Vf 



+ 2m|2(.'-.)f=0 (III). 



Again, At-T^U 



dt drdt' 
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And if X, /i, 1/ be the direction cosines of miniy or r, drawn from 
m to m', a?' — a? = Xr &c., and 

dr _ daf A\f dsf 

di'^Tt^^'di'^'''di 

(dx dv dz\ 
and the last line of III. becomes 

That is, since X/x = Xi/ = /ai/ = on average, and 
on average, the last line of III. becomes, 

= 2/ by the conditions for /, if -^ = 0.* 



dr 



Again, m -^ = X, &c. 



if X, F, ^ be the components of all the forces, whether external 
or intermolecular, acting on m. 

* If the molecules be scattered uniformly ^=0 on average. If not 

uniformly, it will be generally negative, because we take a molecule for centre, 
and there are more molecules where p is greater. But (see Art. 106) it will be 
small. If negative it helps my present argument, but I shall treat it as zero. 
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Therefore m -t— 2/a?' — a? = X^fiso - x), 



Making these substitutions, and arranging the terms, we find 

= -X/..{(|)V(|)V(|)] 

Vac ar ar ar at at/ 
+ 2Z2/(a/ - «) + 2F2/(y' - y) + 2Z2/(y - ^) 

mi ^ • ^ v^ /. /d^ dx' dy dy' dz ds!\ 

= (IV). 

66. Now if there are no intermolecular forces, 
2Z2/(a/ - a?) + S F2/(y' - y) + 2^/(/ - ^) 

is for an infinite system negligible, and equation IV. is satisfied, 
whatever/ may be, by making 

dx daf dy dy' dz dsf _^ 

dt dt dt dt dt dt 

« 

J L dx ^ dy , dz\ A dx' . dy" dsf\ ^ 

on average for each value of r, and therefore for each value 
of /. This expresses condition A, which thus appears as a 
consequence of there being no intermolecular forces. 

B. 4 
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66. Bu t if the r e are intermolecular forces, we shall find 
that condition A cannot ex ist. 

In the limiting case when the encounters are binary, or no 
molecule is within the sphere of action of more than one other 
molecule at the same time, we have 

Xtf{x' -x)+ Tlfiy -y) + Zlf{^ -z)^- Rrf. 

In the general case consider a spherical shell of radii 
r ...r-V dr described about m. Let R be the intermolecular 
force between m and a molecule m' within the shell, the angle 
made with x by the line mm' drawn from m to m'. Let fr be 
the value of / for the radius r, and let 2^ denote summation 
for all molecules m within the shell. Then if X, &c. relate 
only to the forces exerted on m by molecules within the shell, 

Xy{x' - a?) = "frRr (2r cos Of, 
the mean value of which, if the molecules be distributed through 
the shell at haphazard, is —^Sr^^. 

Again X includes, not only the forces acting on m due to 
the molecules in the shell r ... r 4- dr, but also all other inter- 
molecular forces acting on m, e,g. forces due to molecules 
within the shell r^ ...r +dr\ which forces we will denote by 
R. The part of Xlf(x' — x) due to these forces is 

—frRr%r COS ^Sr' COS ffy 

where ff corresponds to above. If we assume 

2r COS ^Sr' cos ^ = (1), 

on average, we shall have 

But whether X%f{x' — x) be equal to - jSriZ^or not, it is 
nevertheless negative, and not zero. What has been proved 
for X is equally true for T and for Z, So that 

nf{x' -x)-{- TXfiy -y) + Ztf{z! - z) 

is in any case negative, and in the limiting case of binary 
encounters, or wherever we can assume (1), is equal to — ^Rrf, 

Since, as we have seen, in equation IV, the term 

SX%f{x - a?) -f &c. 
is negative and not zero, it follows that the sum of the remaining 
two terms, that is, 
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^"^^Kdt dt ^ dt dt ^dt dt) 

. ^ x* ^f 1^ ^ . dy , dz\ A dx' . dy' . dz\ 

is positive and not zero. But this is inconsistent with condition 
A, It follows that we must discard condition A when inter- 
molecular forces exist. 

67. First solution of equation IV. Let for any value of r 

/dx ^ da! dy ^ di/ dz^ dz'\ 
"^U^'^+d«^- rf« + d«^'- d-J = *'^^''* on average, 

where ^ is a function of r such that 47rr'^, and therefore also 

46in*^fy vanishes at either limit, when r is zero or infinite. Also 

djSc dvf 
assume all terms of the form -77 -^ , &c. to be zero on average. 

Then as in Art. 54, 

K> df l^dx ^ dy dz\ (^ dx\ dy^ ^ ds!\ 

= J2<^ -i- , because V = /t^^ _. ^ _ ^ 

= -/ 47rr></i/ar - ^ j 4i7n^p -^fdr, 

and using this in IV, we find 

2ZS/ (x -x) + &c. - JS r47rrV ^fdr = 0, 

which is satisfied by assuming for every r 

Ir {X (x' -x)+ 7(y' -y) + Z(z -z)}=i*m*p ^ ; 



dr 



also I 4nn^p^dr = — i I 47rr*p -^ dVy 

/•OO 

whence I inrr^pdydr = ^Rr, 

in all cases in which we can assume (1). 

Second solution of equation IV. Since / is arbitrary we 

may equate separately to zero the term in IV. involving r -^ , 

4—2 
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and the sum of the two terms involving / That is, we make 
for each value of r, 

^ f^ dx , dy ^ dz\/^daf, dy' d/\ ^ ,^^ . 
J [dx^ dd dv ^ dvf dz ^ dz'\ 

^^ "^Vt^^-d^^dt^^-dt^lJt^W^ 

+ Z2,(a?'-a?)+F2.(y'-y) + ^2,(/-^)=0...(IVb), 

in which X includes the intermolecular forces acting on m due 
to molecules at whatever distance, but a! — on, &c. relate to a 
particular distance r. 

The result (IV b) and the first solution show that the 
velocities of m' molecules distant r from m are correlated 
with that of m, or 

fdx ^ dx dy^ dy' . dz ^ dz'\ ^ 

And the corrdaticm is a function of r. 

In this solution 4t7rr^p<l>, and in the former one — J47rr*p ^ , 

is a determinate part of the Virial of the intermolecular forces 
acting on m, though not necessarily equal to Br, except in 
cases where (1) can be assumed. 

When r exceeds the " radius of action " of the molecule m, 
and 80 R = Oy X nevertheless includes the x force due to 
molecules distant r' from m, for which R is not zero. But 
ultimately as r increases, S^ cos 6 Sr* cos ^ = for all values of 
r' for which R is sensible. So that as r increases, 

dtX ujX 
m -IT 2r -tt + &c. becomes in any case zero. . 

But the correlation may possibly retain a sensible value at 
distances from m greater than its radius of action. 

Equation (IV a) expresses that although the velocity of m 
is correlated with that of m\ yet it is not on average correlated 
with the velocity of m' resolved in direction mm\ that is not 

correlated ^^^ X ^^' + /li ^ + 1; J^ . For (IV a) with (IV b) the 

simplest, though not perhaps the best, solution is the following. 
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If mP represent in magnitude and direction the velocity of 
m, and if it be taken for polar axis of a sphere of radius r 
described about m as centre, then the average motion of 
molecules on the surface of that sphere is on the meridian 
towards P, and proportional to the sine of the colatitude. 

68. If the force R be attractive instead of being repulsive, 

dx due 

-77 -^ + &c. will be negative instead of positiva The velocities 

of m and m' may be said in this case to be contrarelated. If 
there be two sets of molecules, A and B, and every A repels 
every other A , every B repels every other 5, but A and B 
attract each other, the system will tend to form at every point 
two streams, the A*8 moving in one direction, and the JBs in the 
opposite, which we may compare, 

* si parva licet componere magnis,' 

with the mutually attracting and repelling electricities supposed 
to constitute an electric current. 

It should be noted that for any molecule m. 



^"'^fed«+**^)' 



the summation including all values of r, is if (1) be assumed 
equal to the Virial, taken as positive when the forces are re- 
pulsive, of the intermolecular forces acting on m. (See Art. 85.) 



doD djSC 
I have treated -^^ -^ as zero in the absence of intermolecular 

dt dt 

forces. Strictly, n the number of molecules in the system being 

finite, and the centre of inertia at rest, it must be negative, but 

it may be neglected when n is great. 

69. I have thus proved that if there be finite intermolecular 
forces, condition A cannot prevail, and 

dxdx' dy d^ dz d/ 
dt dt dt dt dt dt 



54 KINETIC THEORY OF OASES. [CHAP. IV. 

has on average a finite value, a function of r, and is positive if 
the forces are repulsive. The investigation does not apply to 
elastic spheres, in which the intermolecular forces are infinite 
during collision. I shall therefore devote the two following 
chapters to proving the same thing for equal elastic spheres by 
an independent method. It will be, I regret to say, a much 
more laborious method. If however the reader is now satisfied 
that I have established my case for finite intermolecular forces, 
he will have little diflSculty in provisionally taking for granted 
that the same thing can by some means or other be established 
for the limiting case of elastic spheres, in which the forces 
become infinite for very small values of r. And omitting 
Chapter v, proceed at once to Chapter vi. 



CHAPTER V. 

C ORRELATION OF VELOCITIES. ELASTIC SPHERES . 

60. I PROPOSE in this chapter to prove for elastic spheres a 
proposition corresponding to that proved in the last chapter for 
molecules between which finite forces act. The elastic spheres 
shall be each of unit mass and diameter c. The same pro- 
position for spheres of unequal masses or diameters is much 
more laborious, and when proved does not bring us much 
nearer to any probable form of molecule. 

Equal Elastic Spheres. 

Assume the chance that a group of n spheres shall simul- 
taneously have component velocities 

OLi ... flti -+• doii 

A ... A + dA 1- *^r the first 
7i ... 7i +^71 



On ... Oin -^ ddn 

0n'"0n + dfin\ for the nth 

Jn ... Jn + dyn 

to be Cc^^dai . . . dyn* 

If we assume condition A, we virtually assume that 
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Now it will be found that when the diameters are finite 
the motion is not stationary when 

but becomes stationary when 

Here a, a', &c., are component velocities at the instant con- 
sidered of any pair of spheres, and 6 is a function of the distance 
between them at that instant which, except for very small 
distances, is evanescent. I admit of course that given con- 
dition A always existing, the distribution of velocities among 

^ the spheres denoted by €~*^ (*'+^"+>*> is not disturbed by collisions. 

A^ That has been proved many times. But the point is that 
condition A itself cannot continue to exist when collisions 
take place. In order that it may continue to exist, it is 
necessary that the system be continually receiving disturb- 
ances, or else that the diameters of the spheres shall be 
infinitely small, so that no collisions occur. 

Definitions. 

61. Let P be any point, r the distance of any other point 
from P. Let / be a function of r satisfying the following con- 
ditions, namely, 

(1) /is equal to unity at all points within a small sphere 
of radius a described about P as centre. 

(2) /is finite, positive and continuous everywhere. 

(3) / is of negative degree at all points outside of the 'a' 
sphere and such that the series denoted by 2/*, or I 4nn^fdr, 

is ultimately convergent, but contains no high powers of - . 

We may take for 'a' the radius of a sphere which on 
average contains one molecule. 

Let now f = v^ » ^ ^ ^v » ^ ~ "W ' ^^^ summation in- 
cluding all the molecules of the system, and / having for each 
molecule the value proper to its distance from P. Then f , ty, ^ 
are continuous functions of x, y, z, the coordinates of P. 
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LetM„=ffjdadydzay[f^ + '^£j . in which {§ + ^is the 

value of that function at a given instant within the element of 
volume dxdydz, a'y is the value of the product 07 for the 
molecule (if any) which is within that element, or the sum of 
those products for all molecules within the element if there be 
more than one, and the integration includes all space. Then, 
for any element of volume, 07 changes from instant to instant 
by the passage of molecules into or out of the element. 

There are two analogous functions, namely. 



,Jf„-///«j,J,A,(g + |), 



and also three others of the form 



* Maa}=ljldxdydza^-^&c. 



It will be sufficient in this chapter to deal with M^zy and we 
may for the present drop the suffix, and write M for M^z. 



Swnmary of the Argument 

62. It is necessary first to treat of the case in which the 
molecules, though possessing finite mass, are of infinitely small 
dimensions, so -th at collisions between them will not occu r. 
That being the case, the system will be in stationary motion if 
the velocities are , distributed according to Maxwell's law, Q 
being a sum of squares only. This we may admit whatever 
views we may hold as to the possibility of any other law, I 
shall assume Maxwell's distribution to exist. It follows from 
Q being a sum of squares only, that the feet of any molecule, or 
any group of molecules having any given velocities, affords, if 
the whole number of molecules be infinite, no presumption 
whatever with regard to the velocities of any other molecules. 



^ 



\ 



58 KINETIC THEORY OF OASES. [CHAP. 

We then find for this system the mean values of f", ( -i^ J , &c. 

JTLf 

In this system of material points —jt = 0, which gives rise to 

equations A and B below. It is then proved that if, with the 
same masses, the spheres have finite diameter c, so that collisions 

occur, -7- , which was zero, becomes, if f", &c. retain the same 

mean values as before, positive, and the motion not stationary. 
But the motion must become stationary. Therefore in the 
system of finite spheres p, &c. cannot retain the same mean 
values as before. It follows that in this system Q in €"*^ cannot 
be the sum of squares of the velocities only. 



A System of Material Points. 

63. The chance that the molecules shall at any instant 
have for their component velocities in x, ai ...ai+dai, &c., is 
A€~*^dai ••• da», where -4 is a constant and in the system now 
considered 

Q = ai« + Oa' + &c. 

But from Art. 61, 

^^n = f 2/"— Uifi — ... — 0«_i/n_i. 

It conduces to simplicity to assume /n=l> but, as will be 
seen, it in no way affects our conclusions. That being assumed, 

«» = ^y- ai/i - 02/2 - &c. 

Substituting this value for a^ in Q, we obtain 

Q = r (2/)» + (1 +/1O a,' + (1 +/,») a,' + &c. 

+ 2/i/aaiaa + 2/i/,aia8 + &c. 

If we now integrate €~*^ for ai ... Or^i successively with limits 
± 00 , we get the result next stated. 

Let D be the determinant of the coefficients in Q, each of 
the coefficients of products being divided by 2, that is 
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D = 






Let Dpq be the minor formed by omitting the pth row and 
gth column. Similarly let Dp^n be the minor formed by 
omitting rows p and r, and columns q and s. Then by inte- 
grating for tti ... ««-!> Q is reduced to ^y" f*> whence ?*= oryj- 

See Appendix (gr). Also by integrating for Og ... a,^_i only, Q is 
reduced to 



4#^«x'+^«.?+i 



i) 



1221 



D, 



■D« ^ 



'1321 



'1221 -^ 

Again we find (see Appendix {g)) 

Ai = 1 +/i^ 4-/2" + . . . / Vi = S^y « since /„« = 1, >^ 

D^ = 2/» -/l^ 

Ai28 ~y 1/2* 
i)„ = 2)„ = &c. = D = (2/)», 



\ (19)- 



We have then for the chance that a^ and ^ shall lie within 
assigned limits irrespective of the values of aa...a»_i the 
expression Ae^^^daid^ in which 



and 



tp /iS/ , (2/y 

2/» -/i» "" '' 2/' -/i» ■" -^^ 2/» - A' "■ ''• 



We see now that the assumption ^n^ 1 has not affected the 
result, because writing 2^ for/* does not alter 2>, 5, or r. 

64. Let now f be the value of that function at P the 
origin, f + 9f its value at a neighbouring point P'. Then at 

P', 6"*^ becomes 6""*^ ^^^Jh ^~^^* -^^ ^^ f ^® ^®^y small, this 



is equal to 



€-^<i (1 + 2hpa,d^), 
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But if at P, f is not negligible, we must use instead of 
a, fi, 7 the following functions, namely, 

u = a — ^, v^fi-rj, w = 7--f, whence a = t* + f &c., 



and ^ = ?(ti + f)»-2pt4 + f f + rp, 

or Q' = qu^ + 2q-p^ + (? - 2p + r) f*, 

and let c"*^ = <^ (u f ). 

Then if ^ (t^ f ) be the value of this function at P, ^ (u f ) its 
value at P', 

Again, let if>(u)du be the chance that at P w shall lie 
between u and u + dw, whatever f may be, that is 

J —00 

and the same chance at P' is 

Also ^^(„f) = (^)^(«f)_|-^(„|), 

if f -^1 denotes differentiation with u constant. Therefore 
^, (u) -*(«) = d^fj^ (^) <^ („ f ) _ af I'^df 1^ ^ (u f ). 

because ^ (u^) vanishes at either limit. Therefore 
^ (tt) - <^ («) = - 3? |"^df J^ <^ (« ?) 

J -00 

And this is proportional to 
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if h'^h 



wr 



whence <^ (w) = <^ (u) (1 + ^h'ud^). 

65. Now with the above definition of M we may write 

-37- = 0, that IS 
ait 

the summation being for all the elements of volume, each with 
its own a7 and its own ( 77 + 7p ) • 

Now we shall find that on average -r- arf is of the opposite 

sign ^o ( 7p + j^ ) » ^^^ 77/ ( 77 "*" 7F") ^® ^^ *^® same sign with 

urf. So that -37 consists of two series neither of which is zero, 

dt 

but which are equal and opposite. 

66. Let us consider the plane of xz, and ^ as varying 
in direction z, f as varying in direction x. Let us first deal 

dp 

with -^ only. Consider the two infinite planes z==0 and 

di 
z = dzy and suppose for a moment -^ constant between them. 

The number of molecules at or near the plane z = dz for which 
tt, or a — f, lies between u and u + du exceeds the correspond- 
ing number at or near the plane -sr = by the quantity 

Ae-^'^duih'u^dz. 

az 

That follows from Art. 64. To fix the ideas let u be positive. 
Let us separate from this class of molecules those whose z 
velocity is 7, and suppose 7 negative. Let us then call the 
separated class the class (v/y\ Then the number of molecules 
of the class (^7) which pass through the plane z=^dz per unit 
of area and time exceeds the number of the same class which 
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pass through the plane z = per unit of area and time by the 
quantity -4 €"*'*** e--*^' du2h'wy-^dz. And the number of mole- 
cules of the class (t^) which enter 8 the space between the 
planes per unit of time exceeds the number of the same class 
which pass out of S per unit of time by 

which is positive because v/y is negative. 

If SNdvdy be the number of the class (tuy) within S, we have 

^ dudy = - Ae'^'"^ e-'^y* dudyih'v^ ^ (20). 

But if ^7 be the mean value of vr/ for all the molecules within 
S, u and 7 now taking all possible values, 

uy = Ij Nv/ydvdy / 1 1 Ndudy, 
and therefore 

because -7- / \Ndudy = on average, 

= - ^ jJ€-^''''€-^y'2h'uYdudy ^ by (20) 

2hd'z' 
I have assumed in this proposition that the space 8 within 
which -^ is constant, is the continuous space between two 

infinite parallel planes. That will not occur in fact. But it is 

not essential to the argument that S should have that shape, 

or that it should consist of one closed surface. The proposition 

holds equally if S consists of all those elements of volume for 

dp 
which at a given instant -^ has the same value. And it is true 

on average for each element of volume. 
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d— d 



67. Again, ^^^^(a-f)^, 

d — d *r- 

Now a molecule moving within 8 with positive 7 moves 

from less to greater f, -^ being positive; if with negative 
7, it moves from greater to less f. For all molecules moving 

within /S, -^ 7^ is positive and equal to 7* -1^ . If /S be divided 

into parallelepipeds of base dxdy and altitude PP* 

///t* f dxdydz =jji' i^p - ^p) dxdy, 

^p and ^iK being the values of f at the upper and lower surfaces 
of the parallelepiped. By the motion of all these molecules 
7^ is increased. 

But at the upper end where f is greater, molecules with 
positive 7 pass out of the space, at the lower end they pass in. 
For negative 7 the converse is the case. By these entries and 

exits 7^ is therefore diminished at the rate / 17* (fp - ^jy) dxdy, 

or by the same quantity by which it is increased by the motion 
of the molecules within S, On the whole therefore if f alone 
vary 

^"'^ dt'^^dt'^-^dt'^^ 

d — 
^dt"^- 

We have then from Art. 66, 

d — __ l^df 
d< "^ ~ 2hdz' 

Similarly, if f only vary, ^^ S^ = - ir g , 

and therefore if f and f both vary, 

d ,-, 1 (d^^d^\ ,^,, 
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and (A) becomes 

///"^i(S-S)i«^-i(i*sr}-« -(B^ 

It follows that «7 ji ( 3^ + ;p] must be on average positive, 
and equal to ^(g + g)V 

68. We must however consider the value of -r: (-r- + ^^ ) 

at \dz dxj 

more carefully thus : 

Let P, P' be two points on a sphere 8 whose coordinates 
referred to the centre G are respectively x, y, z, and x, y, — z, so 
that PP' is an edge of the parallelepiped, whose base is dxdy 
at X, y. Let fp, fp' be the values of f at P, P' respectively. 

Then if S be the volume of the sphere, and -^ the mean value 

df 
of -^ within the sphere, 

the integration dxdydz being throughout the space 8, Let the 

radius of S be * a ' (Art. 61). Then we may take -^ as equal to 

^atO. 
dz 

Suppose a molecule at C whose velocities are a, fi, y. The 

d^ 
change with the time of 8-^ due to the motion of that molecule 

^ V dx"^ ^dy'^'^ di) jj ^^^ ^^^ ^ ^^^' "^^^ change due to 
the motion of all external molecules is zero on average. 

Let now / be the value of / referred to P as centre. Then 

fp = -^ . And let /' be the value of / referred to P' as 

centre, so that f p* = -^ . Then due to the motion of the 
molecule at C 

d ,^ f. X 1 f d ^ ^d ^ d\[df df\ 
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because on average -i-2/=0 &c. The minus sign is here 

required because / is referred to P, and therefore -^ has the 

opposite sign to that which it would have if referred to (7. 

Thatis, f=-?f. 
dz rdr 

Now to form equation (B) we shall have to multiply by a^. 
And that will make the terms in a, /S disappear, because 
a^ary = a)3a7 = 0, on average and therefore on average 

And therefore 

-»^si-«v|^/Mf-f)' 

the integration being over so much of the plane of anf as is 
within S, that is "• *7 ;^ ;j^ = *V s?"^ ' ^ being the value 

of^whenr=a. 
ar 

By symmetry 

'dtdx '^ 'Sfada' 

It follows from (A) and (B) that in stationary motion of 
the medium of material points on average 



, Case of Spheres with Finite Diameters. 

69. We now pass to the case in which our molecules, 
instead of being material points incapable of colliding with each 
other, become equal elastic spheres of finite diameter c. It will 
be shown that as the result of their collisions with each other, 

-J- , which in the system of material points was zero, becomes, 
long at least as f^, f ^ J , &c. retain the same values as before. 



so 
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positive, and the motion not stationary. It is worth while to 
show first in a general way by the aid of a diagram how this 
effect is produced by the collisions. 



$i 





Let us take the plane of the paper for the plane of xz and 
suppose f alone to vary in z, -r- being at a given instant 
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positive neai' a certain point. That being the case, there is at 
the instant in question near that point a mean motion, or 
stream, of the molecules in direction x increasing as z increases. 
We may suppose this instantaneous stream to be positive for 
positive values of z, negative for negative values of z. It is 
indicated in the diagram by the long arrows. 

Let us consider two kinds of collisions, (1) and (2), the 
directions of motion of either sphere before and after collision 
being shown by the small arrows. For simplicity these directions 
before collision are here drawn parallel to a?. 

If we examine collision (1), we see that the molecule which 
before collision was moving towards positive x is deflected 
upwards, i.e. towards the positive stream, and the molecule 
which before collision was moving towards negative x is deflected 
towards the negative stream. In either case the efiect is to 
increase the stream pro tanto. 

If we examine collision (2) we see that the reverse is the 
case. The effect of the collision is to diminish the stream. But 

now, -^ being positive, there are more collisions of the kind 

(1) than of the kind (2) per unit of volume and time. And 
therefore collisions tend on the whole to increase the stream. 

For the same reason, it will be found, they make —rr positive. 

It is true that in a rare medium the effect here indicated is 
very small. But the question is only of its sign, not of its 
magnitude. 

70. Since every collision changes the direction of motion 
of the colliding spheres, the expression in Art. 66 derived from 

(d^ dK\ d(ay) 
{dz'^dxj dt 

contains, in addition to the result calculated in that article, 
a new term derived from collisions which we have now to 

calculate. 

5—2 
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Let V be the half relative velocity of two molecules, X, fi, v 
its direction cosines before collision, V, /i', v after collision. 
Let \F= Vx, vV=Vg, The number of collisions per unit of 
volume and time for given V is 27rc*pF. Hence, for the change 
with the time of VxVg due to collisions with given F, we have 



27rcVF.F»(\V-\i;), 

the bar denoting mean values. We have now to find the 
values of W and \v for all collisions, given F. 

At the instant of collision let the centres of the two colliding 
molecules he A, B; so that AB is the line of centres. Let BD 
be the relative velocity. Or, if DE = BE, one molecule has 
velocity DE, and the other BE, in addition to the velocity of 
their common centre of inertia, whatever that may ba Let 
the angle DBA = ft If -40 = BO, is the point of contact 

Let BX, BZ be the directions of the axes of x and z, and let 




the angle between the planes DBA and DBZ be ^, and the 
angle between the planes DBA and DBX be ^'. 

Then, if x, y, z be the coordinates of A referred to 
as origin, and if the direction DE be that of X, /i, v. 



c ^ . c 



z = — v^ cos + n Vl — i>®sin cos <f>, 



(22). 



X 



= — X 5 COS ^ -f o ^/l ~ X*sin cos ^' 



The complete definition of ^ shall be this : — 
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When the plane DBA, turning round DB, contains the axis 
of 2^, ^ = for that position of BA which makes the least angle 
with the axis of positive z. And ^' has a corresponding defi- 
nition for the axis of x (see Art. 72, post). 

71. If fo> So be the values of f and f at 0, then at A 
f - f = 5 ( - ^ cos 5 + Vl-i;* sin 0GO8 6)^ 

: ,— ft ^ <^>- 

f-f, = |(-.\cos5 + Vr^sin^cos^')^ 

Then by the same reasoning as employed in Art. 67 the 
number per unit of volume of pairs of molecules whose com- 
ponents of half relative velocity are 

r„...ra^+dv^, Vy...Vy'\-dVy, v,...r:,+dVz 

is at il 
2h 




IT 

Here owing to the variation of f in -?, V^ is involved asym- 
metrically with respect to Vy and Vg, Let us suppose now V 
constant, but the angles and (f> to vary. The general mean 
value of VxVg is of course zero. But the mean value of VxVg, 
given V, in the asymmetrical system is to be found by integra- 
ting the above expression according to and ^, and is not zero. 
Let us denote the mean value of VxVg, or \vV\ for all values 
of and if), V being constant, by (VxVg), the mean value of the 

same function when False varies being denoted by V^Vg. We 
have then, treating f alone as varying, 

( r^Vg) = — f *' 2 sin 0co8 0d0 f d<^ (1 + 4A' F^ (f - f o)) 
ff J Jo 

in which f — f © ^^ ^^^ value given in (23), and the general 
means of functions of Vx and Vg, or XF and vV, are taken 
after integration according to and ^. That is, writing Vx/V 
for \ &c. 
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Jo 



+ 

TT 



I ^'^' ^r ^^j'^'f' ^ ^'^ ^ ^^ ^ *^' I ( "" T ^^® ^ 



^v:-^vj . 



H 2_ !L sin ^cos ^ 1 



4A'c F.»F.» d^ 



(24), 



3 F d^ 
because the general mean value of Fa^F^, or \i/F', is zero. 

We have thus proved that, on the average of all colliding 
pairs of spheres, (VxVg) is before collision of the opposite 

sign to j^ or -1^ . This does not depend on the result of 

collision. It is therefore true if, instead of being elastic spheres, 
the molecules be centres of force, or whatever be their form. 

72. We have next to calculate the mean value for all 
collisions of \V, or the value of \v after collision. This will 
depend on the form of the molecule, and we shall now deal only 

with elastic spheres. It will be found that XV is of the same 



sign with (^ + ^). 



We use the same notation as before. When the plane DBA 
turns round DB until it contains the axis of z, let Ag, Ag, which 
are on opposite sides of DB, be the two positions in which that 
plane is cut by the circle which A describes. Let Ag be nearer 
to the positive axis of z than AJ. Then, for 4„ ^ = 0; for 
Agj <f> = 7r. Similarly, when the plane turning round DB con- 
tains the axis of w, A^, A^ are the two positions of A in that 
plane, and, if A^ be nearer than A^' to the positive axis otw, ^' = 
for iljc, ^' = TT for Ag. 

Let € be the angle between the plane of DB and z and the 
plane of DB and a?, so that <f> =<f>—€, <f> = ^' + €. It will be 
seen that, if \v is negative, cos e is positive ; if Xi; is positive, 
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COS € is negative. In either case 

— Xv 



COS€ = 



sin€ = 



A* 



\/l-X«Vl-i;«' 



.(25). 



we have then, to find X' and Vy 

v' = - i; COS 2^ + Vl - v" sin 2^ COS ^, 

X' = -Xcos 20 + Vl-Vsin 20 cos ^' 

Whence, substituting VJ for VF, and F/ for vV, we derive 
the two symmetrical systems 

F; = -FrCos2^ + VF/ + F/sin2^cos^ 



/ 1 



F,F. 



F«' = - F, cos 2^ - ^^^^^ sin 2^ cos 



F/ = - F, cos 2d - 



F«F 



...(26). 



a?»^ « 



- sin 20 cos <f)' 



VF,« + F? 

F„'= - Frcos 2^ + VF/ + F/ sin 2^ cos f j 

The term involving sin e will disappear in the subsequent inte- 
gration, and is omitted. The above systems hold for all values 
of Fa; and Vg. 

73. We now find that for given F, taking into considera- 
ble 
tion only the variation of f , or -p , 

(V^TJ) = i I^d0 rd6 2 sin cos tf (1 + ^h'r^f^,)V^'r;. 
frJo Jo 

But with the above values of V, v\ov -~ , -r^, 



f d^rd(^2sin5cos«F/F; = 0, 

JO </o 



as is easily seen. Therefore 

1 C^ 



(VJV;)~f d0{ #2sin^cosd4A'F«f-foF«'T^; 
^yo Jo 

= -4ih'c[ doTd^ An cos 0V^ 

TT Jo Jo 



x(-;- 



COS0 + 



^v:+v* 



s^sindcos^^JF^F,'^. 
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by (23). In this expression we have to substitute from (26) 
the values of F/ and F/, and integrate according to 6 and ^. 

To effect the integration according to ^, we reject odd 



powers of ^, and for cos^ ^ write ^. We have then 

(F/F/) = - ^4A'c ^'de rd68m0coB0r^ 
' IT az Jo Jo 

X |-^cosd+ ^^^-L+J^sindcos^ r^'V; 

ih'c I^d0 r d6 sin cos 0VJ'QH, 

Jo Jo 

( F. V'F'TF* ) 

where F^ j— ^ cos ^H ^= — ^ sin cos ^j- , 

© = (« y^ cos 2^ + VT7TT7 sin 25 cos <^), 
H^ f- F^cos 25 - --=^^= sin 25 cos 6] , 






or(F«'F/) = -24A'c-?-y-^ sin 5 cos» 5 cos» 25 d5 

+ g 4A'c-!-^ i / sin 5 cos« 5 sin^ 25 d0 

d^ FVF*4-F*^ r*» 
-^*h'c ""^ y^^'^ ^j sin>5cos5sin25cos25d5 

+ ^4A'c -^p-?- i / sin« 5co8 5sin 25 cos 25 d5 

^^'^ih'cUp^j sin5cos>5cos«25d5 

^1: F*F* ri» 

+ ^4A'c-!^* ij sin5cos>5sin«25d5 

- J Wc V.'V^ j^sin^ cos 5 sin 25 cos 25 d0 



+ ^4A'c-yiJj sin« 5 cos 5 sin 25 cos 25 d5 



^^ OF THK ^ 

UNIVERSI r :>! 
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+ ^4A'c-^f^i sin' ^ cos ^ sin 25 cos 25 d^ 
"d^*^^ F 3.5,7 +d5*''''3.5.7*' ^' 



but ^^"-^ ^/^.'- 1 



F» 3' V, 3.5 



or F"F=5. 






V ' 
Hence the expression becomes 

df .., V*r; -11 + 12-4 + 10 
^*Ac— p 37577 

By symmetry if we treat -p as varying instead of -~, we 
obtain the expression 

3.6 da; ^'^'^ V ' 
and therefore, given V, 

And therefore comparing this with (24) 

(r.-F.,-(r,n)=4»-«Z^(|+jL^)(g^§) (27). 

74! To find the rate of change with the time of F^F^ 
given F, we multiply the last expression by 27rc'pF(F being 
the half relative velocity), that is by the number of collisions 
which take place with given F per unit of volume and time. 

Let us denote the change due to collisions by rr . 
The result is, given F, 

|(F.F,)=8*W,(l + 3lj)K.-F.-(f + D. 
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And then we introduce the factor 



v/f 



2A' 2* ^-Av.' e-*(F/+0 dV^V^V,. 

IT IT . 



and integrate for all values of V^y Fy, Vg and so we obtain 

^ 

76. This being the value of -^ VxVg, what is the value 

■ ot 

If Ux, Uy, Ut be the components of the common velocity of 
the two spheres, a, fi, y the velocities of one, a', ^, y those of 
the other, 

fi-^Uy+Vy, P=Uy-Vy, 

Therefore ay + aV = 2 ( ir,[7, + F,F,). 

Hence since Ux, Uy, Z7, do not vary, by collision, . 



^ (07 + aV) = 2 ^ F«F, (28). 



It follows that on average of all collisions 



= 2A '''^^ (3 + 375) U ■*■ J) ' 
and therefore on average 

(S-i)ic^)-A«'K^3^,)(S-iy <^)- 

76. We have next to consider the change in the term 
a7-i^f-i^ + -ir) due to the variation of a, or 7 by collisions. 

As before let us treat f only as varying. Let us use the 
construction of Art. 68. Consider two spheres in collision 
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their point of contact being at G the centre of the a sphere. 
The z ordinates of their centres are 



^ (— i; cos d + V 1 — v* sin 0co&<f>), 



c 



and ^ (— 1/ cos + Vl — i/^sind cos ^). 

To find the value of ^ fp due to the collision, let us refer/ 

to the point P. Its value at the point of contact being /, 
its value at the centres of the two spheres respectively is 

f +^(-v cos + ^1- v" sin cos ^) ^ , 

and /— ^ (— 1/ cos ^ + Vl — v^ sin cos <f>) -4- , 

with corresponding terms in -/- and 4- ^vhich will disappear 
on integration. 

Now for the first sphere a is increased by the collision by 
Vx — ^fl?> for the second it is diminished by the same quantity. 
Therefore as the result of this collision 

9fp = 9 ^= "" ^^ J^("~ " ^"^^ ^ + ^^^^^' ^^ * ^^^ *^)( '^^'^ ^*)- 
For the minus sign here introduced see Art. 68. 

Substitute for VJ — F^ its value fi:om (26), multiply by 
2 sin cos 0, and integrate according to and ^. That gives 

2/af p = - ^ c ^ vl^dol^d^ 2 sin 5 cos ^P(?, 

where -F={- i; cos 5+ Vl — i^* sin dcos^}, 

© = f-\ 2 cos' d - , ^^ sin 25 cos <f)^ . 

That is 2/3fp =f - ^ c ^\i;F. 

Now multiply by ^irt^pVy the number of collisions given V 
per unit of volume and time. That gives 

a£p_ 16 3 rnrl^ 
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Similarly ^ = - j| ■^'^pV.V, ^^ , 

and |(f.-?.) = -^TC'pnF,^(f-f). 

and ljjd-dy(S^-^^) = -^^^^pr.r.jjd.dy^^[%-%) 

where the integrations are throughout the section of the ' a ' 
sphere made by the plane of xy, 8 is the volume of the a sphere, 

and - ^ is the value of - ^ when r = a, 
a da r dr 

But ^jj dady (f p - M = ^ > 

as in Art. 68, therefore 

This expresses the change ^f j^ + j^ due to collision between 

two spheres. Let a, )8, 7 be the velocities of one of the two 
colliding spheres, a, 0, 7' those of the other. 

Then '^:J contains the two terms 
at 

But a = IT, + Fa, a'= IT,- F, 

as in Art. 75. 

Therefore (a7 + a 7') ^x^z^^ VJ F,' on average, 
because UxUzVa-Vg = on average. 
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Therefore for each sphere separately 

because ay = UJU,* + ¥,*¥* + U*V* + Z7/F,* 

= 4iV*V* on average. 

We have then on average, remembering Art. 68 (C), 

Adding this to Cj, we find for the term in -jr due to collisions, 



dt 

dM 
dt 



-§«.,///«,* i (§+§)• (O. 



We interpret this as follows. As the result of collisions, 
molecules moving in a given direction, as that of a?, are on 
Average deflected towards that side where the average motion 
of the molecules, that is the stream, is for the time being in x. 
And so collisions tend to increase the stream. Compare Art. (57). 

77. Now referring to equation (A), we see that ^, which 

in the medium of material points was zero, and so the motion 
stationary, has, when the spheres have finite diameter, become 
positive, namely 



-///«^U|(«4f) 



-^($+S\ 



2A \dz ' dx) 



2 , I (dp df\») 



of which the first two terms are together zero, as in (A), and the 
third is positive. Therefore if in the system of spheres of finit.e 
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diameter ( -jt ) ^- have the same values as in the system of 
material points, -^ is positive. 

But in stationary motion -^ must be zero. Therefore in 

stationary motion of the system of finite spheres (;p) &c., 

cannot have the same mean values which they have in the 
system of material points. 

78. It might be suggested that the explanation of the 
difficulty consists in the fact that when the molecules have 
finite diameters the quantity of momentum transferred across 
any plane per unit of area and time is increased in a certain 
ratio, 1 ; 1 + J)c, where, in the case of elastic spheres, k = fTrc'p. 
But this increases both terms in (A) in the same ratio, and 

therefore, so far as this property is concerned, -r- remains 

zero. This effect arises from the usual convention concerning 
elastic bodies, according to which there is on every collision 
an instantaneous transfer of momentum through a certain space. 
It has nothing to do with the change of direction of the relative 
velocity, which also occurs on collision, and which we now have 
considered. The first effect would take place in precisely the 
same way, if after every collision the direction of the relative 
velocity were restored by Maxwell's corps of demons to what it 
was before collision. 



CHAPTER VI. 

CORRELATION OF VELOCITIES. ELASTIC SPHEREa 

79. The solution of the problem of the last chapter is 
this. In the system of spheres of finite diameter f has become 

That gives by arranging the terms in (D) 

We may assume that *7 j>(j^ + 3^) ^ ^o* altered by the 

substitution of f + f ' for f , or a7 -^ f -p + ^ j = 0. 

The first two terms of this equation are together zero 
by (A) Art. 65. The remaining two terms can be made zero by 
suitably choosing the ratio of the means 

\dz dec) ' \dz dx) 
We have, namely, writing k for f ttc"/?, 

]\dz^ dx) ^\de^ dxj] \dz^ dx) ' 
whence 

if AC*, x* &c. are negligible. 
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80. Instead of using Mxg= I jj dxdydzay(-^ + ^] , or 
the analogous terras Mgey or Myg, we may make use of 

-W«r = jjjdxdydza* ^ , 

or the analogous terms Myy or Mggy and we obtain the same 
result. I have set out the mathematical work in the Appendix 
by reference to Arts. 66, 68. We may therefore now use the 
complete expression for Jf, namely 



Af=///c^.,d.{«.g + ^| + ^§ 



\dy dx) \dz dx) 

-^(g+i)}- « 

and we then reduce -^ to zero by making 

/5ry__^/^Y /dry, ^ ('^^w. 

Now so long as the distribution of velocities is represented 
by €~*^, and Q is the sum of squares of the velocities, ^ &c. and 

(-p ] &c. must have either the same or less values than those 

found by the method of Art. 63, for material points. In 
stationary motion therefore in the system of finite spheres Q 
cannot be the sum of squares only. It must contain products 
of the form 6 (aa' + )S/8' + 77'). And the question now is, what 
must be the form of the coefficients 6, or as I shall call them 

(dW\^ fdi \' 
-p J : f -pi 

may have the required value. 

81. By way of illustration let us suppose n molecules within 
a sphere of radius -K, their velocities being distributed according 
to Maxwell's law e"*^, where Q is the sum of squares of the 



VI.] CORRELATION OF VELOCITIES. ELASTIC SPHERES. 81 

velocities, and let us impart to each molecule, in. addition to the 
velocity which it has in that motion, the small velocity X' in x. 

This being done, the chance that one of these molecules 
shall have velocity represented by lines drawn from the origin 
to some point or other within an element of volume dfl, whose 
centre is distant p from the origin in direction making the angle 
-^ with the axis of a?, is 

Ae-^^ (1 + 2ApZ' cos ^) dn. 

The chance that another molecule shall have velocity repre- 
sented by lines drawn from the origin to some point or other 
with an element of volume dSi' distant q from the origin, in 
direction making the angle ^ with p, is 

A^-^ dXi' {1 + IhjqX' (cos -^ cos <^ + sin i|r sin ^ cos e)}, 

where e is the angle between the plane of p and q and the plane 
of ^ and X, The chance that both molecules shall have the 
velocities aforesaid is 

^ag-A(p«+fl«) dOdO! (1 + 2ApZ' cos 1^) {1 + 2AgZ' (cos i|r cos 

+ sin i|r sin sine)}. 

Multiply this by \ sin -^ d-^, and integrate for all values of 
€ and '^, and we obtain 

^2g-A(pH9«) dftdft' (l + ^'^>3 cos^) , 

expressing now the chance that the velocities p and g, of two 
molecules, in whatever direction either is, shall make the angle 
^ with each other. Now let u^ v, w be the components of p ; 
u\ v\ w' those of q. Then 

dUdSX = dudvdwdu' M dv/ y 
pq cos <f> = uu' + vv* + ^uv/. 

Therefore, since, with these substitutions, the expression does 
not contain <f>, it gives the chance that two molecules, both 
within the R sphere, shall respectively have velocities u...u+du, 
&c., in the form 

B. 6 
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that is, assuming X\ &c., to be small, 

in which 

4A 
o 

« 

Evidently, if we take Y' and Z' into account, we shall have to 
substitute Z'« + F« + Z'^ for Z'« in the last expression. 

We see then that in the case now treated the coefficient of 
correlation b has the form 

Q 

or, since X^+Y^ + Z^ = 



2nA' 



X, F, Z denoting the velocities of the centre of inertia of the n 
molecules, 

b = --(X'^ + T^ + Z^)/(X^+Y'^ + Z') (31). 

82. Next let us employ the converse method, namely, the 

coefficients 6 being supposed given, to find the ratio f'*:f^ 
And first let us consider a finite number of molecules, and all 
the 6 coefficients equal, and very small, so that 6^ 6*, &c., may 
be neglected. 

We have generally, referring to the determinant of Art. 63, 
IP = ^r^ , and therefore 



= r + 9 



2A2)' 



9f ' being the variation of f * due to the introduction of the 6 
coefficients. Now it is proved in the Appendix (n) that in this 
case 9X) = 0. Therefore 

Ai aAi 



f«=a 



ViD 2AD • 
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Also /)„=!+/.•, /^ /^„... 

/1/2, 1 +^*, /a/s, . . . 



and (see Art. 83) 9 (1 +/«) = - 6/,, 3 (l +/,«) = _ ft/^, &c., 

3/i/. = i(6-6/;-6/,),&c. 

If we were now to make every /= 1, we find that every 
constituent in Ai as it originally stood is in the varied form of 

2 — 6 
Ai multiplied by — g— . Hence 

But on the same assumption that every y= 1, 

^ = _L 

And r = aA./2Ai) = -2-l^'^6, 
whence Wi + V^+n= -^VlZ^ i, 



n being large, or 






6=-f(r+v'+n. 



as we obtained in a different way in Art. 81. 

83. Next, let us assume b, the coefficient of correlation 
between the velocities of two molecules, to be a function of the 
distance r between the molecules to which it relates. And let 
it be required to find, on this hypothesis, the ratio 



r--r--m-{i)- 



For the system of finite spheres we shall have 



6—2 
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where bpq is the value of b for the two molecules whose velocities 
are o^, a,, &c., whence, substituting for a„ as in Art. 63, 

-2a,f(l-J6i„2/)-&c. 

+ 2 (/,/, + iftij - iftin/a - ifew/i) «!«» + &C. . . .(32). 

Then f' + fi = ^?5g|^ = f ' + 26 ^'Z 2Ai), 

every 6 being small, so that its square and higher powers may 
be neglected, and, for the same reason, the vaiiation of D may 
be neglected, as will be easily found. Appendix (w). 
Now in this case [see Appendix (o)] 

in which p and q respectively take all values from 1 to n, there 
being n molecules in the system. Also if we assume the 
molecules to be scattered through space without regard to the 
position of p, with density p, SSip^/p/^ may be represented by 
the integral 

4i7rr^pdrbpgfg ; 

and if bpq is very small except for very small values of rpq (see 
Appendix (p)), we find 

tfpr^Trr^pbpJqdr^y;! ^irr^pb^^dr (33), 

J c J c 

= S/p^Si, if 26 = f A^irr^pbdr, 

J c 

and therefore p = ^^ 

because '^fp = Ai \ and if 6 be so chosen that ^26 ^ — /c, 

r=^r (34). 

84. Again, it can be shown that under the same conditions 

(f)='W <^^>' 

and this is the required ratio. 



J c 
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f f 

For since f = «! ^+ Og ^+ &c., 

whence ^""T '^"^^ -«2^-&c., 

and Q = ttj' + a,' + ... + a/ 



20,0, 
0. 



+ -aV' «i«a + &c. 



Now the determinant of this expression differs from that of 
Art. 63 only in the form of the coefficients. Therefore 



\dz) ~2hD~2h^^- 



Again, when the b coefficients are introduced, we have, 
substituting for «» in ^imaiOn, &c., 

«-»T-(S)"-^('^?-'"r:)<+*«- 

and by the same process as before 

U^ / W ~2hD~ 2h~ \dz) ' 
It follows from the above that any form of b which makes 

J2i = — «, or^l 4nrr'pbdr = — K, 
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and satisfies the equation (33), gives a solution of 

and the solution is independent of the form of the function/, 
so long as that function satisfies the conditions of Art. 61. 

The value of 26 is thus determinate. But the particular 
form of b which will make ^ / 4nrr^pbdr = — /c is not determi- 
nate. For instance, let 

3 r** 

That is a solution of (33) and (35) for all values of /i greater 
than 3. 

Again, let 6 = — — when r < a, 

, /A - 3 c» a** , 

6 = — ^— 1 — when r> a. 

fji a^ r*^ 

That also is a solution for all values of /i greater than 3. 
Here a may be the radius of a sphere which on average contains 
one molecule, as in Art. 61. 

85. The law of distribution of velocities being €~*^, where 

let a, /8, 7 be the velocities of a molecule of mass m, a, ^, 7' 
those of any other molecule. Then, 6 being small, it is easy 
to prove that 

36 



= -rs6, 

= 2^:7, 

y^ that is, ^S (aa + fiP' + 77') is equal to the Virial (taken as 
positive) of the intermolecular forces, represented in case of 
elastic spheres by collisions. In Art. 58, dealing with finite 
intermolecular forces, we advanced a step farther, proving name- 
ly that the value of |m (aa' + P^ + 77') for any given r is equal 
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to a particular part of the Virial specially related to that value 
of n If we are to make a similar step in advance in the case of 
elastic spheres, it must be by dividing the Virial, if that be 
possible, into distinct portions, each specially related to a 
particular r. Instead of the Virial of the definite force which 
w! exerts on m at distance r, let us calculate the chance that 
m', so situated, shall be the next sphere to collide with m, 
multiplied by the mean Virial of that collision if it occurs. 
That gives us a part of the whole Virial which stands in 
special relation to r. 

Consider a molecule m. The number of collisions between 
w. and m with relative velocity V in time dt is proportional to 
irc^pVdt, Suppose the whole system to have attributed to it a 
velocity equal and opposite to that of m for the time being, so 
as to reduce m to rest. Then the sphere mf which collides with 
m with relative velocity V must have described a certain free 
path, which must have commenced at some distance r ,..r + dr 
from m. 

The chance that it shall have commenced at distance 

r..,r + dr is of the form qe^^dr, where q is independent 

of r, but a function of F. But mf starting from collision at 

distance r from m, is as likely to be moving in any direction as 

in any other; we must therefore divide e"^ by 47rr*. Again, 

. 2 Vc 
the Virial of the collision if it takes place is h ^ -r by Art. 28. 

Let us then introducing the fiEictor — , assume for given V 

, , irc'F c „ 2 Vc 

and therefore, given V, 

i26 = I inrr^pbdr = ^ tc*/» —j— j — €~^dr 

2 ^ mV* 

and taking means for all values of V 

\U = -lir<^pT=-KT. 
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To express 6 as a function of r, we must perform the 
integration according to V, and not according to r. That is 



1,_2 TTC* 



irc\'2^"mF*€-«^dF, 



2 3 47rr* 
in which 9 is a function of V, and C is the usual constant. 

That appears to satisfy all necessary conditions for 6. But 
it is not a complete solution of the problem, because the intro- 

duction of the factor — ^ is not proved to be necessary, except 

to satisfy (33). It is possible that, as in Art. 57, the corre- 
lation, that is, 6, may depend on the angles which the velocities 
of two molecules make with r, the line joining their positions, 
as well as on r. 



CHAPTER VII. 

FINITE INTERMOLECULAR FORCES. BINARY ENCOUNTERS. 

86. It is worth while to apply the method of Chapter v. to 
molecules between which finite forces act. Nothing is assumed 
concerning the force acting between two molecules, except (1) 
that it acts in the line joining their centres, (2) is a continuous 
function of the distance between their centres, (3) becomes 
repulsive and infinite as that distance is indefinitely diminished, 
(4) is evanescent at all distances greater than a certain line 
c which is very small compared with the dimensions of the 
space in which the system of molecules is moving. A sphere of 
radius c described about a molecule shall be called its sphere of 
action. 

It seems necessary to assume (3) that for sufficiently small 
distances the force becomes repulsive and infinite, if we are to 
hold that there is a very small distance within which two 
molecules cannot approach each other. 

87. With this definition of a molecule, let us assume that 
the number of molecules scattered through a given space is so 
small, or the aggregate volume of their spheres of action bears 
so small a ratio to the space, that no molecule is ever within 
the sphere of action of more than one other molecule at the 
same instant. In other words, the encounters are binary. Or 
that the number of such complex encounters which take place 
in unit of volume and time bears so small a ratio to the 
number of binary encounters, that we may without appreciable 
error assume all the encounters to be binary. 
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88. Suppose first a molecule P fixed in space. Let PX 
be any line through P. 

Let Q be a point just outside of the sphere of action of P 
and distant I from PX, And let another particle at Q of mass 
m move with velocity q parallel to XP, An encounter takes 



Q 




place, and since the encounters are assumed to be binary, 
m describes a symmetrical curved path QRQ, R being the 
apse. Let the angle XPR^yjr. The whole eflFect of the 
encounter, so far as regards the change of direction of the 
motion of m, is the same as if m received an impulse 2mq cos '^ 
in direction PR, Then also P experiences an impulse 2mgrcos'^ 
in direction RP, 

Let the number per unit of volume of molecules which at 
or near Q have velocity q ,,,q -\-dqm direction XP hef{q) dq. 
Then the sum of the impulses which P experiences in unit of 
time from all the molecules which so move is 2m/{q) q^dq cos yfr, 
and the mean force per unit of time on P due to all molecules 
moving from Q in direction XP is 



2m lf(q) (f cos -^ dq. 



'^ being a function of q, and the integration including all values 
of q. If I also varies '^ is a function of q and L 

Now suppose a circle described about X with radius Z in a 
plane at right angles to XP. Let m, moving parallel to XP, 
be before encounter anjrwhere on that circle, and let the angle 
made by the plane QXP with the plane of the paper be e. R 
will evidently lie on a circle parallel to the circle described 
about X. 
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If PZ be drawiL at right angles to PX in the plane of thie 
paper, the mean force on P resolved in direction PZ due to the 
molecules moving from Q in direction parallel to XP is 



2m cos € I f{q) q^ cos yfr sin ^^q. 



And if now I also vary, the mean force on P in direction PZ 
due to all molecules moving in direction XP, whatever the 
value of I, is 

2m I cos edej ^irldl \f{q) q^ cos -^ sin y^dq, 

'^ being now a function of I and q. It is of course zero. 

89. If P, instead of being at rest, be moving with velocity 
u in direction PX, we must write u-\-q ioT qm. the above 
expressions. The force on P for any given value of e will 
then be in direction RP, and equal to 

2m I 2irldl j dqf(q) (q + w)» cos '^. 

And the mean force on P in direction PZ will be 

2m / cos ede I 2irldl \f{q) {q + u)^ cos -^ sin '^dq. 

Now let the plane of the paper be the plane of XZ, XP the 

direction of x. Let the function f have at P the value fo> aiid 

d^ 
at any point whose z ordinate is dz the value ^o + -f^ dz. Then 

the velocity of P relative to Q will he q + u-\-~^l cos e. 

And the mean force on P in direction PZ is now 
4wm I cos ede j Idl I dq/(q) (? + «^ + ^ t^ cos e ) cos yjr sin '^. 

That is, if u, v,w be the component velocities of P, 
— = — 4flrmi coseeJei Wnd5/(g)f j + «^ + i-r^cosej cos'^siu'^ 

= — 47rml WZl (ig'/(g)(gr4-t*) i-j^cos-^sin-^. 
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And u being negative let us write — w for w. Then 

u-jr — + 4s7rmj Idl I dqf(q)(qu + u^)l^cGS^Ir&iayp^ 

= + 47rm / PcU I dqf(q) v? coa '^jr sin yfr -^ , 

m 

on average of all values of i^ ; or writing ^ for nm^ 

dw _ df 1 

where K = 4nr I Ml I dqf(q) cos yfr sin ^fr. 

By symmetry w-ir has the same value miitcUis mutandis, 

and therefore -j-uw^ k (^ + ^ I , 

dt \dz dx) 

(S-g)^--'(S-S)' ••■•«■ 

90. This agrees in form with the result obtained for 
elastic spheres in Art. 76. But /e, instead of being a known 
constant, now depends on the unknown law of force and is a 

function of (f. Further, in the expression for /c, cos-^ is 
necessarily positive, because -^ lies between zero and ^ , and 
also cos ^ sin -^ is necessarily positive, and k is necessarily 

TT 

positive. Now if ^ > ^ , cos ^ sin -^ diminishes as -^ increases ; 

that is as q increases, because -^ is positive ; if -^ < j , cos ^ 
sin '^ increases as q increases. 

Also the condition that '^ shall be less than -j is that I be 

less than a certain magnitude, which we may call io, and the 

greater q is the smaller must l^ be. It follows that as 5*, or 
the mean kinetic energy, increases, it becomes less probable 

TT 

that -^ shall be less than - , and therefore more probable that 
cos '^ sin '^ shall diminish as q increases. 



I 
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Therefore if gr" be great enough, cos -^ sin y^ taken on average 
of all values of I and q, will ultimately diminish as q^ increases. 
And therefore k must for suflBciently high values of 5', that is 
of the mean kinetic energy, or temperature, if the two are 
identical, either diminish as the temperature increases, or in- 
crease less rapidly than the temperature. 

J TLf 

91. We have now to consider the other term in -7- , 

at 

namely, 



07 



dt\dz dxj' 



Let us again assume the encounters to be binary. As in 
Art. 68 describe a sphere of radius a about any point 0, and let 
PP" be a double ordinate parallel to z, P and P' being on the 
surface of the sphere. 



Then g = |J"(irdy(fp-M- 



4tira^ 



Suppose two molecules to encounter one another, the point 
of contact of their spheres of action being at the centre of the 
sphere. Let F^, Fy, F^ be the components of their relative 
velocity before encounter. We may suppose their common centre 
of inertia at rest, since its motion will not alBFect fp on average. 



Let us consider the change of -^ due to this encounter 
only, that is 



//«»(f-t)- 



dp 
First to find -^ , let / be referred to P as centre, and let 

f^ denote its value at the centre of inertia of the two spheres of 
action. Let / be the value of / at the centre of either sphere 
when the encounter begins. After the encounter has ceased / 
will have become/', corresponding to the final relative position 
of the two spheres of action. Let F^., Vg have become V^y F/. 
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If d^p be the change in fp due to the encounter, d^p is, 
omitting the factor ^ which is on average constant, propor- 
tional to, and has the same sign as 

that is, if we write Vx = W, Fa.' = VF, &c., the same sign as 

x/'-V". 

Now let be the angle between the relative velocity, whose 
direction cosines are \, /i, v, and the line of centres at the 
commencement of encounter. 

Then 2-^ — ^ is the corresponding angle at the end of the 
encounter, -^ being the angle so denoted in Art. 88. 

Then we have, ^ being the radius of the sphere of action, 



/=/o+ |(- I/COS ^+ Vl - i^sin ^ cos 0) ^, 
/'=/, + 1 (- 1/ cos 2^-^ + Vl - v^ sin 2>^-^ cos if>) 



df. 



also as shown in Art. 72, 

\' = — Xcos2'>!r — , sin 2'>Ircos0, 

and therefore, omitting cos if> and writing ^ir for cos* <^ as the 
result of integrating according to <^, 



\/' = V/o + 1 Xi/ . (cos 2^ cos 2i/r- ^- isin 2^ sin 2^^ - ^)^ 

= Xy; + |\l/. ^{c0S^(4c0S*'>;r-4c0S»'>;r+l) 

— COS 2 cos* '^ sin* '^ 

+ sin ^ 2 sin yfr cos yfr (2 cos* -^ — 1) 

+ sin sin yfr cos yfr (2 cos* yjr — 1)}. 
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Also X/*= X/o 

+ \ ^ (— 1/ COS ^ + Vi — 1^ sin ^ COS <t>)^ 



= \f^-\v%coB0-f 



dz 
on average. Also for the pair of spheres (X' — X)/© = 0, and 

I 2 cos ^ sin 0d0 (X7 - \f) 
= \r|^[ '2cos»5sin^d^(6 cos*^- 6 cos2'>;r+ 2) 

+ \i;|^f 28in«^co8^c?^3sin>^cos>^(2co8«^-l)...(36), 

Z azj ^ 

which can be easily shown to be positive, being less than yjr, 
and y^ less than ^ir. Therefore ^y 'ji [;r "^ ;r') ^ positive. 

92. If instead of being repulsive as hitherto assumed, the 
force between the molecules P and m be attractive, it is possible 
that the approach of the two may result in a permanent union, 
each molecule revolving in a closed curve round the common 
centre of inertia. Leaving this case out of consideration, the 
form of the path described by m with reference to P regarded 
as fixed will for attractive force be as in this figure. 




TT 



The angle yfr is now between ^ and ir instead of between 

zero and ^ as in the case of the repulsive force. The impulse 

on m is represented as before by 2mq cos y]r, cos yfr being now 
negative. The mean force on P due to a succession of en- 
counters is in direction PR instead of RP, The reasoning of 
Art. 89 will apply equally well to this case, except that in the 
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result u -IT will be of the opposite sign to -^ instead of the 

-1^ + -1- j -1- t/w; will be negative. 

That leads in stationary motion to a positive sign for the 
coefficients 6. It is conceivable that the force between two 
molecules, being a function of the distance r between them, 
may change sign at some value of r. If however there be a 
limit of distance within which the centres of inertia of two 
molecules cannot approach each other, the force must ulti- 
mately, as r is indefinitely diminished, become repulsive. And 
as the density increases the near approaches becoming more 
frequent, the repulsive forces must ultimately predominate. 
The effect of the attractive forces, if such exist, will be, not to 
make the 6*s positive, but to make them have less negative 
value than they otherwise would have. 



93. So far we have assumed the encounters to be binary, 
so that for each encounter either molecule describes a complete 
symmetrical orbit with reference to the other, complete that is 
in the sense that after the two have separated to the very small 
distance c from each other, the motion of either is sensibly 
rectilinear in the asymptote. If we remove the restriction to 
binary encounters, we may reason as follows. 

dp 
Assume -^ to be positive at the origin P. Then a molecule 

at P moving with u negative has greater velocity relative to 
molecules on the side of positive z than to those on the side of 
negative z. On the average therefore the potential of the 
mutual action of P and the other molecules will be greater 
with z positive than with z negative. Let p be that potential. 

Then -^ is positive. 



Therefore if u, v, w h^ the component velocities of P, 

iW dp . , . A J XT. i» ^^ • 

-Y- = — -r- IS on average negative. And therefore u -rr is on 
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average positive, or has the same sign with -^ . Ako its sign 

has no relation to ^ . 

ax 

Similarly w -i- is on average of the same sign as -^ , and 

dp 
independent of the sign of -p . And therefore 

/d^ dt\ d— (d^ df \ / dw du\ 

is necessarily positive. 

And therefore in a system of molecules of the kind now 

considered -j- is positive, and the motion not stationary if 

f", ly^ (p* of Art. 61 have the same mean values as they would 
have in the medium of material points. To effect stationary 
motion, f ", rjj^y f* must be increased as we found to be the case 
with elastic spheres. 



B. 



CHAPTER VIII. 

GENERAL THEORY OF THE STATIONARY MOTION. 

94 Having established the results contained in Chaptera 
IV. and v., I shall now assume for the law of distribution of the 
coordinates and velocities among an infinite number of mole- 
cules in stationary motion the following, namely : The chance 
that at any instant the coordinates of the molecules, n in 
number, shall lie respectively between the limits 

Xi and Xi + dxi, 

^1 » Zi-\rCLZxy 







• • • • 


and 


Xn + dXn, 






Vn 


» 


Vn + dyny 






Zn 


» 


Zn + dZn, 


and their 


component 


velocities between the limits 






«! 


and 


«! + duu 






A 


» 


I3i + dl3,, 






7i 
• • • 


iy 


7i + dyu 



n 

2 



On and On + dOny 
Pn « fin + dfint 



(h\ 
-j ^A€''^^dxi...dzndai.., djn. 
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Here Q = ai(ai« + /8i« + 7i«) + aj(aa» + /3a'+72') + &c. 

+ 612 (aiOa + /SA + 7172) + &c 

The ' a ' coefficients are independent of the positions of the 
molecules. Every 6, as bpg, is a function of the distance, rpg, at 
the instant in question between the molecules to whose 
velocities the suflSxes pq relate, such functions satisfying the 
conditions of Chap. v. D is the determinant of the coefficients 
inQ. 

Again, A may be explicitly a function of the coordinates 

X, y, z, &c., or may be a constant, V5 is a function of the 
coordinates, only as contained in the coefficients 6. Then 
evidently 

n 

Instead of using the word chance, we may say that the 
above expression represents the time during which on the 
average of any very long time the coordinates and component 
velocities will be found respectively within the limits aforesaid. 

96. The function Q may conveniently be put in the form 

+ ttaOa* + iftjaOatti + ib^sfil^ + &C. 
+ &C., 

in which each line is appropriated to a particular molecule. 

96. According to the results obtained in Chapters iv. 

— VI. it is necessary for stationary motion that f '* + rj'^ + f '^ 
(Art. 79) shall have a certain value. That necessary value it 
will have if suitable values be given to the coefficients b as 
functions of r. The coefficients 6 being so determined, the 
motion is stationary, so far as not to be disturbed by the 
collisions or mutual forces of the molecules. I assume now 
that in Q the b coefficients have the values so determined. 

7—2 
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97. We can now deduce the fundamental property of this 
distribution, namely, that for any given set of positions of the 
molecules, or as we shall call it, for any given configuration, 

dQ dQ . 

«! -=— = Ota -f^ = &c., on average. 

For -^ = 2a^ai + hj^ + ftwOs + &c. 

If, with «! constant, Os, a, &c., assume all possible values 
consistent with the conservation of energy, we find (see 
Appendix (A)) that the mean value of 

-^ , or 2aiai + h^d^ + &c., 

is YT ai, D being the determinant of the coeflScients a, 6, &c. 

in Q, and Ai its coaxial minor obtained by striking out the row 
and column containing a^. Therefore also 

dQ D , 
"^^d^^A;^^ on average. 

But it is also shown (Appendix (gr)) that 

> hD' 
whence on average 






Similarly "> T^ ~ r ^ 



and OjT^ = Oj J =&c I. 

dUi OMi 



y 



98. Ina smuch as Jhe syste m is infinite, and there is no direct 
actio n between any molecule a nd any other molecule e xc e p t at 
distances very small compared with the dimensions of the 
system, not only is the sum of the potential and kinetic 
energies constant, but the kinetic energy of the whole system 
is itself separately constant. Let the kinetic energy be 

r=Jmi(ai8+i9i« + 7i«)+iwa(a2*+/82«+72')+&c. 
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Here if «! be velocity of translation m is mass, if a be 
velocity^f rotation m denotes moment of inertia, and so on. 
Then/we have by the constancy of T, 

^1^1 -J7 + ^«2-rr +&C. =0 II. 

at at 

For the same reason Q is separately constant, or 

^^ + |«§ + &c = o III. 

aai at da^ at 

Now II. and III. suggest the inference 

An obvious solution of this is Q = jSm (a^ + /8^ + ff). But, 
as proved in Chapters iv. and v., this does not give stationary 
motion when the molecules have finite dimensions, except in ^ 
the limiting case of infinite rarity. 




99. Concerning the Maoowell-Boltzmann Law 

rOyPL* = T^Oj' = &c. 

This follows at once from Aj^.. 07, and T V.^ if TV. l>e accepted. 
It must however be noted that IV., although a consistent, is not . 
a necessary, consequence of II. and III. For while ai, ota, &c. are yr^ 
unchanged, let all the molecules undergo small displacements. 

Then -^^ , -^ &c., depending as they do on the intermolecular 

forces, will receive variations, which may be arbitrary, as the 
displacements are arbitrary. If under these circumstances 

JQ 

-p , &c. remained unchanged, IV. would be a necessary conse- 

quence of II. and III. But -j^ &c. do not remain unchanged. 

So IV. is not proved by II. and III. 

If we assume condition A, we assume in effect that the 
6 coefiicients are all zero, and Q = 'Zma\ It follows then at 

once from the law of distribution €~*^ that wiitti* = m^a^^ = &Dcr^f^ 
The law is therefore true whenever condition A can legitimately 
be assumed. It is true, that is, in the limiting case of infinite 
rarity . 
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/ I am not aware that it has ever been proved in any other 
case, or by any assumption not equivalent to condition A. 

If it were true in the general case, we should have, comparing 
II. with IIL, 



and therefore by Art. 97, 

D D 



= &c. 



4 



^ 



or whAi = ^WaDffl = &c. 

But Du D22, &c. are functions, not only of the masses of the 
molecules rrii mj, &c., but (if they be elastic spheres) also of 
their diameters, (if they be centres of force) of their eflFective 
voluines. It seems therefore to follow that the law 

mitti* = maOa*, &c. 

cannot hold universally. It can be accepted only on the 
authority of the great physicists by whose name it is known. 

100. If therefore we write 

dQ 



g^/^ = AMX2, &c., 

fii, fi2, &c. will in general have different values for different 
kinds of molecules in our system. Only in the case where the 

molecules are all of the same kind, may we write -^^ I m = fia, 

where fi is constant. 

In all cases however in which the b coeflScients are very 
small /ii, /ia, &c. will be very approximately equal to each other. 
For in the limiting case when the Va are all zero 

Q = Oifli^ + a^^ + &c., 

and we know that Oi = mi, 03 = ^, &c. When the 6*s are not 
zero, every coaxial minor, as Dn, consists of the product of all 
the axial constituents miTn^ . . . , and of other terms each of which 
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contains the product of not less than two &'s. The quantities 
nfhPiit w^aa^ &c. will therefore, as long as the 6*8 are small, differ 
from each other only by small quantities of the second order. 

Field of external force. 

101. We may now equate to zero the time differential of 
A'J D€~^^, and so obtain, applying Boltzmann's general method, 

\dt dx dt dy dt dz) 
T» X dx dy r, dz 

^"* d«=«' i='^' di^'f- 

Again let the forces acting on the system have a potential X' 
Then 

^ dai dx 
dt dxi * 

jn.^ ^ &c 

And the equation becomes, arranging the terms, 

\ dx m dx da J \ dy m dy dp) 

y dz m dzdrfj 
1 ^^ dDdb , s 

■^^U^^^Hbdi''^- 

Let us first consider the terms which involve differentiation 
of the 6's. Since every 6 is a function of the corresponding 
r, and of that alone, 

db ^dbdr ^db x^af « 
dx^ dr dx'' dr r *' 

where r»=« (a; - a/)> + (y - i/y + {z-/y. 
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mt /. K«>c dDdb ^^ dDdbdr 

Therefore S2a jt- -7- = SSa jt zr 1- 

do dx do dr ax 

dD db^ Xi—x^ dD dftu x^ — x^ ^ 
dD dbj2 Ola ^^2 o^ 

1 dTi2 dTi2 Q 

because -r- = ~T~ ®c. 

The above expression can by taking the terms in pairs be 
resolved into a series of terms of the form 

dDdb {a-a'){x''x) 
db dr r ' 

and is therefore zero on average. 
Similarly 

can be resolved into a series of terms of the form 

//yV . zjyQ/ . ^\dh (a'-a)(x'-x) 
(aa' + ^^ +^7)^ y 

and is therefore also zero on average. 
102. We have then 

\ cw? m dx day 

+ &c. = 0, 
or dividing by 4, 

/ dlog^ h dxdQ\ ^ ( dlogA h dxdQ\ ^ ^^^^ 
\ * dxi mi dxi daj \ dx^ m^ dx^ da^J 

And taking mean values of -^ from Art. 97, 

doL 

^ f dlogA ^dx -P y . (dlogA ,dx D \ ^ _q 
^\ dxi dxi miDiJ \ dx^ dx^ m^D^) • — • 
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And since — =r- = — =r- = &c. = /Lt, 

This is satisfied by assuming, if the molecules be all of 
the same kind, 

^ = €-*'*x (37). 

We have then for the law of distribution of coordinates and 
velocities in a field of external force whose potential is x 

n 

(-] V5 e-'^x e-A« d-i?i . . . dznda,. ..dyn (38). 



Now mitti^ = 



hD h/jL 



and nr=i2m(a« + i8^ + 7'^) = 2^, 

the system consisting of n molecules of the same kind. 



But Q = OiOa' + ftMaiOg + &c. 

_ cRi Ai + biJDu + 032)22 + &c . 

hD 

_3n 

Therefore Q = finT (39). 

103. With regard to the inclusion of the intermolecular 
forces under the potential x ^^^ same considerations apply as 
in Chap. ill. I shall assume that we may use the form 

As an example let us now again consider the system of 
equal elastic spheres in vertical column. Arts. 29 and 45. In 
that case x is the potential of the constant vertical force. Let 
us suppose K, and therefore all the b coefiicients, very small. If 
the X velocity of any molecule, as i^, be given, the mean value 
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of Wj is — y~ Ui, where D is the determinant 



2) = 



U^ ^ Dj23 • • • 



for in this case the a coefficients may each be replaced by unity. 

Now Di2 contains the term —612 multiplied by axial coefficients 
only, and other terms each of which contains the product of at 
least two 6 s. Therefore if each 6 be small enough, the mean 
value of u^ is ultimately — ^&ut^. 

Similarly u, = — J h^y^ , &€., 

and therefore j^ = ym,{Uir^ = ( 2 — J St*-*) 2^1 , 

or 7ni/i = (2 - \W). 

It thus appears that e-^^ diflFers from e""**^, when k is small, 
only by small quantities of the second order. Now the fiwjtor 
€~<*"*^, or €~*, in the expression for the density in Art. 29 
corresponds, not to e"*^^, but to e""^*. That is, if '^ be the 
potential of the intermolecular forces, €~* is what €~^ or e"^* 
becomes in the limit when the forces become infinite for the 
distance r^c. 

And K being small the solution agrees with that of Art. 29. 

It does not however agree with that of Art. 29 when k^^ &c. 

cannot be neglected. Which result then are we to accept 

in this case? Evidently the result €~*'*x, because Art. 29 is 

based on the tadt assumption of condition A, and therefore 

holds only for small values of k. To sum up the. discussions of 

this problem. 

3 

If ic be neglected, the law of density is €""^>f, with 2^ = tt • 
If /c^, &c., but not Ky be neglected, the law is €~**x€~*, with 

If «*, &c., are not to be neglected, the law is e***€~*, with 



I 
1 
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Boltzmann's Mini mum function. 

104. The following theorem is not to be compared in 
egance with the theorem founded on condition A, Art. 38, 
but has the advantage that condition A is not assumed. 

Let /(i^ ... w«a?i ... Zn), or fy be any positive function of the 
coordinates Xi,..Zn, and velocities Ui...Un, of a system of 
molecules, and of certain parameters a, 6, &c. Let 3/" be a 
small variation of/, a, 6, &c., being constant. Let us write 

/+a/=yT+?,ora/=gr/ 

Let the variation 3/) or g/, be made subject to the following 
conditions, viz. 



\\...f\-{-qdxi.., dzndui.,.dun 

^ \\'» •fdsci . . . dZfijdui . . . duny 
or writing dsda for the continued product of differentials, 

jj.../qd8d(r = (1). 

Also lr"f^ -{- qCdsda =^ 1 1 , . .fOdsda; 

jL.fqOdsda^O (2), 

where C is any function the constancy of which is prescribed. 

Again, let 

J7 = JJ.../(log/-l)(fod<r, 

H + dH=^jj...flTq {Iog(/lTg) - 1} dsda. 



or 
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We have then 

i5r+ aff = ||.../(log/- 1) dsda 

+ \\...f\ogl-{-qdsda 
+ jj...fqlog/d8da 
+ 1 1 ...fqlog 1 + q dsda 
— Il,..fqd8da 
= H + jj.../^dsda 

+ jj-'/{q'-'^)dsda 

+ jj' . ./3 log/ dsda 
= H+ijj...fq^dsda 

because q^ is to be neglected. 

In order that dH may be necessarily positive, or H may be 
least possible, subject to (1) and (2), it is necessary, and it 

is sufficient that [[ fq log f dsda = 0, that is by (2) that log / 

shall be a function of C, that is that /= ile"*^, where 

lL..fGdsda I 

is constant. ^ 

I 
We may write (^V A^/D for A and Q for G, The actual 

value of H when minimum is 

|log^ + log4+|logi>-| (40). 
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Physical effect of the b coefficients. 

106. From the law of distribution at which we have now 
arrived we can find the mean value of any of the products of 

the form aa', &c. For instance, ajfla = Ag/AZ), D being the 
determinant of the coefficients. And all the Vq being negative, 
Di2 is necessarily positive. If we have only two variables x 
and y, and if 

6 being negative, we find Dy^ = — 6, which is positive. 
So with three variables, a?, y, z^ and 

we find A2 = — Zajftia + 6136281 

which, every 6 being negative, is positive. (See Appendix (6) 
as to the sign of these minors.) 

The conclusion to be drawn from this is that if two mole- 
cules are so near each other that 6 is not negligible, they are 
on average moving in the saine direction. There is at every 
point on average a determinate energy of stream motion. It 

follows further fi'om this uniformity of direction that if a^, ay^ 

fiy, be the mean values of those functions for all molecules 

within any small finite space, (aj8)^ &c. have greater value 
than they would have if the molecules were material points. 

Another consequence of the h coefficients. 

106. But we might also, treating the velocities as con- 
stants, proceed as follows. 

Let for a moment the chance that two molecules distant r 
from each other shall have velocities in x, u .,, v -{■ du and 
u' ...u' -{■ du' be proportional to 

in which h may be 5 — , and n >3. 
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It is convenient here to treat b as positive, and prefix to it 
the minus sign in the index. That being the case, if 6 be not 
negligible, that is if r be small, uu' is more likely to be positive 
than negative. If r be great, and therefore b negligible, it is 
indifferent. It follows by the theory of inverse probabilities 
that, it being given that uu' is positive, r is more likely to be 
small than if nothing were known of the sign of uu\ If on the 
other hand it be given that uu' is negative, r is less likely to be 
small. 

Suppose a molecule at P with velocity u, and an element of 
volume dcD at P', distant r from P. Let the chance that there 
shall be within da) a molecule moving with velocity u' be y^do), 
when PF", or r, is great, and 6 = 0. Then when 6 + 0, the 
chance becomes 

which is greater or less than y^dco, according as uu' is positive 
or negative. 

Further, the chance that there shall be a molecule within 
da) with velocity either u or - u is, since a priori + u' and — u' 
are equally probable, and the system is supposed infinite, 

W being now treated as essentially positive. That expression 
is equal to y^dco if 6^ and higher powers of 6 are neglected. But 
if 6^ &c. be taken into account, it is greater than y^dco. 

107. It follows from this result that under the law of 
distribution of coordinates and velocities at which we have 
arrived in this chapter, if at any instant we calculated for each 
of the molecules the value of S6, or the sum of the 6 coefficients 
connecting the u velocities of that molecule with the velocities 
of each of the others, the average value of 26 would be 
greater than it would be were all the molecules scattered 
through space at random. That signifies that the molecules of 
the entire system would be collected into denser and rarer 
masses — the effect of which evidently is to increase the meaa 
value of 26 for all molecules. Further the ratio f : f* would 
have a greater value for the denser masses than for the rarer 
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ones, or the stream motion would be greater for the denser 
masses than for the rarer ones. So that the system would tend 
mOre and more, with increasing number of molecules in a given 
\ y^pace, to assume the form of a number of denser aggregates, say 
clouds, moving through a comparatively rare medium. 

108. It is not essential to my argument, but I will here 
give what appears to me the best possible definition of density 
at a given point at a given instant, for molecules of finite 
dimensions. 

I would define namely as follows. The density at any 
point P at any instant is — 26, where 26 is the sum of the 
b coefficients, all being negative, of all the molecules of the 
system referred to the point P — that is, b for any molecule is 
a function of the distance of that molecule from P. 

109. It is assumed throughout this investigation that the 
quadratic function Q is positive, that is, that the coefficients 
a, by are such that Q cannot be made negative, whatever values 
we assign to the velocities. It is assumed that the distribution 
of the molecules in space, and therefore the coefficients b, are 
on the average such as would be calculated by the method above 
indicated, the whole number of molecules of the system in a 
given finite space, i.e. p, being known. Should that condition 
ever feil, the motion of the system can no longer be represented 
by the law considered in this chapter. It is reasonable to 
expect that the change of the mathematical formulae will 
coincide with a change in the physical state of the system. 

110. The condition that Q shall be necessarily positive, 
whatever values be assigned to the velocities, is that the deter- 
minant D, of the coefficients in Q and all its coaxial minors must 
be positive. See Appendix (p). As the b coefficients increase in 
absolute magnitude, D, being positive, diminishes. But as the 
density p increases cceteris paribus, the 6 coefficients do generally 
increase in absolute magnitude, and D, with all its coaxial 
minors, generally diminishes. The b coefficients are generally 
also functions of the mean kinetic energy, or temperature, 
as well as of the density, and if the temperature be high 
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enough, probably diminish as that increases (Art. 90). As the 
density increases, or temperature diminishes, or both, the 
system may ultimately arrive at a state, in which D or some of 
its coaxial minors becomes zero or negative, and then the change 
of physical state ensues, whatever its true description be. Now 
to any group of molecules we choose to select, and therefore to 
all the molecules forming one of the denser masses or clouds 
above mentioned, belongs a particular coaxial minor, that 
namely obtained by striking out all the axial constituents of D 
except those which belong to molecules of the group. Gener- 
ally therefore the condition may be expected to fail for some of 
these coaxial minors before it fails for the complete determinant 
D. So that as density and temperature vary continuously, the 
change of state, or liquefaction if such it be, will take place 
partially and gradually, some portions of the system having 
passed into the changed state, while other portions remain in 
the original state. Further, the higher the temperature, at all 
events after a certain point is reached, the greater the conden- 
sation necessary to make D, or any coaxial minor vanish, and 
therefore to produce the physical change in question. 

Such is the process which our analysis leads us to expect. 
Physicists may consider how far it corresponds with what is 
known to take place in gases under condensation, or on what 
(if any) farther hjrpothesis it may be made to correspond 
with it. 



CHAPTER IX. 



OF MOLECULES AS CARRIERS. 



\/lll. The molecules being elastic spheres, and condition A 
being assumed, to find the mean free path for a sphere issuing 
from collision with velocity o). As the mean free path has 
received more than one definition, I define it as follows. Sup- 
pose a sphere m to start from a collision with velocity o), and to 
undergo n successive collisions. After each collision let each of 
the colliding molecules have restored to it the velocity which it 
had before collision. Then our sphere will describe n paths 
Z], 2a ... Zn between its collisions, each with the same velocity w. 
The mean free path for velocity &> is then 

^ n "t" ^ "t" « « ■ "h ^n 

n 
when n becomes indefinitely great. 

If y^ be the velocity before collision of a sphere mf which 
collides with m, E the angle between the directions of motion 




of the two spheres before collision as in the figure, the relative 
velocity of approach is 

iJ= Vft)' + ^a - 2a)^cos JP, 

Let c be the sum of the radii of m and m\ Let/('^) dy^ be 
the number of m' spheres in unit of volume. Then the number 

B. 8 
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of collisions which m undergoes with the mf spheres per unit of 
time is 

iff fW dyfrrrc^R sin EdE=^ Trd^Ni suppose, 
J oJ 

and the mean free path for an m moving with velocity o), if no 

collisions occur except with spheres of the m' class, is \« = — t^^jt . 

If there be two classes rrii and ma, and m can collide with another 
sphere of either class, we shall have a corresponding number of 
collisions per unit of time between m and ma, which we will 
denote by irc^N^* Then tlje mean free path for velocity o) is 

^ = — rTTr irT and so on. To calculate the mean free path 

is to calculate Ni + N^y &c. 

112. Another method of obtaining this result is as follows. 

Let S be the whole space in which our spheres are moving. 
Consider a sphere m moving with velocity o), and another sphere 
TOi whose velocity and direction of motion are comprised re- 
spectively within the limits -^ ... -^ + d-^ and E ,.,E + dE, as 
in the figure. 

Now mi may be anjrwhere within S. If its centre be at this 
instant within the cylinder whose base is a circular area of centre 
m and radius c, the sum of the radii of m and m^, at right 
angles to iJ, and whose height is Rdt, then, but not otherwise, 
a collision will occur between m and mj within the time dt after 
this instant. 

The chance that such collision shall not occur is then 

ird'Rdt 



1- 



S 



Now the number per unit of volume of spheres mi whose 
velocity and direction of motion are comprised within the limits 
aforesaid is f(y^)dy^^sin.EdE. And therefore the number 
within 8 is Sf(y^) dy^ J sin EdE. Let us now assume that the 
aggregate volume of all the mi spheres within S is negligible 
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compared with S, Then the chance that during dt m shall not 
collide with any one of the class rui is 

raised to the power of 

8f{ylr)dylr^8mEdE, 

that is, in the limit when dt becomes infinitely small, 

1 - irc^Rdtf(it) dyfr J sin EdE, 

and the chance that m shall not during dt collide with any mj, 
whatever be the values of y^ or E, is by the same reasoning 

1 - irc^dt I jfiir) dyft^R sin EdE = 1 - ir<^dtN^ suppose. 

Similarly if there be other classes of spheres, w,, wij &c., 
with which m may collide, the chance that it shall undergo no 
collision in the time dt is 

1 - ir<^dt (Ni + i^2 + &c.) = 1 - irc^dtN, 

By the same reasoning the chance that m shall undergo no 
collision in n successive intervals of time dt is 

1 - irc^Nndt^ 1 -irc^Nt, 
if ^ = ndt 

Also if X be the space described in time*^ 

X = 0)^, or ^ = - . 

Therefore the chance that m shall survive for a distance X 
without undergoing any collision is 1 — irc^N - . 

Let l-7rc^JV- = ^(X) = ^. 

Then ^ is the proportion of the whole number of spheres m 
which, starting with velocity «, traverse the space X without 
collision. Also ^ (X + dK) = ^ (X) ^ (c?X). It follows that 

_r=_<^ =-% if &= , 

aX ^ tt) ^ ft) 

8—2 
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and therefore <f> = C€-*\ Also </> = 1 if X = 0. Therefore (7=1 
and <l> = 6""*\ Whence if X be the mean value of X given o> 



-/: 



6-*^XdX=l= « 



which is the same result as we obtained above (111). 
The expression 

TTC" I d^ / sin EdEf{-^) d^ ^lo>^ + -^P - 2ft)'</r cos E 

Jo Jo 

can be calculated numerically in terms of-Trc^ and h if 

The fundamental equation. 

113. Let be any quality, e.g. mass, colour, charge of 
electricity, momentum in given direction, &c. which a molecule 
may possess in greater or less degree, and which it can carry 
with it unchanged by its own free motion, and which, as 
regards the aggregate possessed by two colliding molecules 
is not changed by collision. Suppose a quantity of gas con- 
tained in a vertical cylinder of which the height is very small 
compared with the diameter of the base. Every molecule that 
enters the cylinder through its upper flat surface shall enter 
charged with the quantity Gi on average of our supposed quality 
Q. Every molecule that enters the cylinder through the base 
shall enter charged on average with the quantity 6^. When 
the motion has become stationary, every molecule within the 
cylinder at a height z above the base will be on average charged 
with a quantity of the thing in question between (tq and (ti. 
Let it be called G{z). G{z) in stationary motion will be a 
function of z. 

114. Now consider a layer Pi of molecules in the cylinder 
between the planes Zi and z-^ + dzi. If the molecules within 
that layer were enclosed between two elastic planes Zi and 
z^ + dzi for a finite time, the distribution of velocities among 
the molecules within it would become, under the circumstances 
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supposed, completely independent of the quantity of which 
they happened to carry. It would become in the strictest sense 
of the term " molecular lingeordnet." In other words condition 
A would prevail. 

Collisions would occur within that layer, and any molecule 
emerging from such collision would be as likely to be moving in 
any one direction as in any other. The chance of its having 
any particular direction would be independent of its charge G. 

Consequently if any such molecule, having undergone col- 
lision within the layer Pi arrives without further collision at 
the layer P, the mean quantity of which it brings with it is 
the mean quantity of per molecule for all the molecules 
within the layer at Pi. This is the fundamental assumption. 

116. Now let Gz be the mean value of that quantity for 
all the molecules within the layer z ... z + dz a,t P, Suppose a 
molecule to arrive at that layer from above with velocity 
between v and v + dv, and in direction making angle between 
and 0+ dO with the vertical. Let X be the distance which 
on its arrival at P it has traversed since its last collision. Then 
such last collision must have occurred at a height \ cos above 
P. Therefore the quantity of G with which the molecule is 
charged is the mean quantity of G for the height z + \cos0, 
that is 

G {z -{-Xcose)^ Gz + \cos0^ , 

because X cos is small. 

The number of molecules which, having velocity and direction 
between those limits, pass through the horizontal plane P per 
unit of area and time iBvf(v)^ sin cos 0d0dv. 

Hence the quantity of G cajried through the layer P per 
unit of area and time by molecules coming from above is 

w w 

vf(v) dvij sin cos 0G {z)d0-\-il X cos» sin 0d0 ^/(«^) vdv 

=^{iO(z) + iX^)f(v)vdv, 
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and taking the mean of all values of \ and t; it is 
Here X = / f(v) vXdv. 



= I /(f) 

Jo 



By the same reasoning the quantity of G carried through 
the layer P per unit of area and time by molecules coming 

from below is J(? (^r) — J \ -y- . 

And the quantity passing through P downwards exceeds 
the quantity passing through P upwards per unit of area and 

time by JX -T- . 

116. This is the fundamental equation from which Boltz- 
mann calculates the rate of diffusion, viscosity &c. in gases as- 
suming condition A to prevail. It would represent with complete 
accuracy the initial motion of the system, if the quality G were 
suddenly attributed to the molecules in the degree proper to z, 
but without regard to their velocities. 

It takes no account, as he points out, of any variation in the 
distribution of velocities which may be caused by the quality 
itself, for if for instance G be momentum in given direction 
at right angles to the axis of our cylinder, such a change in 
distribution must exist. It is assumed that the effect of this 
change in distribution on the mean free path X is so small as to 
be negligible, a condition which must be secured if the given 
terminal conditions, e.g. Go and Gi, do not differ very widely 
per unit of distance. 

117. The method also takes no account of the stream. The 
effect is to set up a stream, so to speak, of G from the upper to 
the lower parts of the cylinder. And that, as will be found, 
causes a small error in the result. The error may indeed be 
negligible so long as the stream velocity is very small compared 
with molecular velocities, which will be the case if the terminal 
conditions do not differ widely per unit of distance. 
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It therefore does not sensibly affect the results calculated 
from Boltzmann's formula in any case to which as I understand 
he proposes to apply that formula. It is worth while however 
to consider the sign and order of magnitude of this deviation 
from complete accuracy. 

We are asked to suppose that in stationary motion molecules 
carrying positive (?, have a mean velocity u downwards ; mole- 
cules carrying negative have a mean velocity u\ u and u 
being functions of z. Hydrodynamical considerations make Ou 
and Ou' sensibly constant throughout the cylinder in stationary 
motion. 

But when collisions take place in the layer Pi, condition A 
is supposed to prevail so that the molecules issuing from 
collision are as likely to be moving in any one direction as in 
any other, and to carry with them, in whatever direction they 
do move, the quantity of due to the height z. 

In fact molecules with positive G, having before collision 
the mean velocity u downwards, will issue from collision with 
some mean velocity downwards, less indeed than u but com- 
parable with it, and therefore will not be moving in all directions 
indifferently, but on the average downwards. And molecules 
moving downwards will possess on average a greater quantity 
of than that due to the layer ia which their last collision 
took place. For let o> be the absolute velocity of a molecule 
m before collision. Let y^ be the velocity before collision 
of the other molecule m concerned, E the angle between 
their directions as in the figure (Art. 111). Then the velocity of 
m after collision is the resultant of (I) the common velocity of 
m and m', (2) their relative velocity, in whatever direction it be, 
after collision. But for the relative velocity after collision all 
directions are equally probable. The relative velocity therefore 
contributes on average nothing to the resultant. 

The resultant velocity after collision is then on average 
the common velocity. And the velocity after collision re- 
solved in its direction before collision is the common velo- 
city resolved in direction «, that is ii m = mf co + y^ cos E. 
We have then for the mean value of the velocity after 
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collision resolved in the direction before collision the ex- 
pression 

A I dV^e-**" '^ \ Va)» + ^ - 2®^ cos E (o) -4- >^ cos E) sin EdE, 

Jo Jo 

That can be calculated as a function of a>. As however we 
intend to treat it as negligible, it is not worth while to calculate 
it here. But we can easily prove that it is necessarily positive. 
We may therefore conclude that the mean velocity after col- 
lision of a molecule carrying positive Q resolved in its direction 
before collision is/® on average, where/is some positive quantity 
less than unity. 

It follows that the mean velocity downwards of the mole- 
cules carrying positive O, which before collision was u, is after 
collision fu. 

We have thus proved that the molecules carrying positive 
issue from collision in the layer Pi with mean velocity fu 
down the cylinder. For the same reason the molecules carrying 
negative issue from collision with mean velocity /tt© upwards. 
Therefore of the molecules which reach the layer P, having had 
their last collision in the layer Pj, the proportion which carry 
positive is not the same as the proportion of all the molecules 
within the layer Pj which carry positive G, as it was assumed 
to be in forming the fundamental equation, but is greater. 

118. The effect of the irdroduction of the b coefficients in 
this case. 

It has been shown in Chapters iv. — ^vi. that molecules near to 
one another are on average moving in the same direction. The 
energy of their relative motion is therefore less than it would 
be were condition A assumed, and therefore their encounters 
pro tanto less frequent. That tends to increase the mean free 
path. 

On the other hand under assumption B the molecules tend 
to become aggregated into denser and rarer masses. That 
tends to increase the frequency of their encounters and there- 
fore to diminish the mean free path. 
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It is probable that for ordinary gases under ordinary con- 
ditions the b coefficients will be very small. And if so the 
effect of the b coefficients on the phenomena of diffusion &c. 
will be of the same order of magnitude as those sources of error 
mentioned in Arts. 116, 117, which as we have said Boltzmann 
and others have agreed to neglect in calculations concerning 
these phenomena. I shall not therefore here further consider 
the effect of the use of assumption B instead of A on such 
calculations. 



CHAPTER X. 

ON THE HYPOTHESIS THAT HEAT OR TEMPERATURE IS R EPRE- 
SE NTED BY THE KINETIC ENERGY OF MOLECULAR MOTI ON. 

Natanson^s Theorem. 

/ 119. M. Ladislas Natanson (Interpretation cin^tique de la 
foQction de dissipation — Bvlletin de VAcadSmie dee Sciences de 
Gracovie, D6cembre 1893) considers a medium composed of 
molecules in motion. The components of the molecular velocity 
of a molecule are u, v, w, those of the apparent or stream 
velocity of an element of volume are f , i;, f. I have here inter- 
changed Natanson's symbols to make them agree as far as 
possible with my own, and so avoid confusion. Evidently, taking 

mean values, 

u = v = w = 0. 

Let K = ijjjpi^^-{-rf^+^)dxdydz, 

p denoting density, and 

E = \ j j j p (u"" -{■ v^ + tu^) dxdydz. 

The integrations are throughout all space occupied by the 
system. 

120. M. Natanson gives no definition of f , i;, f> but is, as I 
understand, willing to accept Boltzmann*s definition — that f at 
-^ P is the momentum in x of the centre of inertia of all the 
molecules contained in a sphere of radius r described about P 
as centre divided by the volume of that sphere, in the limit 
when r becomes infinitely small. 

I should propose myself to define it as in Art. 61 by the 
function /, so that f , i;, f, Natanson's components of stream 



/■ 
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velocity, agree nearly with the f, 17, f, of Chapter v. But 
M. NatansoD does not on his own assumption require any 
further definition. 

The normal pressures per unit of surface are 
The tangential pressures are 

Pxy^pUV, Pxz^pUW, Pyz^pVW. 

121. He then gives the following fundamental equation, 
viz. q being any function of i^ 4- f , v + 1;, w 4- f, the time varia- 
tion of q is found from 

Here X, F, Z are the components of external force, and -^ 

da ^ 

the change of q due to encounters, -^ the total change of q 

with the time. 

Writing i^ + f for q, he finds 

And again, putting g = (u + f )' + (v + rj/f + (m; + ?)*, he obtains, 
neglecting terms of the third order, 

J 

/>^ (f ' + 1;'* + S^ + w' + 1^ + w') 

d — _ 

+ -J- (2^pu^ + 2rfpuv + 2f/» ww) 

+ J- (2^fyuv + 2rfpv^ + 25y> tw;) 

J 

+ ^ (2^fmw + 2i7/»tw; + 2^pw^) 

= 2p(fZ+i;F+?^) 

+/»^(?'+'?"+r'+M''+»'+w) (5)- 
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Combining (4) and (5) we get 

^ 

= Pg^(P + ';' + ?' + «» + t^ + Mi») (6). 

"-' i-«- I-'- §-'■ 

Add to the first member of (6) 

which is zero. Also write 

d , ,- - — 
^^ ip («' + »* + M^) 

And then in the equation so obtained integrate for x, y^ z, 
and make l^ + mri + nl^— at the bounding surface, I, m, n 
being direction cosines of the normal to that surface. Assume 
further that X = F= ^ = 0, or no external forces act. 

The result is 

dE 
dt ■*" 



/ / / (p^*^ + P^^ + P^c + fyowA + puwB + puvG) dxdydz = 0. 



Treating equation (5) in the same way, we obtain 
-jT "■ 1 1 kpw'flt + /)t^6 + fyuj^c + pvw-^ + puwB + puvG) dadydz — 0. 

And therefore, there being no external forces, -j— = — -i- , 

or the kinetic energy of stream motion can increase or diminish 
only at the expense of molecular energy. 
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122. Let us now write 

F^{p- pvJ") a + {p-p^)b + (p- pvf) c 
— pvwA — puwB — puvC. 

By a further use of the fundamental equation, putting 
q=^(u + ^y, and then q^v + rj x w + ?, we find 

+ 2 (pvviy + 2 (pmuy + 2 (pwvy}. 

And finally 

dE 
dt 



^'^^fff(F-p0)da:dydz, 



where = a + b + c = -^ + ^ + ~. 

cut ay az 

123. I t appears then that E, the molecular or heat e ner gy , 
g ains or loses at the expense of the energy of stream moti on, 
according to the law of encounters " loi des chocs mol^culaires." 
If, that is, as the result of encounters the tangential pressures 
are increasing, E diminishes, and K increases. Now this is 
precisely what happens according to Chap. V. ante, when, 
for instance, instead of being material points, the molecules 
become elastic spheres of finite diameter. As the result of 
collisions in that case |* in my notation becomes |* + f , &c. 
And as we have seen, Art. 105, (a7)',(a^)«, and (fiyY increase, 
but only up to a certain point, at which namely the ratio f 'VP 
has the proper value for stationary motion. The relation 
between the energy of visible motion K and that of molecular 
motion £ is in Natanson's theorem the same as that between 
K and heat in Thermodynamics. 

Messrs Bryan and Boltzmann's method. 

124. In order to confirm the hypothesis that the temperature 
of a gas or other substance is proportional to the mean kinetic 
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energy of its molecules, Messrs Bryan and Boltzmann (Vienna 
Sitzungsberickte math, naturw. Classe, Band III., Abtheil II. 
/Dec. 1894, also referred to by Boltzmann, Oastheorie, p. 136, et 
seq.) treat the following case. X and Y are two infinite parallel 
planes. In the space 8i to the left of X is a gas A. In the 
space 82 to the right of Y another gas J8. In the space 8 
between X and Y there acts on molecules of gas A a force 
directed from Y whose potential becomes infinite on Y, and 
there acts on molecules of gas B a force directed firom X whose 
potential becomes infinite on X. These conditions insure that 
no molecule of gas A will be found in 82, and no molecule of gas 
B will be found in 81, But in the space 8 there will be a 
mixture of A and B, and encounters will take place between 
their respective molecules. 

Messrs Bryan and Boltzmann use generalised coordinates, 
but it will be sufiicient in this notice to deal only with velocities 
of translation, as though both A and B molecules were elastic 
spheres. 

It is assumed that the mean kinetic energy of A molecules 
at any point within 8 is the same as their mean kinetic energy 
within 81, and that the mean kinetic energy of B molecules 
within 8 is the same as their mean kinetic energy within 8^. 

It is assumed further that as the result of encounters taking 
place within S, between A and B, the molecules of A acquire, if 
they have it not already, the same mean kinetic energy as those 
of B, That result seems to me to be questionable or to be only 
^^'VB.pproximately true, if the molecules have finite diameters or 
spheres of action. But I state this with difiidence when such 
authorities are against me. 

Making however these two assumptions, the conclusion 
follows that the mean kinetic energy for the A molecules in 81 
must become equal to that of the B molecules in 82* And 
therefore that if the mean kinetic energy be originally different 
for A and for B, a finite quantity of energy will pass from 81 to 
82, or vice versa. 

The condition that no energy shall pass is that the mean 
kinetic energy shall be the same for A as for B. But we may 
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regard the two gases as two bodies in contact, which condition 
is more nearly realised as ^S is diminished indefinitely. The 
known experimental condition that no energy shall pass between 
Ufio bodies in contact is that both shall be at the same tempera- 

I Xure. Hence in the case of the two gases at all events kinetic 

^ energy plays the same part as temperature. 

P rofessor J. J, Thomson's method. 

126. The same hypothesis is tested in a different way by 
Professor J. J. Thomson (Application of Dynamics to Physics and S/ 
Chemistry, p. 91). He puts the following case. Suppose a tube in 
which is moveable a piston of very great mass M, On one side 
of it is a set A of material particles each of mass m^ in motion. 
On the other side a set B of material particles each of mass m^ in 
motion. Considering the two systems of particles as two sub- 
stances, and the piston as a conductor of heat, the condition 
that no heat shall pass through the piston from one substance 
to the other is that the two substtuaces shall be at the same 
temperature. If temperature be represented by kinetic energy 
of translation, then the condition that neither set of particles 
shall on average either lose or gain energy by collisions with 
the piston ought to be that the mean kinetic energy of transla- 
tion is the same for the A as for the B set. 

Let U be the velocity of the piston, v^ the velocity 
normal to the surface of the piston of a particle mi before 
impact. It is assumed that particles striking the piston re- 
bound from it as perfectly elastic bodiea Let U\ Vrl denote , 
the values which CT, v^ assume after impact. Then we have / 

m^Ui + MU= mii^i' + MU\ 

From which deduce 

. „ ,, ^M^m^U^ - 4if7?ii (if- mO Uu^ - 4ifmiV / x 
m,(i^>-u,«)^ ^M+m^y ^""^^ 

Now add together the equations of this form for all collisions 
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which take place in unit of time. Let Ni be the number of 
such collisions. Also let = ^MU^ be the kinetic energy of the 
piston, and 0^ the mean kinetic energy of the particles of the 
set A which strike the piston. As I understand Professor 
Thomson, the mean is taken for all collisions, not for all 
particles. 

It is assumed that the velocity U of the piston is zero 
on average of time, because the piston is being struck on the 
other side by particles of the B set and the pressure is supposed 
to be the same on either side. Hence we may neglect in the 
expression (a) the term which contains U in the first degree. 

The expression then takes the form 

Nfie.^NU20-m(^^ (1). 

Note. This should, I think, be §^i instead of \B^ because as MU^=^ 23^ 
so ^it^'^l^i, Si denoting the whole mean kinetic energy of all molecules 
which collide with the piston or with other portions of the elastic surface 
bounding the A molecules, i.e. 

2^1=^1 (V + v + V)- 

But this does not affect Professor Thomson's argument. 

Similarly for the B set we have 

NM = N,i20-m^^!^ (2). 

and since 3^ = — Nid0i — Njd02 = on average 



a 



This gives 26-^0, = ^i (0, - 0,). 
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Substituting in (1) and (2) 

and therefore ' d0i = dO^ = 0, if 0^ = 02. 

It follows that if 0^ > 0i the eflFect of all the collisions is 
to increase 0i at the expense of 02, and vice versa. And the 
condition that no energy shall pass from A to B or from JB to -4 
is that 02 = 01, or the mean kinetic energy of translation for all ^^fr 
par ticles colliding with the piston is the same for A as for B. 

126. Now let/(i^) du be the number per unit of volume of 
particles of the A set whose velocities normal to the piston lie 
between u and u + du. Then the number of such particles 

which strike the piston per unit of area and time is / f(u) vdu. 

Jo 

The sum of the kinetic energies of their motion normal to the 

piston is ^mi I f(u)u^du. The mean kinetic energy of motion 

Jo 

normal to the piston for all A particles irrespective of their 

striking the piston is ^i / f(u) u^du. Let f{u) be so chosen 

Jo 
that the ratio 



I f{u) u^du I f(u) u^ du 
Jo , ./o 

/•CD • fOO > 

I f(u)udu I f(u)du 

Jo Jo 

is an absolute constant, and therefore the same for A as for B, 
That being the case, when 0i = 02, or the mean kinetic energy 
per collision is the same for A as for B, the mean kinetic energy 
per particle is also the same for A as for B, And therefore the 
condition that on the whole there shall be no transfer of energy 
through the piston from -4. to £ or vite versa is that the mean 
kinetic energy per particle is the same for A as for B, 



B. 



^ 
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ISlt. If we assume for A particles 

and for B particles f(ug) = C^-^'^^\ 

we satisfy the condition . 

Po ssibly there may be other forms of f(u) which satisfy it. 
But however this may be, some assumption with regard to the 
relative frequency of diflferent values otu appears to be necessary. 

Further, f(u) being properly chosen, we satisfy at the same 
time the two conditions (1) that the pressure shall be the same 
on either side of the piston, the number of particles per unit of 
volume being the same for A as for JB, and (2) that there shall 
be no transfer of energy. 

C The law of distribution of velocities thus appears to be 

JL of the essence of the matter, if the molecular motion is to 
y represent heat. 

The Second Law of Thermodynamics. 

128. Our gas system being in stationary motion according 
to the laws above investigated, we might impart to it a certain 
quantity of heat SE, That being done, it would assume a new 
form of stationary motion, in which the mean kinetic energy T, 
the volume ©, and the mean potential energy ;^, will be gener- 
ally diflFerent from what they were in the original motion. If 
P be the external pressure which in the original motion was 
just suflScient to prevent expansion, the change of volume da 
can be eflfected only by doing an amount of work Pd(o against 
the external force P. A part of the energy dE suppUed to the 
system is spent in doing this work. The remainder is spent in 
increasing either (1) the mean kinetic energy T, or (2) the 

mean potential energy, ^^ of the intermolecular forces, so that 

dE = XdT + dx + P3», X being the number of moleculea 

129. If the external forces P also have a potential, there is 
in general a determinate relation between dm and dT, and a 
determinate relation between dta and 3^. So that in this case 
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there is only one independent variable, and no question can 

arise whether or not -^ is a complete differentiaL Any theorem 

proved on the hypothesis of P having a potential is not the 
Second law of Thermodynamics as usually understood. 

But the external force P is in Professor Thomson's language 
a "controllable" force. It is possible therefore by suitably 
varying P to maintain the system in stationary motion with 
any arbitrarily assigned values of T and to. But dT and 
9© being given, 9^ is determinate. There are then, P being 
controllable, two independent variables. It is convenient to 
use P and &> for independent variables. 

130. If in rectangular coordinates the abscissa represents 
6), and the ordinate represents P, any point in the diagram 
represents a determinate state of the system, in which it is in 
stationary motion with the values of P and co corresponding to 
the point, it being understood that the controllable force P is to 
have such value given to it as will be necessary to maintain 
that motion stationary. So also any curve drawn on the plane 
of the diagram represents a series of states through which the 
system might be made to pass successively, P being made to 
vary as required. But it is understood that the change of 
state, that is the passage from one point in the curve to 
another, is effected so slowly that stationary motion is always 
attained, and therefore the kinetic energy corresponding to the 
controllable coordinate o) is for our purpose to be taken as zero. 

If the curve be a closed curve, the system having passed 
through all the states denoted by the curve, has at the end 
of the process the same values of P and 6>, and is therefore 
for all experimental purposes in the same condition as it was 
at the beginning. Contemplating the system itself alone, 
nothing that we can observe by experiment has happened to it. 
Now if denote the absolute temperature, the second law of 

Thermodynamics asserts that in this complete cycle I — = 0, 

or -^ is a complete differential of some function of and «. 
u 

9—2 
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That function was called by some writers the Thermodynamic 
function, by others the Entropy of the system. 

dE 
131. As I have said, -^ being a complete differential, no 

change takes place in the system itself as the result of any 

complete cycle of operations of the kind described. But IdE^O, 

and jPdco^O, Generally a quantity of heat denoted by jdE, 

drawn from some external source, has been converted into 

mechanical work jPdca outside of the system, and a further 

quantity has been transferred from a hotter to a colder external 
body, as shown in treatises on Thermodynamics. But if our 
system be enclosed in a non-conducting envelope, so that no 
heat or other form of energy is allowed to pass into it from 
without, or out of it into external space, it may by suitably 
varying the controllable P be made to expand or be compressed, 
and in so doing it does external work, positive in case of 
expansion, negative in case of compression. This work is done 
at the expense either of the absolute temperature or of the 
potential x ^^ *^^ system itself The line traced on the 
diagram in such a process, a line that is for which 3^=0, is 
called an adiabatic line. If the cycle be A BCD on the diagram, 
AB denoting expansion at constant absolute temperature 0, BG 
adiabatic, CD compression at constant absolute temperature ff, 
and DA adiabatic, it is found that, the heat which must be 
supplied from without to effect the expansion AB is to the 
heat withdrawn during the compression CD, both measured in 
mechanical units, as : 0\ If temperature is represented by 

the kinetic energy, we ought to find T : T' :: : ff, or -^ \^ 

/a^omplete differential. 
132. In the Philosophical Magazine, January, 1876, a 
proof of this proposition was given by the present writer. 
This was employed afterwards in an improved form by 
Dr Watson in his Kinetic Theory of Oa^es, and has been 
accepted by Professor Bryan. It was founded, as will be seen. 
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on Boltzmann's theorem, that the chance of a group of \ 
molecules having their coordinates between the limits 



Zx ... ZK + dZx, 

(in which position their potential is ;^) is proportional to 

€~^^dxidyi ... dzx» 

This is established only on the assumption of condition A. 
I will here give that proof of the second law in order that we 
may subsequently see what modifications are required in it 
when we abandon condition A, and employ the more general 
method of this work. I here use h instead of 2h of Art. 44. 

We have dE=^'KdT + dx+ Pdfo. 

3 
Let T=^— the mean kinetic energy of a molecule. 

1 2A 
Ihen T'^'S' 

and jr"^^3~^ "^"3 ^ "3 ^^' 

Now substitute for Pdfo from Clausius' equation, Art. 24, 

f Pft) = Sim (w^ + 1;^ + 1(;») + iSS iS^, 

j\ j\ 

or Pda>=^i^m(u^-\'V' + vf')— + iX^Rr — . 

Now if there be a general cubical expansion, « becoming 
CD + 3a), every line r in the system becomes 

r + i-r, 

or ^or = i — r. 

^ ft) 

Therefore the last equation becomes 

P3ft) = iSm (t^« + i;« + !(;») — + 22 i»r. 
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Again, for each configuration 

da 



22i2ar=-ax = -^aco. 



Hence taking mean values 

Here -^ is the mean for all configurations of -^ , that is 

the mean of a variation, while d% is the variation of a mean. 

The distinction between dx ai^d ^ 9<» is of the essence of the 
matter. 
Again, 

ltm(u^ + i^ + v^)— = T — , 

^ ^ ^ (O h €0 

if \ be the number of molecules of the system. 
We have then 

and -s--r9<!D = *X ttj 9®- 

3 ^ 6> 3 {2o> 

And therefore 

r"^3"^^^3-^x+t^^-^di^« <">• 

133. So far our equation is derived from physical con- 
siderations. 

Now consider the auxiliary function 

w = log / / / e~*x (foji. . .dzK . 

Then 9t^ = 37 3A + -j— 9g), 
and jT 3A = — 91 1 1 e'^^x^i ••• ^^a 

= -x3A, 



X.] THERMODYNAMICAL RELATIONS. 135 

because, since the chance of the configuration ^i ... ^x is 
evidently 

Again, -j-da) = — 777 /// e'^^h^ dx^.,. dzx 



+ 






the last term being necessary because every element of length 

da?i... is increased in the ratio 1:1 + j^ — by the general 
cubical expansion above mentioned. 

That is J— 9g) = — A j^3ai + - 9ai, 

because diday. . . . dzf^ — \ — , 

ft) 

and 9i^ = — y9A — h^-\ , 

^ dft) ft) ' 



or S>^^--SA^aft) = §91^ + 1x9* (III). 

Substitute this in (II), and we have 

dE 

which is a complete differential. 

134. Let us now consider what modification has to be made 
in the above proof when Q, instead of being the sum of squares 
of the velocities, has the value given to it in this treatise. 

With the form of Q containing the h coefficients, the 
relative frequency of the several configurations is different 
from what it was when Q was a sum of squares only, but the 
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values of x ^^^ -j are not altered for any given configuration. 

Hence x ^^^ ;j have diflFerent values from those which they 

had in the former case, x ^^^ j^ *"*^ ^he same as in the former 

case. Also 7 = « i = sr7 • 

Zfiti Zh 

Nevertheless the physical equation 

dE = XdT + dx + JPd(o (I) 

remains unchanged in form with the altered value of x- 
Further when we substitute Clausius' expression (Art. 24) for 
Pdco, the equation 

remains unchanged in form with the altered values of h, %, 
and ^. 

Now let us consider the auxiliary function «. We might 
in the proof above given have put it in the form 



U = log jJ[(J:)' €-^ U+2(aH^Hy«)) ^X, . . . dz^du, ... dyx- 

In our present problem it would then take the form 



u 

^JJJ \'i 

A 



= logjj€~'^'x(^)' ^De'^'idsda, 



where ds stands for dxi,.,dz\ and da for dai.,.clyx. Then 
A' and o) are the two independent variables. Let 



A 

A\2 



f) \/56-*Q = JPl 



Then 9^ = 377 9A' + j— da>. 

ah ato 




X.J 
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Now -rn dh' 
ah 



-dh' 



'^'^Fdsda 



Ih 



-//.-. (i'^,(»-^)-^<W. 



-dhf 



-*'x Fdada 



//- 



Ih-H^^-S)'^''^ 



= - x^h 



.' + 77 ffe-'^m' ^ dsdi 



(1). 



Again -r- Oft) 



--9g> 



+ 



jj €-*'x Fdsda 



II 



e-*'x Fdsdar 






the last term referring to the variations of all the elements of 
length dx, dy, dz, consequent on the expansion denoted by 9®. 



That is -T- da) 
dm 

--9ft) 



'^'x Fdsda 



Ih 



II •-H'"t-'£)^ 



9ft) 



Oft) 



ft) 



+ 



9ft) 



//•- 



X Fdsdar 



II 



€-^'^ -J- dsdar 



■in 
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Adding together (1) and (2) we find 

+ 77 - 11'-"''' i^h' % + 3« ^ dsda- 

But now I / Fda = — is constant whatever be the values 

of h and to or h! and ©. It is therefore not affected by the 
variations of h' and 6>. 

These variations affect neither the quantities ai...7A, nor 
the limits of the da integration. 

We have then 

^^=xaiog2'+|3(A?^)+iatt, 

a result of precisely the same form as in the former case. 
136. In any case in which 

the proof of the second law for the two independent variables T 
and 6) can be presented in a simple form. 

For we have in that case from II 

Now let i = ^ , that is i is the time in which a particle 

moving with the velocity of mean square would describe a 
distance proportional to the linear dimensions of the system. 

Then |31ogft) = 2aiogi + 91ogT, 

and therefore 

= 2\a log {if). 
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We might define \2iT as the Action of the system par excel- 
lence. With that definition the Action is the same for all 
states of the system defined by dE=^ 0, that is for all states of 
the system on the same adiabatic line. And we may say that 

— = 9^ j either is, or is propor- 
tional to, the Action. 

136. But ^X~^^^ cannot exist for any system such as 

we have treated as a gas, except for a so-called perfect gas in 
which x = 0. If there be a potential x ** *11> ^^^ relative 
frequency of different configurations of the molecules must 
generally vary when T varies. And therefore although for any 
configuration the change of ^ as ««> varies is independent of jT, 
yet the mean of all values of x will not be independent of 
jT, because the variation of T causes some configurations, i.e. 
some values, to be represented with greater, and others with 
less frequency. We may for instance conceive a case in which 
the molecules placed at equal distances from each other would 
exert no sensible forces on each other by reason of their 
distances, and so x = 0. If we set them in motion they will 

in some cases approach so near to each other that x^^- 

137. Professor J. J. Thomson in the work above referred 
to gives a proof, not, as he says, of the second law of Thermo- 
dynamics, but of a proposition analogous thereto. He does not 
confine himself to the two variables T and c», but employs T 
with any number of controllable coordinates ^. Then he 
assumes that the mean potential energy Xy ^^ ^^ ^^ notation 
F, is completely fixed by the controllable coordinates, that is, 
is not a function of T. . That restriction prevents Professor 
Thomson's proof from being applicable to gases, as a proof that 

dE , , 

-^ is a function of temperature and volume, except in the 

limiting case of a perfect gas in which V is non-existent. Nor 
does Professor Thomson propose to apply it to such cases. 
It may well be that in the more general class of problems 
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treated of in his work the condition that V shall be a function 
of the controllable coordinates is not a serious restriction. 

138. The following is another form of solution of this 
problem. Let there be, instead of the one controllable coordi- 
nate 0), any number of them as coi ... a>n* We will now use w 
as a type controllable coordinate. It is assumed that the rate 
of change of any c» with the time is so slow that we may leave 
the velocities d)i ...(»„ out of account. Let Fi . . . Fn be the 
generalised components of momentum corresponding to a>i . . . ©»• 
Then it is a characteristic of the stationary motion that every 

dV 

-y. is zero. And therefore i( pi,,,pn be the components of 

external force required to maintain the motion stationary 

dT dxr u 

p = -^ 7^ for each c», 

^ da) da) 

by Lagrange's equations. See Art. 24. 

Generally let F be the function denoting frequency, that is 
^ is a function of the coordinates a?i ... ^n» and the velocities 
ai...7n, of the molecules of the system such that the time 
during which on average of any very long time, they lie between 
the limits 



Zfi . . . Zf^ -J- (^Zfi , 

tti . . . «! + dai , 



is Fdxi ... dzndai ... dr/n, or as we may write for brevity Fdsda, 
Then the mean value of any function as ;^ is 






day 



And dx = [[ (X^f'+^X) <^^^ I jj ^dsda. 
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In the same case we have, if t represent the kinetic energy 
in any particular state, T the mean value of t, 

dT=jj TdFdsda- +jJFdrd8d<r 
= jj TdFdsda +jJF^ da>d8da. 

139. Let us at this point assume the function ^ to be a 
product of two functions, F=ff\ where/ is a function of the 
coordinates only. In that case we have also 

or ^^~£ da>^jx^fds/jfds. 

And now 

dE=dT-\-dx+pdfo 

= 2ar + f x¥ds - [[ TdFdsda, 
because -;— 3 w = 1 1 -F t— dcodsda 



therefore 



da) 

= ar - /j TdFdsda, 



^-^^2dlogT + 2 jpfds-jj'^dFdsda. 

f)F 
And now we make -^ a complete differential by assuming 

where if> and yft are any functional symbols. 
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A solution of these equations is 

where C, C are the known constants. 

This solution is the one we should take under the assumption 
of condition A. Another solution is 

if Q = fiT. This is the solution corresponding to Chap. VIII. 

The subject of the Second Law of Thermodynamics has 
been very fiilly investigated by Professor G. H. Bryan in his 
Report on the present state of our knowledge of Thermo- 
dynamics, Part II., British Association, 1894 
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MATHEMATICAL PROPOSITIONS. 



In the preceding chapters I have assumed without proof 
certain propositions relating chiefly to determinants. They are 
for the most part very simple, but as some readers may think 
they require proof, I will devote this chapter to proving them. 
I shall require to make use of the following elementary 
property of determinants. If each of the constituents in one 
row is equal to the corresponding constituent in a second row 
multiplied by a constant factor, the determinant is zero. The 
two rows for which this is the case are said to be similar. 
Whatever general proposition is proved of rows is of course 
equally true of columns. 

(a) Let us take for a type determinant of n^ constituents 



i) = 



^11 > ^12 ••• ^in> 



^ni> ^n2 ••• ^1 



nn> 



in which Cpq is the constituent (whatever its value) in the pth 
row and qth column, and Oja is not necessarily equal to C^, &c. 
This can be made to represent any given determinant by 
assigning proper values to the constituents (7. 

« 

Let Dpq denote the minor of D formed by omitting the pth 
row and qth column, that is by striking out Cpq and the row 
and column in which it is ; Dqp the minor formed by omitting 
Cqp. In the determinant or any of its minors the diagonal 
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from the first constituent in the top row to the last in the 
bottom row shall be called the aoois, and its constituents aodcd 
constituents. 

(6) The sign of D^q shall be such that 

D = On Ai + (^iJ^u + . . . + GinDin 
= C/ja-Da + C/ja/^aa + • • . + G^D^ 

= &C. 

In order that this may be the case, we must in forming 
any minor Dpq apply the following rule of signs, namely to the 
product of all the axial constituents of that minor the positive 
or negative sign is to be attributed according as p H- 5^ is even 
or odd, and we deduce therefrom the sign of the product of any 
other set of constituents of the minor according to the usual 
rule in expanding determinants. In like manner Dpr^ shall 
denote the minor formed by omitting Gpq and G^g and the rows 
and columns in which they are. It follows from the above 
convention concerning signs that 

p=nq=n 

nD— X 2 GpqDpq. 
p=l q=l 

Then evidently -rp^ = Dpq, whatever values p and q may 

dLpq 

have. 

(c) Now consider the new determinant 

Gil — OiXi, Oia — 61X2 . . . Gin — OiXn, 
Ca — 0^1, C22 — 0^2 • • • G^n ~~ Ov^ny 



D' = 



If all the (7s in any two rows, as the pth and q\h, were zero, 
U would be zero, because the pth and gth rows would then be 
similar. It follows that every term in D contains at least 
n — \G factors, and therefore that U contains no products of 
any two or more of the new partial constituents 6x, For such 
products necessarily enter in pairs equal and of opposite sign, 
as for instance 6iX^^ — 61X16^, 

Therefore i)' = D - tXOpflc^pq. 
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(d) Now form a new determinant of (n + 1)* constituents 
by adding to Z) an (w + l)th row and (w + l)th column. The 
(n + l)th row shall be 

P(n+i) i> P(n+i) 2 • • • P(n+i) n> 2a,i^i, 

and the (n+ l)th column shall be symmetrical with it, contain- 
ing ^i(tt+i) in the first row, jSqin+D in the second, and so on. 
Also let I3pq = l3qpy for all values ofp and q. 

Let A be that new determinant, namely 

^11 > ^12 • • • ^in > Pi (n+i) > 
A= Cai, 02S*.. 02ni P2(n+i)> 



(a) Next in X) write 



/) ^ Piln+i)' z) «. _ Pi (n+i) P2 (n +i) 

^1^1 = "o > *'l^2 *= A > 



Then we have 






''n+1 



and by the ordinary expansion 

A = 2a„+iX) — 2,^Pp (n+i) fig (n+i) -Dp^. 



Therefore 



ly 



2a 



A. 



n+i 



(/) In this last deduction we have attributed no particular 
values to the C constituents. 

It would therefore hold for all the minors of 1/ not contain- 
ing the suflSx w + 1 as well as for jD' itself, namely 



and therefore 



1 — 


^^+1 


"11 » 


/ 
2 — 


1 

2a„+, 


Ai^; 


A.' 


Au 




jy 


A 


> 


A, 
ly 


~ A 


,&c. 



B. 



10 
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{g) Now consider the fimetion €~*^ in which h is any con- 
stant, and Q is a quadratic function of n variables u^yU^.^.Uny 
namely 

Q = Oi V + 612^1^2 + (h^^ + &c. 

Let us eflFect the integration 



/ 



€"*^ dun between limits + 00 . 



The terms in Q which contain u^ are 
which may be put in the form 



\Un^n 



+ 



2Va^ / 



4an 



Let us for a moment write 



^\n^ + fc2 n^2 + &C. _ 



Then 



and 



2Va„ 

J 00 

= €~*^ multiplied by a numerical factor. 
Comparing these results with (e), (/) we see that if 



D = 



Oi2f 2i(l^y 0^,..02ni 



D containing (/i-h 1)^ constituents, 

6 « 



and U = 



2a, - ^'"^ 



b^o — 



'12 



2an' 

2an 



'12 



1 PinP2n 



618 — 



6mfe 



in^'sn 



'18 



2an ' 



2a,- 



'2n 



2a„' • 



iy containing w' constituents ; 
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A/ Ai A/ A2 



then 



jy" D' ly" D' 



That is, these ratios are not altered by the change in the 
coefficients Oi, 61a, &c., consequent on the integration accord- 
ing to iin. Evidently the same law holds for every subsequent 
integration. Let us then perform the integrations 



/// 



and let the result be 



///■ 



If €"* Uit*i'+fiit«*i««+^«««') dAhdu^ is proportional to the chance 
that two variables u^ and v^ shall lie between the limits 
t*i . . . 1*1 + duy^ for the one, and iig . . . t^ + du^ for the other, we 
easily find the mean of u^, thus 

the limits of integration being for each u±qo 

" h 4iA^A^ - Aa' hd ' 
if d be the determinant 

2-4-1, Aa* 

Aaj ^ilg. 

But ^ = ^ by the last article. 

— Ai 

Therefore ^^ hD' 

A2 

Similarly UiU^ = ^ , 

— -^ 



and so on. 



10—2 
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(h) Under the same circumstances to find the mean value 
of -^ , or 2aiUi + butt^ + bi^v^ + &c., v^ being given (see Art. 97). 

That is to effect the integration 

III. ..du2 . . . dunc"^^ (2aiUi + bi^v^ + &c.). 



Let 



2aiUi + bi^Ut + ftuWg + &c. = P. 



The coeflScient of i^i in P is twice the coefficient of Mj* in Q. 
The coeflScient of every other ii, as i^, in P is equal to the 
coeflScient of Wjii, in Q. 

After integrating according to Un the coeflScient of Ui in P 

is 2ai — ^ while that of t*i* in Q is Oi — j~- . 

The coeflScient of t/j in P is 612 — ^-^, ai^d the coeflScient 

2an 

of i^it^ in Q is the same. 

The same relation between the coeflScients in P and those 
in Q which existed before integration according to Un exists 
also after that integration. The same rule holds for every 
subsequent integration to n^ inclusive. But when all the 
integrations have been eflfected to Wa inclusive, Q has become 

Wi*. Therefore P has become ^^ Wj which is the mean 



value sought. 



A. 



(t) Let us now consider the detenninant 



2) = 



Q2 



-M /i/s, A/i, (l+/s")>.- 



treated of in Chapter v, Art. 63. Here S is written for 2/ 
of that Article. 
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Comparing it with the type determinant, we find 

Ci2 = (?a = —fiS, &C., Ogg = Oja =/i/2i &C. 

Let it be required to find the values of D and its first 
minors. 

(j) We observe that in D every row except the top row con- 
tains one constituent of the form 1 -\-/^, e,g, the (p + l)th row 
contains 1 H-y^*. If it were /p^ instead of 1 +//, the {p + l)th 
row would be similar to the top row and therefore D would be 
zero. 

If we strike out any column other than the first, say the ^th 
column, we deprive the qth row of its constituent 1 •\-f^^. 
Hence the gth row is now similar to the top row. 

And if, to form a minor, we strike out any other row, neither 
the first nor the ^h, the determinant, i,e. the minor so formed, 
is zero. We see then that if /j > 1, and 5^ > 1, and p^q, Dpq = 0. 
But if we expand D, the constituent Cpq appears only in the 
form DpqCpq, and Dpg = 0, Therefore D when expanded does 
Dot contain any anaxial constituent except those in the first 
row or first column. 

From (i) and (j) we see that Z) = /S*. For the same reason 
B^=S^, and so on for every coaxial minor except D^, but as we 
shall see later Ai = S/* + 1. 

(A) We have next to find i)^, D^, &c. Striking out from 
B the top row and second column to form Aj, we find that in 
D12 every row except the top row has one constituent of the 
form 1 +/*. Therefore by the same reasoning as in (j) and by 
the rule as to the sign of the product of the axial constituents, 
-Oia = jDja = +/iS. When we form -Du, the product of the axial 
constituents is by our rule to be positive. But in order to 
reduce this minor to one in which the top row does not 
contain a constituent of the form 1 +/', we must interchange 
the two top rows, and thereby change the sign, of the minor. 
Therefore As = + /a^ aiid so on. In forming any minor Dip we 
have if p is odd, an odd number of changes of sign. 
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(0 Next let us deal with Ai* 



Let A = Ai = 






or 



All products of the form fp^f^ evidently disappear in 
expansion, and therefore all terms involving /* except in the 
first degree disappear, and therefore 



1 



n 



= S/', if/„=l. 

1 

t 

(m) Again, if we strike out any column except the first, say 
the second, we thereby deprive the second row of its constituent 
1 H-Za^ And therefore if we now strike out the top row we 
shall have 



A,,= 



42 



Jlj2» J3JS> J^J^t ••• 

Aft. (1+/.'), /./4. ... 
AU AA. (1+/A. 



= —fifi by the rule as to the sign of the axial constituents. 
Similarly Ai, = -/i/j. 

Also An = 2/»-/i'> 

A« = 2/^-/2S&c. 

(n) In the type determinant let now every C become 
C + 9(7, and let all the dC's be small, so that we may neglect 
products of them. That being the case 



dC 

P = l g = l dCpq 



da 



PQ' 
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Let US apply this to the determinant of (i) in which 

Cu = -s*, c„ = -M c„ — fs, &c., 

Here as we have seen -j^- = 0, unless either p = l, or q = l, 

dVpg 

OT p — q. 
Therefore 

dD = DadCj, +fiSdG^ +f^Ga + &c. 

+ S^C„ + S^G„ + &c. 
Now as in Art. 83 let us make 

aOu = 0, ac„ = ib,^, dc^ = it^s, &c. 

aC!ji = J6i^, dGa = ibin8, &c., 

Then as products of the b's are to be neglected, we have 
dD = + 6,„/.S» + b^f^' + b^f^" - &c. 
- Wi'S" - 6»./^' - &c. 
= 0. 
(o) Let us treat i)u, or A, in the same way. Here 

C'ii = l+/i', 0„=/y„ &c. 
And aAi = 3A, 

"^^aCu+^-ae^+fec. 



And a(7u = - 6in/i , 9^23 = - fean /a , &C., 

Therefore, since each anaxial term/i/a* &c. appears twice, 
a A = - (S/* -/i') 6i„/, - (S/^ -/.O h^f, ~ &c. 

-/i/a (^12 - ftm/j - &2n/i) -/i/s (^is - h.nfz - hnfi) - &C. 
~ "" ^2/ ly 2 "~ ^13 J if » — • • • ~- ^pqjpjq 

-^ib,nA + b^A + &c.) 

+/A«/x +/A„/x + . . . +/.'6a./. +/,'6^/, + &c. 
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That is aA = - fcu/i/a - 6u/i/» - &c. 

In this result every hpqfpfq is to be counted once only, that 
is, we do not distinguish 6^ from hqp. If hpqfpfq and bqpfqfp 
are counted separately, we must write 

and this form is the one required in Art. 83. 

(p) To prove that for certain values of 6, 2/j,26p^/<, = 2^*26 
see (33) of Art. 83. Let /(Art. 61) = -^^ when r>a. 

Taking for origin, let p be any point Op =p. About p 
describe a sphere of radius r less than p, and let j be a point 

on its surface. If w + 2= 1, or /= - , the mean value oi fq for 

all points q on the sphere is jr-y ovfp. And equation (33) is 
accurately satisfied, so long as r < p. 

If n + 2 > 1, we have in usual notation 

sin 0dd 



Jo 



W+2 

{p^-\-i^'' 2pr cos 6) 2 



1 (p-\- ry - (p - ry 



2npr (p* — r^y 



,^_ n.n— In— 2 ,. , . . « 

np^h- H 5—5 p»-»r» + &c. 

1 Z , o . 



npr (p* - r^Y 

pnrH\ 2.3 p2 A p^ p^ J 

=fp + terms containing/p multiplied by -^&c. 

And if now 6 = — f - j , where /a > 3, 

/p r 4>7rr^pbf^r = f^^ f ^irr^pbdr =//S6, 
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C 

because, - being negligible, all the terms in the integral which 
follow the term fp^^b are negligible compared with/p'^Sft. 
(q) To find the condition that the quadratic function (q) 

Q = a^Ui^ + buUjU^ + OaWa* + . . . 

shall be always positive whatever values be attributed to 

Evidently since Wj ... Un may be all zero, Ui must be positive. 
Similarly ag, as, &c., must be positive. 

Also, given Un, Q is minimum when 

JQ 

i— = 2anUn + bmUi + h^u^ + &c. = 0. 

And, when minimum, Q must be positive. Substitute in Q 
for Un the value found for it from 



that is 

Q then becomes 



2anUn + binVri + tan^a + &C. = 0, 
^in^ I ^2n^^ "T . . . 



Q^ = a/i^i' + ftiaX^ + a5'^2'* + &c., 



/) 3 h ^ 



where ai' = ai--j^, aa' = a2 — 7^ , &c.. 



and every a' must be positive. Also 612 = ^12 — ^-^ &c. 

Proceeding in the same way Q is ultimately reduced to 
Q=:-^Yr ^^ and as this must be positive, ^r must be positive. 

Similarly jr- must be positive, and so on. Therefore D and all 

its first coaxial minors have the same sign. But if we only 
reduce Q to a quadratic function of i^i and v^ 

Q-^ll -. 2 I -^12 ^. ^. ,1 -^82 « 

^-L/iaa -^laai -^71221 
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From which and the corresponding equations we see that 
all the second coaxial minors have the same sign as the first 
coaxial minors, and therefore the same sign as Z>. 

Similarly we can prove that the third coaxial minors must 
have the same sign as the second and so on, so that D and all 
its coaxial minors must have the same sign. And this sign 
must be positive because Oiyfh^ &c., are coaxial minors. This 
is the required condition. 

(r) Referring to Chapter v, in which it is proved that with 
(r) spheres of finite diameter c, 



where JIfr* or 



vz 






to prove the same thing for 



= \\\ dxdydk 



M^ = III dxdydza^ 



dec 



Consider the space between two infinite planes a; = and 
x=^dx, between which -^ is supposed constant and negative. 

Let ii = a — f . Then it is proved, as in Art. 64, that the number 
per unit of volume of molecules for which a — f lies between u 
and u + du, when a? = 0, exceeds the corresponding number when 
x=dxhy the quantity 

ax 

Therefore the number which enter the space between the planes 
exceeds the number which leave it per unit of time by 



C Ijdyd 



'ze-^'^^duih'ua^ . 

ax 



Hence as in Art. 66, it follows that, for each element of 
volume dxdydz, 
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ax 

Now u and a are not independent, as u and 7 were in 
Art. 66. Hence 

2A - 2A 3 3 



«*'«' = ^/ a* - 



2A' " 2K 4A» 2ft'2A ' 
and 2A'w«a3 = -^ . 

^^''''^ dxiA'dxH 2AWy ^ ^' 

which corresponds to 

\dz dxj dt 2 A \dz dx) ' 

of Art. 66. 

It follows, the motion being stationary, that 

didx 2h \dx) ' 
or by the process of Art. 68 

, 1 d/, , 3 /d^V /o\ 

-«*7:5;:(^ = «) = 2aU) (2)- 

When the material points are replaced by spheres of finite 

do^ 
diameter c, -^ acquires a new term due to collisions, which we 

will denote by ^ , and we have 

d^d^^d^djl 
dx dt dx dt * 

^xy y^y* y^zy being the components of the half relative 
velocity of two colliding spheres. 
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Also ^* = 27rc«/> F(\'» - \«) F», 

where >-F= F«, fiV=-Vy, i/F= F„ 

and using the same notation as in Arts. 70 — 71, we find 



dt 



= i 2ird'pV4h sf I ^*^ fdej'd<l>2 sin 5 cos ^ 

X (-\cos^ + Vl-\'sin d cos <^') |V(4cos* d- 4 cos'd) 
+ ( I - A.') sin'^ 20 cos" <^' - 2X Vl-V cos 2^ sin 2^ cos <^'} 

4 

whence introducing the factor Ce-^^'^^F'dF where 0=-j- (2A)* 

and integrating from F = to F = oo , 

aF/ 8 ^ 1 df 



or 






We have next to calculate the term 

due to collisions. That gives by the same method as that 
employed in Art. 76 

n 

ii^=-- 2ir<fp V j*d0J'd^2 sin ^ cos ^ 
X (- \ cos ^ + Vr^Oi' sin cos <f>') 
X (- 2\ COS* ff + Vl-V sin 20 cos <]>') -^ir= a). 

Whence we obtain for a pair of colliding spheres 
,9<if 4 ^ Tr. 44^,1 d/. , 
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and for each sphere separately 

.9 <^f A ^ T7i 22^,1 df, . 



=fl-^''40^y(2) w- 



Adding this to (3) we find for the change of Mja with the 
time due to collisions 

dM, 



30 . 1 /df\' 
Using Mgfg for Ma^g or M in Art. 79, we find 



dt dx 2h \dx) 



-^-4{0Hf)]-4(S)'-»- 

The first line is zero. The second can be made zero by 
making f;^) : i'T'] '-'^ k \ 1 — k b& before. 
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